stata循环命令问题

现在有变量X_1,x_2,x_3,x_4,x_5,x_6，要生成5个新变量y_1=x_1*x_2*x_3*x_4*x_5*x_6;y_2=x_2*x_3*x_4*x_5*x_6;…；y5=x_5*x_6,如何用一个循环命令解决上述问题?感谢大神点！！！

Nick Cox 建议

1. * Example generated by -dataex-. To install: ssc install dataex
   
2. clear
   
3. input float(x_1 x_2 x_3 x_4 x_5 x_6)
   
4. 3.6439085  .8855875  2.178382    .13567242  9.159781  .5004264
   
5. 4.6231413 .02498507  3.880392     .7207576  8.755236  .5516323
   
6.   5.80438 .12259632   5.36664     .0965251  4.605244  .1663715
   
7. 7.822571  .5383945  .8250254 .00018472524   5.31311  .4234573
   
8. .6603088 .29715487 1.0873007    .35380825  3.255852 .50436085
   
9. 2.589655 .05635693  7.182124    .27086043  5.187715  .6142209
   
10. 5.666703  .7410993 2.4030285     .5832268  4.603475  .4479784
    
11. 1.876804  .6852789  7.749197     .8524442  7.455486 .14944874
    
12. 4.2674637  .9836587  9.832918     .8806648 1.1855963  .8706222
    
13. .8688191  .6941437  .1651141     .0421818  5.836682 .27148572
    
14. end
    
15. 
    
16. * code
    
17. forval j = 1/5 {
    
18.   gen y_`j' = x_6
    
19.   quietly forval k = 5(-1)`j' {
    
20.     replace y_`j' = y_`j' * x_`k'
    
21.   }
    
22. }

复制代码

非常感谢，[em20][em20][em20][em20][em20]

黃河泉 发表于 2017-3-6 17:26
Nick Cox 建议

黄老师您好！
请帮忙解决循环只执行一次，不能向下执行的问题。谢谢您！到第二个g,drop后的 (g==`t') 就不能变化执行了。试了，不是drop不满足条件中断。
请问这个命令因为什么问题，g=1执行之后就不能执行了呢？g=2更改为符合条件的也不能执行？谢谢

解决的问题：同一公司年度有多个s0701a==1或s0701a==2的seperation值，保留唯一seperation值。
第一步，s0701a==1，seperation不缺失，删除s0701a=2的所有行观测值；然后取s0701a==1，seperation的均值，保留一个seperation值。
第二步，s0701a==1，seperation缺失，删除s0701a==1的所有行观测值；然后取s0701a==2，seperation的均值，保留一个seperation值。

数据如下：

1. * Example generated by -dataex-. To install: ssc install dataex
   
2. clear
   
3. input str10 codeyear byte s0701a double seperation byte r_tag double(g seperation0 seperation1)
   
4. "0000012003" 1       0 1 1       . .
   
5. "0000012003" 2       0 1 1       0 .
   
6. "0000012004" 1       . 4 2       . .
   
7. "0000012004" 2       . 4 2       . .
   
8. "0000012004" 2       . 4 2       . .
   
9. "0000012004" 2       . 4 2       . .
   
10. "0000012004" 2       . 4 2       . .
    
11. "0000012016" 1       . 1 3       . .
    
12. "0000012016" 2  4.0318 1 3  4.0318 .
    
13. "0000012017" 2  4.0696 1 4  4.0696 .
    
14. "0000042020" 1  6.8294 2 5       . .
    
15. "0000042020" 2  3.9358 2 5  3.4147 .
    
16. "0000042020" 2  2.8936 2 5  3.4147 .
    
17. "6030002013" 1       0 0 6       . .
    
18. "6038832015" 1       0 3 7       . .
    
19. "6038832015" 1   .2434 3 7       . .
    
20. "6038832015" 1       0 3 7       . .
    
21. "6038832015" 2 10.6729 3 7 10.6729 .
    
22. end

复制代码

不完整的循环命令，不能执行一次后不能执行：

egen g=group(codeyear )          // g分组
sum g
bys codeyear:  egen  seperation0=mean(seperation)  if  s0701a==2
gen seperation1=.
  cap label variable year "年度"
  cap label variable codeyear "代码年度，唯一标识"
  cap label variable seperation "两权分离率(%)_没有单位"
  cap label variable  s0701a     "1=年报中公布，2=根据股权控制链计算所得"
  cap label variable  r_tag         " codeyear 重复次数"
    dataex
forvalues t=1/7  {
set trace on
   if  (g==`t')  &  (s0701a==1)  & (seperation!=. ){
    drop  if  (g==`t') & (s0701a==2)
     }

你看看这是不是你要的结果

1. * Example generated by -dataex-. For more info, type help dataex
   
2. 
   
3. collapse seperation,by(codeyear s0701a)
   
4. bys codeyear:gen wanted=seperation[1] if !mi(seperation[1])
   
5. bys codeyear:replace wanted=seperation[2] if mi(seperation[1])
   
6. collapse wanted,by(codeyear)
   
7. 
   
8. list
   
9. 
   
10.     +------------------------+
    
11.      |   codeyear      wanted |
    
12.      |------------------------|
    
13.   1. | 0000012003           0 |
    
14.   2. | 0000012004           . |
    
15.   3. | 0000012016      4.0318 |
    
16.   4. | 0000012017      4.0696 |
    
17.   5. | 0000042020      6.8294 |
    
18.      |------------------------|
    
19.   6. | 6030002013           0 |
    
20.   7. | 6038832015   .08113333 |
    
21.      +------------------------+

复制代码

zdlspace 发表于 2021-8-10 03:17
你看看这是不是你要的结果

谢谢黄老师！帮助很大！
不过，还有两个问题请教：（1）我的循环为什么只能执行一次，之后，如上图，组变量g不根据循环变化？（2）collapse后，seperation只剩下一个均值，如何保留原来的变量和其他许多变量？

wanliliu 发表于 2021-8-12 19:36
谢谢黄老师！帮助很大！
不过，还有两个问题请教：（1）我的循环为什么只能执行一次，之后，如上图，组变 ...

我不是给你写了代码了吗？不用请教黄老师了，黄老师已经退出论坛了

wanliliu 发表于 2021-8-12 19:36
谢谢黄老师！帮助很大！
不过，还有两个问题请教：（1）我的循环为什么只能执行一次，之后，如上图，组变 ...

你写的循环有问题，你只告诉forval如果满足 (s0701a==1)  & (seperation!=. )时怎么运算，并没有告诉它，如果不满足怎么算。关于collapse如果你要保持原来的所有数据，可以使用gcollapse。代入如下：

1. clear
   
2. input str10 codeyear byte s0701a double seperation byte r_tag double(g seperation0 seperation1)
   
3. "0000012003" 1       0 1 1       . .
   
4. "0000012003" 2       0 1 1       0 .
   
5. "0000012004" 1       . 4 2       . .
   
6. "0000012004" 2       . 4 2       . .
   
7. "0000012004" 2       . 4 2       . .
   
8. "0000012004" 2       . 4 2       . .
   
9. "0000012004" 2       . 4 2       . .
   
10. "0000012016" 1       . 1 3       . .
    
11. "0000012016" 2  4.0318 1 3  4.0318 .
    
12. "0000012017" 2  4.0696 1 4  4.0696 .
    
13. "0000042020" 1  6.8294 2 5       . .
    
14. "0000042020" 2  3.9358 2 5  3.4147 .
    
15. "0000042020" 2  2.8936 2 5  3.4147 .
    
16. "6030002013" 1       0 0 6       . .
    
17. "6038832015" 1       0 3 7       . .
    
18. "6038832015" 1   .2434 3 7       . .
    
19. "6038832015" 1       0 3 7       . .
    
20. "6038832015" 2 10.6729 3 7 10.6729 .
    
21. end
    
22. 
    
23. 
    
24. gcollapse seperation,by(codeyear s0701a) merge replace
    
25. drop if  s0701a==1 & mi(seperation)
    
26. gen d1=seperation if s0701a==1
    
27. egen tag=min(d1),by(codeyear)
    
28. drop if mi(d1) &!mi(tag)
    
29. 
    
30. list ,noobs
    
31. 
    
32. +-------------------------------------------------------------------+
    
33.   |   codeyear   s0701a   seperat~n   r_tag   g   sepera~0   sepera~1 |
    
34.   |-------------------------------------------------------------------|
    
35.   | 0000012003        1           0       1   1          .          . |
    
36.   | 0000012004        2           .       4   2          .          . |
    
37.   | 0000012004        2           .       4   2          .          . |
    
38.   | 0000012004        2           .       4   2          .          . |
    
39.   | 0000012004        2           .       4   2          .          . |
    
40.   |-------------------------------------------------------------------|
    
41.   | 0000012016        2      4.0318       1   3     4.0318          . |
    
42.   | 0000012017        2      4.0696       1   4     4.0696          . |
    
43.   | 0000042020        1      6.8294       2   5          .          . |
    
44.   | 6030002013        1           0       0   6          .          . |
    
45.   | 6038832015        1   .08113333       3   7          .          . |
    
46.   |-------------------------------------------------------------------|
    
47.   | 6038832015        1   .08113333       3   7          .          . |
    
48.   | 6038832015        1   .08113333       3   7          .          . |
    
49.   +-------------------------------------------------------------------+

复制代码

zdlspace 发表于 2021-8-13 14:58
你写的循环有问题，你只告诉forval如果满足 (s0701a==1)  & (seperation!=. )时怎么运算，并没有告诉它， ...

非常感谢您！大牛！
如果不满足条件，什么都不用做，命令是不是只写else,后面什么都不用写？

zdlspace 发表于 2021-8-13 14:58
你写的循环有问题，你只告诉forval如果满足 (s0701a==1)  & (seperation!=. )时怎么运算，并没有告诉它， ...

谢谢您！ gcollapse seperation,by(codeyear s0701a) merge replace
该命令运行后只有一行值了，后续没法做