Stata:因子变量全攻略

发布时间:2020-06-04 阅读 373

作者:连玉君 | 杨柳

Email: arlionn@163.com

注:该文已发表: 连玉君, 杨柳.《郑州航空工业管理学院学报》, 2018, 36(2): 90-103. 【 -点击下载-】 || 本文PDF版本2

Stata连享会   主页 || 视频 || 推文

连享会 - Stata 暑期班

线上直播 9 天:2020.7.28-8.7
主讲嘉宾:连玉君 (中山大学) | 江艇 (中国人民大学)
课程主页https://gitee.com/arlionn/PX | 微信版


目录


事倍功半 vs. 事半功倍

当需要控制公司个体效应、面对成千上万家公司的数据资料时,你会如何处理?或许你会首先想到对每家公司生成虚拟变量,但是这种做法的工作量实在是太大了!那么,有没有事半功倍的方法呢?下面,小编就带你学习Stata软件中因子变量的使用方法。

1. 问题背景

实证分析中,我们经常需要在模型中加入反映类别的虚拟变量,以便控制不可观测的组间差异。而在另一些分析中,为了刻画调节效应,尚需在模型中加入变量的交乘项平方项。传统的做法是,预先生成虚拟变量或交乘项,进而将它们加入模型。然而,当虚拟变量或交乘项的数目较多时,上述方法就显得尤为不便,不但浪费计算机内存,也会严重降低我们的工作效率。

在Stata中,我们可以使用因子变量(Factor Variable)简化操作步骤、快捷地在回归模型中加入虚拟变量交乘项平方项高次项。更为重要的是,由于引入交乘项或平方项后,解释变量对被解释变量的边际影响不再是常数,而是某个变量(调节变量)的函数,在有些模型设定下,这种关系可能是非线性的。此时,若使用因子变量,并配合Stata中的marginsmarginsplot命令,可以非常便捷、直观地分析关键变量的边际效应

2. 什么是因子变量

因子变量(Factor Variable)是对现有变量的延伸,是从类别变量中生成虚拟变量、设定类别变量之间的交乘项、类别变量与连续型变量之间的交乘项或连续变量之间的交乘项(或多项式)。在Stata中的大多数回归命令和回归后的估计命令中都可以使用这些因子变量^[1]^。

因子变量的五种运算符及其含义如下表1所示:

表1 因子变量的运算符及含义 表1 因子变量的运算符及含义.png

以研究妇女工资的决定因素为例,使用Stata软件自带的数据文件 nlsw88.dta。该数据包含了1988年采集的2246个美国妇女的资料,包括:小时工资wage、每周工作时数hours、种族race、职业occupation、年龄age、是否大学毕业collgrad、当前职业的工作年限tenure、是否结婚married、是否居住在南部地区south、合计工作年限ttl_exp等变量。其中,小时工资wage、每周工作时数hours、年龄age、当前职业的工作年限tenure、合计工作年限ttl_exp为连续型变量;种族race为类别变量(1代表白种人white,2代表黑种人black,3代表其他人种other)、职业occupation为类别变量(13个职业类别);是否大学毕业collgrad、是否结婚married、是否居住在南部地区south为虚拟变量。

在这份数据中有一个表示种族的类别变量race,取值为1、2、3,分别对应白人黑人其他人种。假设我们想在模型中加入一个反映种族的虚拟变量black,当某个妇女是黑人时,black取值为1,否则为0。则传统的做法如下^[3]^:

  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
. sysuse "nlsw88", clear. gen black=1. replace black=0 if race!=2. reg wage black

若延续这一思路,但使用因子变量来生成black变量,则命令为:

  • ounter(line
. gen black = 2.race

只需要一条命令,而且含义非常明确。这里2.race本质上是一个条件判断语句:判断某一行观察值中的race变量取值是否为2,若是,则返回1到变量black中,否则返回0

然而,在多数情况下,我们的目的只是希望得到虚拟变量black的估计系数,而不希望生成或存储这个变量^[4]^。Stata中的因子变量语法完全注意到了这个问题。使用因子变量的标准做法如下:

  • ounter(line
  • ounter(line
. sysuse "nlsw88.dta", clear. reg wage 2.race

注意,我们无需预先生成black变量,而是直接在回归模型中加入了2.race因子变量。有些读者注意到,race变量有三个取值,因此,我们可以在模型中放入两个虚拟变量,此时可以书写如下命令:

  • ounter(line
. reg wage i.race

回归结果如下所示:

  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
(部分回归结果省略)------------------------------------------------------------------------------        wage |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]-------------+----------------------------------------------------------------        race |      black  |     -1.238      0.276    -4.48   0.000       -1.781      -0.696      other  |      0.468      1.133     0.41   0.680       -1.754       2.690             |       _cons |      8.083      0.142    57.06   0.000        7.805       8.361------------------------------------------------------------------------------

从上述结果中可以看到,回归模型中加入了black虚拟变量和other虚拟变量,分别对应race变量的第二个和第三个类别,而第一个类别white被Stata默认作为基准组,目的在于防止完全共线性。black变量的系数值为-1.238,表示黑人的平均工资比白人低1.238个单位,并在统计上显著;other变量的系数值为0.468,表示其他人种的平均工资比白人高0.468个单位,但是在统计上不显著。

在实证分析中,有时会根据研究内容的需要改变基准组的设定,此时可使用ib.b.的前缀,具体写法如下表2中所描述:

表2 基准组的运算符及含义 表2 基准组的运算符及含义.png

在实证分析中,变量的交乘项高次项往往是重要的解释变量。以研究妇女工资的决定因素为例,若我们想在模型中加入黑人每周工作时数的变量,则传统做法是预先生成一个虚拟变量 black 代表是否黑人,再生成一个新变量 black_x_hours 表示黑人与每周工作时数的交乘项,然后再将这个新变量放入回归模型中。Stata命令为:

  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
. gen black=1. replace black=0 if race!=2. gen black_x_hours = black * hours. reg wage black_x_hours

若在Stata中使用因子变量实现上述过程,则命令十分简洁:

  • ounter(line
. reg wage 2.race#c.hours

下面,我们介绍如何在Stata中使用因子变量表示变量的交乘项高次项。我们以研究妇女工资的决定因素为例进行说明。

(1) 两个类别变量的交乘项: 在回归模型中加入种族race和职业类别occupation交乘项,Stata命令为:

  • ounter(line
. reg wage i.race#i.occupation

若在回归模型中既要放入种族race和职业类别occupation虚拟变量,又需要同时放入这两个变量的交乘项,则在回归命令中使用i.race##i.occupation,相应的命令为:

  • ounter(line
. reg wage i.race##i.occupation

(2) 类别变量与连续变量的交乘项: 在回归模型中加入种族race和每周工作时数hours的交乘项,Stata命令为:

  • ounter(line
. reg wage i.race#c.hours

需要注意的是,在上例中,由于我们把hours变量视为连续变量,因此,需要在其前面加上c.符号以便告知Stata该变量是连续变量

(3) 连续变量与连续变量的交乘项(高次项): 在回归模型中加入年龄age变量,以及其平方项,Stata命令如下:

  • ounter(line
. reg wage c.age##c.age

上述命令中c.age表示年龄age变量被当成连续型变量。如果我们在Stata命令中使用i.age,则年龄age变量被当成类别变量处理,此时,类别的个数为年龄age变量中不同取值的个数。

3. 常用回归模型的因子变量表述

3.1 范例 1:邹氏检验

由于不同组别之间可能会存在差异(截距项或斜率项存在差异),因此,我们需要检验这些差异在统计上是否显著。这时,我们可以使用邹氏检验。以研究妇女工资的影响因素为例,我们可以使用chowtest命令检验工会成员与非工会成员两个样本组中工资影响因素是否存在差异(或称之为存在结构变化),Stata命令如下:

  • ounter(line
  • ounter(line
. global xx "hours age tenure ttl_exp married". chowtest wage $xx, group(union) detail

事实上,我们也可以使用因子变量的语法,在回归模型中加入分组变量其他控制变量交乘项,然后再联合检验分组变量的系数以及所有交乘项的系数是否都等于0。此时,即使我们不使用chowtest命令,也可以轻松实现邹氏检验,Stata命令如下:

  • ounter(line
  • ounter(line
  • ounter(line
. global xx "hours age tenure ttl_exp married". reg wage $xx i.union i.union#c.($xx). testparm i.union i.union#c.($xx)

连享会 - 文本分析与爬虫 - 专题视频

主讲嘉宾:司继春 || 游万海

连享会-文本分析与爬虫-专题视频教程
连享会-文本分析与爬虫-专题视频教程

3.2 范例 2:双向固定效应模型

在实证分析过程中,经常需要在模型中加入反映年度、公司或行业特征虚拟变量。当虚拟变量的数目众多时,采用手动输入变量的方式会非常耗时。例如,下述模型(1)是文献中广泛应用的双向固定效应模型

上式等价于

其中,Di 是对应于每家公司的虚拟变量,对于公司 iDi=1,否则 Di=0。对于一个有 N=1000 家公司的面板数据而言,共有 1000 个虚拟变量。Wt 是 T1 个年度虚拟变量^[5]^, 定义方式与 Di 相似。

若使用因子变量,则语法很简单^[6]^:

  • ounter(line
. reg y x1 x2 x3 i.id i.year

其中,i.id 表示 N-1 个公司虚拟变量(Stata 会默认将一家公司作为基准组),即(1)式中的 αi

当然,我们也可以用 xtreg 命令自带的 fe 选项来控制个体效应 ,同时使用 因子变量 来加入年度虚拟变量 ,命令如下:

  • ounter(line
. xtreg y x1 x2 x3 i.year, fe

3.3 范例 3:DID 模型

在实证分析中,常常需要分析政策实施后带来的效果。这时,我们就需要采集实验组控制组两组样本(实验组的样本代表政策实施前后的情况,控制组的样本代表不实施政策的情况),再把这两组样本合并为一份数据后进行回归分析。在建模时,我们需要设定一个是否为实验组或控制组的处理虚拟变量 Treat 放入模型,还需要设定一个反映政策实施时间前与实施后的时间虚拟变量 Time 放入模型,此外,还需要在模型中加入 Treat 变量与 Time 变量的交乘项,该交乘项的系数估计值 a3 就是政策实施后带来的效果,即实验组在政策实施后与假设未实施政策的差异。回归模型如下式(3)所示:

在 Stata 中的命令写法如下:

  • ounter(line
  • ounter(line
  • ounter(line
. reg y Treat Time Treat#Time x1 x2 x3*-或者写为. reg y Treat##Time x1 x2 x3

在多期 DID 分析中,我们常常需要加入年度虚拟变量处理变量交乘项来检验共同趋势假设(common trend)以及政策效果。例如,在 Acemoglu and Angrist(2001) 文中,作者采集了1988年至1997年的人口调查数据,将样本分为残疾人与非残疾人,使用多期DID模型研究了1992年美国对残疾人工作保护法案(ADA)的实施效果。文中模型设定如下(pp.925):

其中,i 为残疾人个体;t 为年份;yit 为工作周数或周平均工资;xit 为一系列控制变量,pit 为时间趋势项;dit 为是否残疾,其系数为 δ; αt 为随年份变化的实施 ADA 与未实施 ADA 对残疾人工作保护效应的系数。当 t1992 年时,αt 量化了实施ADA之后对残疾人工作保护的效果(使用非残疾人作为控制组);当 t<1992 年时,αt 检验了实施 ADA 之前的各年份中残疾人和非残疾人的工作周数或周平均工资是否在统计上具有显著差异,相当于同时进行了DID模型的 “共同趋势假设” 的检验,若 αt 系数在1992年之前不显著,表明存在**“共同趋势”**。Stata 命令如下^[7]^:

  • ounter(line
  • ounter(line
  • ounter(line
. use "ABA_JPE2001.dta", clear. global controls "i.age_G i.edu_G i.race_G i.region". reg wkswork1 i.year##($controls) i.disabled##i.year

Stata回归结果如下:

  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
(部分回归结果省略)-----------------------------------------------------------------------------------         wkswork1 |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]------------------+----------------------------------------------------------------------------------|----------------------部分回归结果省略----------------------------                  |    disabled#year |     Disabl#1988  |  -.6873656   .7197987    -0.95   0.340    -2.098154    .7234227     Disabl#1989  |  -.5385104   .7056609    -0.76   0.445    -1.921589    .8445682     Disabl#1990  |  -2.289956   .6945586    -3.30   0.001    -3.651274   -.9286372     Disabl#1991  |  -2.158817   .6916205    -3.12   0.002    -3.514376    -.803257     Disabl#1992  |  -1.386539   .6873483    -2.02   0.044    -2.733725   -.0393526     Disabl#1993  |  -2.743325   .6953003    -3.95   0.000    -4.106097   -1.380553     Disabl#1994  |  -3.917873   .7198205    -5.44   0.000    -5.328704   -2.507042     Disabl#1995  |  -3.702671    .748066    -4.95   0.000    -5.168863    -2.23648     Disabl#1996  |  -4.472484   .7373381    -6.07   0.000    -5.917649   -3.027319                  |            _cons |   41.05231   .4546278    90.30   0.000     40.16125    41.94337-----------------------------------------------------------------------------------

连享会 - 效率分析专题

已上线:可随时购买学习+全套课件,课程主页已经放置板书和 FAQs
主讲嘉宾:连玉君 | 鲁晓东 | 张宁
课程主页微信版PDF版

连享会-效率分析专题视频
连享会-效率分析专题视频

3.4 范例 4:超越对数生产函数

很多文献使用超越对数生产函数估计生产效率,模型中会包含投入要素 KL 的高阶项,包括平方项和二者的交乘项^[8]^。例如在 Kumbhakar(1989) 中,作者通过在模型中加入二次项来捕捉要素的交互影响和潜在的非线性关系,模型如下式(5)所示。

此时,使用因子变量会让Stata中的命令变得异常简洁,如下所示:

  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
. webuse frontier1.dta, clear  . global y "lnoutput"            // y. global x "lnlaborlncapital"    // x1 x2*-超越对数生产函数. sfcross $y c.($x)##c.($x)      // Eq.(5)

4. 边际效应分析和图形化呈现

一些实证研究的文献中常常加入变量的交乘项以反映可能存在的调节效应,例如 Faulkender and Wang(2006)戴魁早和刘友金(2016)张苏和高扬(2012) ^[9]^,此时最重要的是分析 xy边际效应

x 对 y 的边际影响为:

显然,xy 的边际影响取决于 z 的取值。这种边际效应分析若是采用手动计算会非常繁琐,而后续的图形呈现则更为复杂。然而,若是借助因子变量,并配合Stata中的marginsmarginsplot命令,对于包含交乘项的模型中的边际效应的分析和图形化展示都变得异常轻松。

下面,我们以研究妇女工资的决定因素为例进行说明。使用Stata软件的自带数据nlsw88.dta(1988年美国妇女小时工资),以wage(妇女的小时工资)作为被解释变量,以race(种族类别)collgrad(是否大学毕业)race与collgrad的交乘项作为解释变量建立线性回归模型,Stata中的命令如下:

  • ounter(line
  • ounter(line
. sysuse "nlsw88.dta", clear. reg wage i.race##collgrad

回归结果如下所示:

  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
(部分回归结果省略)-------------------------------------------------------------------------------------               wage |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]--------------------+----------------------------------------------------------------               race |             black  |     -1.442      0.297    -4.85   0.000       -2.026      -0.859             other  |     -0.380      1.348    -0.28   0.778       -3.024       2.263                    |           collgrad |      college grad  |      2.981      0.312     9.54   0.000        2.368       3.593                    |      race#collgrad |black#college grad  |      2.502      0.676     3.70   0.000        1.177       3.827other#college grad  |      1.678      2.297     0.73   0.465       -2.826       6.182                    |              _cons |      7.318      0.158    46.25   0.000        7.008       7.629-------------------------------------------------------------------------------------

我们使用margins命令来计算种族(race)与是否大学毕业(collgrad)的交乘项的各个类别妇女小时工资(wage)预测边际值,Stata命令和结果如下所示:

  • ounter(line
. margins i.race#collgrad

计算结果如下所示:

  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
Adjusted predictions                            Number of obs     =      2,246Model VCE    : OLSExpression   : Linear prediction, predict()-----------------------------------------------------------------------------------------                        |            Delta-method                        |     Margin   Std. Err.      t    P>|t|     [95% Conf. Interval]------------------------+----------------------------------------------------------------          race#collgrad |white#not college grad  |   7.318251   .1582176    46.25   0.000     7.007983     7.62852    white#college grad  |   10.29895   .2693243    38.24   0.000     9.770797     10.8271black#not college grad  |   5.875918   .2519298    23.32   0.000     5.381878    6.369958    black#college grad  |   11.35861    .543853    20.89   0.000      10.2921    12.42512other#not college grad  |   6.938195   1.338677     5.18   0.000     4.313019    9.563371    other#college grad  |   11.59678   1.839835     6.30   0.000     7.988818    15.20474-----------------------------------------------------------------------------------------

使用marginsplot命令将计算结果用图的形式表示,Stata命令和结果如下所示:

  • ounter(line
. marginsplot

图1 种族与是否大学毕业交乘项各类别的妇女小时工资的预测边际值 图1 种族与是否大学毕业交乘项各类别的妇女小时工资的预测边际值.png

从图1中可以直观的看到,不论种族类别,大学毕业的妇女的平均工资高于非大学毕业的妇女,该结果符合我们的一般认知;另外,我们从图中还发现一个有趣的结果,非大学毕业的白人妇女的平均工资高于非大学毕业的黑人妇女,而大学毕业的白人妇女的平均工资低于大学毕业的黑人妇女。

下面,我们使用margins命令附加atmeans选项来计算当其他变量取均值不同种族类别的妇女小时工资的预测边际值,Stata命令和结果如下所示:

  • ounter(line
  • ounter(line
  • ounter(line
. margins i.race, atmeans*-或者写为. margins race, atmeans

计算结果如下所示:

  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
Adjusted predictions                            Number of obs     =      2,246Model VCE    : OLS
Expression : Linear prediction, predict()at : 1.race = .7288513 (mean) 2.race = .2595726 (mean) 3.race = .0115761 (mean) 0.collgrad = .7631345 (mean) 1.collgrad = .2368655 (mean)------------------------------------------------------------------------------ | Delta-method | Margin Std. Err. t P>|t| [95% Conf. Interval]-------------+---------------------------------------------------------------- race | white | 8.024276 .136558 58.76 0.000 7.756482 8.292069 black | 7.174578 .231424 31.00 0.000 6.72075 7.628406 other | 8.041653 1.110659 7.24 0.000 5.863625 10.21968------------------------------------------------------------------------------

使用marginsplot命令将计算结果用图的形式表示,Stata命令和结果如下所示:

  • ounter(line
. marginsplot

图2 不同种族类别的妇女小时工资的预测边际值 图2 不同种族类别的妇女小时工资的预测边际值.png

从图2中可以直观的看到,白人与其他种族的妇女的平均工资高于黑人妇女。

两个连续变量交乘项对被解释变量的边际效应也可以使用margins命令来计算。我们仍以研究妇女工资的决定因素为例进行说明。在回归模型中加入tenure(当前职业的工作年限)及其平方项,并将hours(每周工作时数)age(妇女年龄)married(是否结婚)south(是否居住在南部地区)race(种族类别)作为控制变量,Stata命令如下所示:

  • ounter(line
  • ounter(line
. global xx "hours age married south i.race". reg wage c.tenure##c.tenure $xx

回归结果如下:

  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
(部分回归结果省略)-----------------------------------------------------------------------------------             wage |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]------------------+----------------------------------------------------------------           tenure |   .3015859   .0704148     4.28   0.000        .1635    .4396718                  |c.tenure#c.tenure |  -.0076366   .0037758    -2.02   0.043    -.0150411   -.0002321                  |            hours |   .0749604   .0115722     6.48   0.000     .0522669    .0976539              age |  -.0850009   .0388991    -2.19   0.029    -.1612832   -.0087185          married |  -.4619864   .2544385    -1.82   0.070     -.960949    .0369761            south |  -1.283751   .2493374    -5.15   0.000     -1.77271   -.7947914                  |             race |           black  |  -1.221256   .2870854    -4.25   0.000    -1.784241   -.6582719           other  |   .3129332   1.098312     0.28   0.776    -1.840893     2.46676                  |            _cons |   8.183759   1.632092     5.01   0.000     4.983171    11.38435-----------------------------------------------------------------------------------

由于在模型中加入了 tenure 的平方项,因此,tenurewage 的边际效应将会受到 tenure 取值的影响。上述Stata命令对应的模型设定如下:

tenurewage 的边际影响为:

显然,当 tenure 取值不同时,tenurewage 的边际效应是不相同的。因此,需要先使用描述性统计分析的命令(如sum tenure)查看 tenure 取值的范围,然后再计算当 tenure 取不同的值所对应的边际效应,Stata命令如下:

  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
. preserve. keep if e(sample). sum tenure. restore. margins, dydx(tenure) at(tenure=(0 1(3)25 25.9))

计算结果如下所示:

  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
Average marginal effects                        Number of obs     =      2,227Model VCE    : OLS
Expression : Linear prediction, predict()dy/dx w.r.t. : tenure1._at : tenure = 02._at : tenure = 1... (output omitted) ...11._at : tenure = 25.9------------------------------------------------------------------------------ | Delta-method | dy/dx Std. Err. t P>|t| [95% Conf. Interval]-------------+----------------------------------------------------------------tenure | _at | 1 | .3015859 .0704148 4.28 0.000 .1635 .4396718 2 | .2863127 .0632764 4.52 0.000 .1622256 .4103998... (output omitted) ... 11 | -.0939893 .1304534 -0.72 0.471 -.3498129 .1618343------------------------------------------------------------------------------

使用marginsplot命令将计算结果用图的形式表示,Stata命令和结果如下所示:

  • ounter(line
. marginsplot, xlabel(,format(%3.1f) angle(60))

图3 工作年限对妇女工资的平均边际效应 图3 工作年限对妇女工资的平均边际效应.png

Stata中的多数命令都支持marginsmarginsplot命令。因此,即使对于非线性模型,如logit, tobit等,我们仍然可以借助这两个命令很方便地分析边际效应。

5. 输出回归结果时的问题及解决办法

以研究妇女工资的决定因素为例。wage(妇女的小时工资)作为回归模型的被解释变量,race(种族类别)collgrad(是否大学毕业)race与collgrad的交乘项作为解释变量,并将回归结果输出到Excel中,Stata命令和结果如下:

  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
. sysuse "nlsw88.dta", clear. reg wage i.race##collgrad. est store R1. esttab R1 using D:/Table_factor_1.csv, nogap replace

表3 1988年美国妇女工资模型结果 表3 1988年美国妇女工资模型结果.png

我们发现在表3中有很多变量的系数值为 0,并缺失 t 统计量。造成这种情况的原因有二:一方面是由于有些虚拟变量作为了基准组,例如:1.race0.collgrad1.race#0.collgrad,Stata默认将它们作为基准组,所以就缺失这些基准组的估计系数值和标准误;另一方面是由于有些交乘项中的其中一个因子变量是基准组,其它变量与这个作为基准组的因子变量交乘后的交乘项就被忽略了,所以其估计系数值和标准误就会缺失,例如:1.race#1.collgrad2.race#0.collgrad3.race#0.collgrad。此时,可以使用esttab命令的drop()选项来屏蔽这些系数的显示,还可以使用nobasenoomit的选项,Stata命令和结果如下所示:

  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
  • ounter(line
*-输出结果(不显示基准组和忽略组的系数,使用drop选项). esttab R1 using D:/Table_factor_2.csv, nogap      ///	     drop(1.race 0.collgrad 1.race#0.collgrad   ///         1.race#1.collgrad 2.race#0.collgrad        ///         3.race#0.collgrad) replace*-输出结果(不显示基准组和忽略组的系数,使用nobase与noomit选项). esttab R1 using D:/Table_factor_3.csv, nogap nobase noomit replace

表4 1988年美国妇女工资模型结果(删除缺失估计系数值的变量) 表4 1988年美国妇女工资模型结果(删除缺失估计系数值的变量).png

6. 小结

本文介绍了Stata中因子变量产生虚拟变量交乘项的使用方法,以常用经典回归模型为例,提供了它们的Stata命令,并进一步提供了因子变量marginsmarginsplot命令相配合分析边际效应的示例。在Stata中的大多数命令中都可以使用因子变量的表述方法。该方法可以使Stata命令更加简洁,并能够大幅度提高实证分析的效率,但需要注意分析使用因子变量表述方法后得到的模型设定结构

参考文献:

Acemoglu, D., J. D. Angrist, 2001, Consequences of Employment Protection? The Case of the Americans with Disabilities Act, Journal of Political Economy, 109 (5): 915-957.
Altunbas, Y., M.-H. Liu, P. Molyneux, R. Seth, 2000, Efficiency and Risk in Japanese Banking, Journal of Banking & Finance, 24 (10): 1605-1628.
Faulkender, M., R. Wang, 2006, Corporate Financial Policy and the Value of Cash, Journal of Finance, 61 (4): 1957-1990.
Kumbhakar, S. C., 1989, Estimation of Technical Efficiency Using Flexible Functional Form and Panel Data, Journal of Business and Economic Statistics, 7 (2): 253-258.
Wang, E. C., 2007, R&D Efficiency and Economic Performance: A Cross-Country Analysis Using the Stochastic Frontier Approach, Journal of Policy Modeling, 29 (2): 345-360.
戴魁早, 刘友金, 2016, 要素市场扭曲与创新效率——对中国高技术产业发展的经验分析, 经济研究, (7): 72-86.
王德祥, 李建军, 2009, 我国税收征管效率及其影响因素——基于随机前沿分析(SFA)技术的实证研究, 数量经济技术经济研究, (4): 152-160.
张苏, 高扬, 2012, 大学生学习行为与国家竞争力关联关系的实证研究, 管理世界, (4): 175-176.


注释: [1] 详情参阅Stata帮助文件help fvvarlist
[2] 详情参见Stata帮助文件help varlist
[3] 为了便于说明,后续多数回归命令中都省略了控制变量。
[4] 事实上,只要你的数据中存储了race变量,我们只需要保存好dofile文件,就无需生成black这个中间变量。
[5] 之所以加入T - 1 个年度虚拟变量,是为了防止完全共线性。
[6] 此处设定i.year会自动加入 T 个年度虚拟变量,但Stata会自动删除一个,以防完全共线性。
[7] 文中实证分析所用的原始数据和相关程序可以从作者主页上下载:http://economics.mit.edu/faculty/acemoglu/data/aa2001。 [8] 例如,Altunbas, Liu, Molyneux and Seth(2000) 使用超越对数成本函数估算了日本银行的效率和风险;Wang(2007) 则使用超越对数生产函数研究30个国家R&D效率;王德祥和李建军(2009) 基于超越对数生产函数估算了我国的税收流失率。
[9] Faulkender and Wang(2006) 检验了由公司融资约束对现金持有边际市场价值的影响。戴魁早和刘友金(2016) 研究发现了要素市场扭曲对创新效率的影响存在着企业差异,企业规模在规避要素市场扭曲对创新效率的抑制效应中具有积极作用。张苏和高扬(2012) 的实证研究发现学生来源于城市或农村地区对国家竞争力的影响作用受到每周上网时间的影响,如果上网时间在每周8小时以下,城市大学生的学习行为落入“增进国家竞争力导向”上高效率区域而不是低效率区域的概率比农村大学生要高,否则要低。

相关课程

连享会-直播课 上线了!
http://lianxh.duanshu.com

免费公开课:


课程一览

支持回看,所有课程可以随时购买观看。

专题 嘉宾 直播/回看视频
Stata暑期班 连玉君
江艇
线上直播 9 天
2020.7.28-8.7
效率分析-专题 连玉君
鲁晓东
张 宁
视频-TFP-SFA-DEA
已上线,3天
文本分析/爬虫 游万海
司继春
视频-文本分析与爬虫
已上线,4天
空间计量系列 范巧 空间全局模型, 空间权重矩阵
空间动态面板, 空间DID
研究设计 连玉君 我的特斯拉-实证研究设计-幻灯片-
面板模型 连玉君 动态面板模型-幻灯片-
直击面板数据模型 [免费公开课,2小时]

Note: 部分课程的资料,PPT 等可以前往 连享会-直播课 主页查看,下载。


关于我们

  • Stata连享会 由中山大学连玉君老师团队创办,定期分享实证分析经验。直播间 有很多视频课程,可以随时观看。
  • 连享会-主页知乎专栏,300+ 推文,实证分析不再抓狂。
  • 公众号推文分类: 计量专题 | 分类推文 | 资源工具。推文分成 内生性 | 空间计量 | 时序面板 | 结果输出 | 交乘调节 五类,主流方法介绍一目了然:DID, RDD, IV, GMM, FE, Probit 等。
  • 公众号关键词搜索/回复 功能已经上线。大家可以在公众号左下角点击键盘图标,输入简要关键词,以便快速呈现历史推文,获取工具软件和数据下载。常见关键词:
    • 课程, 直播, 视频, 客服, 模型设定, 研究设计, 暑期班
    • stata, plus,Profile, 手册, SJ, 外部命令, profile, mata, 绘图, 编程, 数据, 可视化
    • DID,RDD, PSM,IV,DID, DDD, 合成控制法,内生性, 事件研究, 交乘, 平方项, 缺失值, 离群值, 缩尾, R2, 乱码, 结果
    • Probit, Logit, tobit, MLE, GMM, DEA, Bootstrap, bs, MC, TFP, 面板, 直击面板数据, 动态面板, VAR, 生存分析, 分位数
    • 空间, 空间计量, 连老师, 直播, 爬虫, 文本, 正则, python
    • Markdown, Markdown幻灯片, marp, 工具, 软件, Sai2, gInk, Annotator, 手写批注, 盈余管理, 特斯拉, 甲壳虫, 论文重现, 易懂教程, 码云, 教程, 知乎

连享会主页  lianxh.cn
连享会主页 lianxh.cn

连享会小程序:扫一扫,看推文,看视频……


扫码加入连享会微信群,提问交流更方便

✏ 连享会学习群-常见问题解答汇总:
https://gitee.com/arlionn/WD