求助：面板数据同年的各截面相除穷举法。

一、数据说明
数据参考如下的官网范例数据【grunfeld.dta】：年份1935-1954年，10个个体，共200个观测值。

1. webuse  grunfeld ,  clear

复制代码

截面ID为company，年份为year，另外就是三个变量：invest、mvalue、kstock。

company	year	invest	mvalue	kstock
1	1935	I1,1935	M1,1935	K1,1935
…	…	…	…	…
1	t	I1,t	M1,t	K1,t
…	…	…	…	…
1	1954	I1,1954	M1,1954	K1,1954
……	……	……	……	……
i	1935	Ii,1935	Mi,1935	Ki,1935
…	…	…	…	…
i	t	Ii,t	Mi,t	Ki,t
…	…	…	…	…
i	1954	Ii,1954	Mi,1954	Ki,1954
……	……	……	……	……
j	1935	Ij,1935	Mj,1935	Kj,1935
…	…	…	…	…
j	t	Ij,t	Mj,t	Kj,t
…	…	…	…	…
j	1954	Ij,1954	Mj,1954	Kj,1954
……	……	……	……	……

二、欲实现的效果
最后想要得到一个【新的数据集】，每年每截面的变量值均为相除得到。也就是说：
如果按照有原本20年数据，10个个体，共计200个观测值；
最后得到的将是10*9=90个新个体/截面（newid字符型），同样20年，共计1800个新观测值。

company	newid （字符型）	year	invest	mvalue	kstock
1	1by2	1935	I1,1935/I2,1935	…	K1,1935/K2,1935
…	…	…	…	…	…
1	1by2	t	I1,t/I2,t	…	K1,t/K2,t
…	…	…	…	…	…
1	1by2	1954	I1,1954/I2,1935	…	K1,1954/K2,1954
……	……	……	……	……	……
…	…	…	…	…	…
1	1byi	t	I1,t/Ii,t	…	M1,t/M2,t
…	…	…	…	…	…
……	……	……	……	……	……
i	ibyj	1935	Ii,1935/Ij,1935	…	Ki,1935/Kj,1935
…	…	…	…	…	…
i	ibyj	t	Ii,t/Ij,t	`	Ki,t/Kj,t
…	…	…	…	…	…
i	ibyj	1954	Ii,1954/Ij,1954	…	Ki,1954/Kj,1954
……	…	……	……	……	……

Question：请教这个新数据集怎么得到呢？

不好意思，错了一个地方：
set obs 1800
gen invest2=.
gen mvalue2=.
gen kstock2=.
gen newid=""
local k=1
forvalue i=1/200{
forvalues j=1/200{
if `i'!=`j'&year[`i']==year[`j']{
replace invest2[`k']=invest[`i']/invest[`j']
replace mvalue2[`k']=mvalue[`i']/mvalue[`j']
replace kstock2[`k']=kstock[`i']/kstock[`j']
replace newid="`i'by`j'"
dis `k++'
}
}
}

你好，我做出一个程序如下：
set obs 1800
gen invest2=.
gen mvalue2=.
gen kstock2=.
gen newid=""
local k=1
forvalue i=1/200{
forvalues j=1/200{
if `i'!=`j'&year[`i']==year[`j']{
replace invest2[`k']=invest[`i']/invest[`j']
replace mvalue2[`k']=invest[`i']/invest[`j']
replace kstock2[`k']=invest[`i']/invest[`j']
replace newid="`i'by`j'"
dis `k++'
}
}
}
不知道可不可以运行，如果不能运行，希望你把数据发给我，我来调试。我的邮箱是812466144@qq.com

1. webuse grunfeld, clear
   
2. 
   
3. set obs 1800
   
4. 
   
5. gen invest2=.
   
6. gen mvalue2=.
   
7. gen kstock2=.
   
8. gen newid=""
   
9. gen cd1=.
   
10. gen cd2=.
    
11. gen company1=.
    
12. gen company2=.
    
13. gen year1=.
    
14. 
    
15. local k=1
    
16. local g=1
    
17. local f=1
    
18. 
    
19. forvalue i=1/200 {
    
20.        forvalues j=1/200 {
    
21.               if `i'!=`j'  &  year[`i']==year[`j'] {
    
22.                      replace invest2=invest[`i']/invest[`j'] in `k'
    
23.                      replace mvalue2=mvalue[`i']/mvalue[`j'] in `k'
    
24.                      replace kstock2=kstock[`i']/kstock[`j'] in `k'
    
25.                      replace cd1=`i' in `k'
    
26.                      replace cd2=`j' in `k'
    
27.                      replace company1=company[`i'] in `k'
    
28.                      replace company2=company[`j'] in `k'
    
29.                      replace year1=year[`i'] in `k'
    
30.                      dis `k++'
    
31.               }
    
32.        }
    
33. }
    
34. 
    
35. tostring company1 company2,replace
    
36. replace newid=company1+"by"+company2
    
37. drop company-time
    
38. gen company=real(company1)
    
39. rename year1 year
    
40. order company year
    
41. drop cd1-company2

复制代码

以上是zhonghuijie提供的代码，经验证准确无误！