如何interchange数据中一部分的obs和vars

数据是从NBER下的一个有关patent所属公司和其他信息的dta，大体结构如下

pdpass        pdpco1        source        begyr1        gvkey1        endyr1        pdpco2        begyr2        gvkey2        endyr2
32040490        65907        m2006        1995        65907        2001        135484        2002        135484        2005

pdpass最多有gvkey5，大部分只有一两个，后面就都missing了

我想要得到的结果大概是这样

pdpass         year           gvkey
32040490     1995          65907
32040490     1996          65907
32040490     1997          65907
............................
32040490     2001          65907
32040490     2002          135484
32040490     2003          135484
32040490     2004          135484
32040490     2005          135484

求高人指点，应该如何写命令，附上原始数据（后缀应为dta，但上传失败，改为txt）。注：有些pdpass的begyr和endyr是不连续的，比如endyr1 是1998， begyr2可能是2003.

dynass.txt (910.98 KB)

2012-3-7 16:53:27 上传

use dynass, clear
reshape long pdpco begyr endyr gvkey, i(pdpass source)
keep if gvkey<.
expand endyr-begyr+1
bys pdpass gvkey: g year=begyr+_n-1
keep pdpass year gvkey

没看太懂，gvkey5和gvkey是什么关系？

aolei 发表于 2012-3-7 17:22
没看太懂，gvkey5和gvkey是什么关系？

同一个公司不同时期用的gvkey，分别称作gvkey1,2,3,...

做出来了，只是方法比较笨，要手动重复五次：
use dynass.dta
keep pdpass-endyr1
gen ny=endyr1-begyr1+1
expand ny
sort pdpass
bysort pdpass: gen year=begyr1-1+_n
rename gvkey1 gvkey
keep pdpass gvkey year

然后再keep begyr2-endyr2，drop missing做一次……最后五个merge到一块儿……

如果哪位高手有简洁的方法，5个币感谢优化

*试一试
use dynass.dta,clear
gen x=max(endyr1,endyr2,endyr3,endyr4,endyr5)
gen n=x-begyr1
expand n+1
bys pdpass:gen year=begyr1+_n-1
gen gvkey=gvkey1
forvalue i=2/5{
replace gvkey=gvkey`i' if begyr`i'<=year & year<=endyr`i'+1
}
keep pdpass year gvkey

aolei 发表于 2012-3-7 21:04
*试一试
use dynass.dta,clear
gen x=max(endyr1,endyr2,endyr3,endyr4,endyr5)

非常感谢帮助

我还有两个问题：
1.用max减begyr1忽略了年份不连续的情况，比如endyr1是1998，begyr2是2003，那么99-02的数据就被凭空造出来了……
2.forvalue过程中的if判断，最后的+1是起什么效果呢？比如时间段是91-93，94-96，97-05.那比如97年，它符合 begyr2(1994)<=1997 & 1997<=endyr2+1(1997)，同时也符合begyr3(1997)<=1997 & 1997<=endyr3+1(2006).但我到data editor去看，97确实被分配了gvkey3，结果是对的。这里面是什么原因呢？还望解惑，多谢

sungmoo 发表于 2012-3-7 22:51
use dynass, clear
reshape long pdpco begyr endyr gvkey, i(pdpass source)
keep if gvkey

哇……这个代码果然简洁，我觉得大概不会有更好的了吧

che_nax 发表于 2012-3-8 10:38
非常感谢帮助

我还有两个问题：

1、没注意到年份不连续，最后总的结果应该可以根据缺失情况keep或者drop一下
2、forvalue是个循环语句，每次只对特定的i起作用
3、崇拜sungmoo