听说你想知道st_data()函数的用法

下面我们通过例子来详细说明函数 st_data() 的用法。首先我们在Stata 中读入系统数据auto，键入如下命令：

sysuse auto,clear

keep in 1/10

接着进入Mata环境，先提取数值变量，输入命令如下：

mata

A1=st_data(1,2) //提取第二个变量的第一个观测值。

A1

A2=st_data(.,2) //提取第二个变量的所有观测值。

A2

A3=st_data((1\2\5),2) //提取第二个变量的第1、2、5个观测值。

A3

A4=st_data((1,5),2) // 提取第二个变量的第1-5个观测值。

A4

A5=st_data((1::5),2) // 同上

A5

A6=st_data((1,5\7,9),2) //提取第二个变量的第1-5和第7-9个观测值。

A6

A7=st_data(.,(3,9)) //提取第三、九个变量的所有观测值 。

A7

A8=st_data(.,("mpg","weight")) //提取变量mpg和weight的所有观测值。

A8

A9=st_data(., "mpg weight") //同上

A9

end

部分结果如下：

除此之外， st_data() 函数还可以直接提取变量的滞后变量和虚拟变量，举例如下：

先构造一个时间变量year，键入如下命令

en year=2000+_n

tsset year //定义变量year为时间变量

然后，进入Mata环境进行提取：

mata

A10=st_data(., ("price", "l.price")) // 提取变量price和变量price的滞后一期的所有观测值。

A10

A11=st_data(., "price l.price") // 同上

A11

A12=st_data(., "mpg i.rep78") //提取变量mpg和虚拟变量rep78的所有观测值，含缺失值。

A12

部分结果如下：

接下来，我们看看例子中加入选项 selectvar 的含义，键入如下命令：

tab rep78, generate(r) //生成变量rep78的虚拟变量

mata

A13=st_data(., "mpg weight")

A13

A14=st_data(., "mpg weight","r2") //虚拟变量r2取值是非0时,提取 变量 mpg 和weight对应的观测值

A14

A15=st_data(., "mpg weight",15) //同上，15表示变量r2在stata内存中变量的序号

A15

A16=st_data(., "mpg i.rep78", ) // 提取变量mpg和虚拟变量rep78的所有   观测值，含缺失值。

A16

A17=st_data(., "mpg i.rep78",.) // 同上

A17

A18=st_data(., "mpg i.rep78",0) //提取变量mpg和虚拟变量rep78非缺失的观测值

A18

end

部分结果如下：