求助~stata多变量回归时循环语句怎么写？谢谢大家！

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

求助~stata多变量回归时循环语句怎么写？
因变量的指标有Y1，Y2……Y6等6个，
自变量的指标有X1，X2……X4等4个;
自变量指标的平方项为：X12，X22....X42;
调节变量有两个M和N，M有3个指标M1，M2，M3;N有3个指标N1，N2，N3;
算出每个自变量指标及平方项与每个调节变量的交互项：X1_M1; X12_M1……;X1_N1;X12_NI......

要求回归每个因变量指标与自变量、自变量平方、两个调节变量交互项、二次交互项的结果。

感觉运算量有点大，有没有大神会用stata循环语句救救我啊？谢谢谢谢！谢大家！

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏1 回帖

关键词：Stata 循环语句 tata 多变量 调节变量 因变量 自变量

有么有人知道啊。。。

fuzixi1125 发表于 2016-8-26 17:39
有么有人知道啊。。。

试试

1. clear
   
2. 
   
3. set obs 1000
   
4. 
   
5. gen y1 = rnormal(5,2.25)
   
6. 
   
7. gen x1 = runiform()
   
8. gen x2 = rbeta(1,2)
   
9. gen x3 = rbinomial(100,0.3)
   
10. gen x4 = runiform()
    
11. 
    
12. gen x12 = x1^2
    
13. gen x22 = x2^2
    
14. gen x32 = x3^2
    
15. gen x42 = x4^2
    
16. 
    
17. gen m1 = rnormal(0,1)
    
18. gen m2 = rnormal(1.5,2)
    
19. gen m3 = rnormal(-0.5,3)
    
20. 
    
21. gen n1 = rnormal(1,3)
    
22. gen n2 = rnormal(1,2.2)
    
23. gen n3 = rnormal(-1,1.2)
    
24. 
    
25. local x "x1 x2 x3 x4"
    
26. local x2 "x12 x22 x32 x42"
    
27. 
    
28. foreach v of varlist `x' `x2' {
    
29.   gen `v'_m1 = `v'*m1
    
30.   gen `v'_m2 = `v'*m2
    
31.   gen `v'_m3 = `v'*m3
    
32.   gen `v'_n1 = `v'*n1
    
33.   gen `v'_n2 = `v'*n2
    
34.   gen `v'_n3 = `v'*n3
    
35. }  
    
36.   
    
37. reg y1 `x' `x2' x1_m1-x42_n3

复制代码

黃河泉 发表于 2016-8-27 08:40
试试

哈哈，谢谢呢。。。

黃河泉 发表于 2016-8-27 08:40
试试

能否再请教个问题，已经用xi i.industry 生成了一个行业虚拟变量，想用xi i.city 再生成一个区域虚拟变量，但是一旦运行xi i.city,原来xi i.industry 生成的的那些虚拟变量没有了，请教是怎么回事？应该怎么修改？谢谢~

fuzixi1125 发表于 2016-8-27 11:04
能否再请教个问题，已经用xi i.industry 生成了一个行业虚拟变量，想用xi i.city 再生成一个区域虚拟变量 ...

会造成这种结果应该就是"线性重合"的问题，当然就是 city 与 industry 之间有一些特殊关系！我看不到资料，所以只能"猜测"到这里！