随笔2：Stata 面板回归控制中的时间段问题【转】

Stata：levelof—区别变量的不同级别

levelof—区别变量的不同级别

语法格式：

levelsof varname [ if] [ in] [, options]

案例

levelof有两种不同的功能。首先，它提供了varname的不同值的列表。更常见的情况是，当你想要循环使用(比如)foreach;参见[P] foreach。levelof会留下一个r(levels)的列表，可以在后续的命令中使用。当需要获得非整数数据的级别时，可以使用matrow(matname)来获得完全精确的级别。

案例1：

. use https://www.stata-press.com/data/r17/auto

(1978 automobile data)

. levelsof rep78

1 2 3 4 5

. display "`r(levels)'"

1 2 3 4 5

. levelsof rep78, miss local(mylevs)

1 2 3 4 5 .

. display "`mylevs'"

1 2 3 4 5 .

. levelsof rep78, sep(,)

1,2,3,4,5

. display "`r(levels)'"

1,2,3,4,5

Showing value labels when defined:

. levelsof factor, local(levels)

. foreach l of locallevels {

. di "-> factor = `: label (factor) `l''"

. whatever iffactor == `l '

. }

案例2：回归合成控制法案例

. use growth, clear

. xtset region time

* Show the unit number of Hong Kong and treatment periods

. label list

. display tq(1997q3)

. display tq(2003q4)

* Replicate results inHsiao et al.(2012) with specified control units and designated post-treatment periods

. rcm gdp, trunit(9) trperiod(150) ctrlunit(4 10 12 13 14 19 20 22 23 25) postperiod(150/175)

如何得到计算得到上面的150/170/175

. xtset region time

panel variable: region (strongly balanced)

time variable: time, 1993q1 to 2008q1

delta: 1 quarter

. label list

region:

1 Australia

2 Austria

3 Canada

4 China

5 Denmark

6 Finland

7 France

8 Germany

9 HongKong

10 Indonesia

11 Italy

12 Japan

13 Korea

14 Malaysia

15 Mexico

16 Netherlands

17 NewZealand

18 Norway

19 Philippines

20 Singapore

21 Switzerland

22 Taiwan

23 Thailand

24 UnitedKingdom

25 UnitedStates

.

. . display tq(1997q3)

150

.

. . display tq(2003q4)

175

为什么是1997q3为150呢？因为time变量为132---192，而150对应的是1997q3，132对应为第一个季度数据

. levelsof time

132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 15

> 3 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 17

> 4 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192