受教了，感谢。
我之前的写法是这样的，我感觉至少在代码长度上要简洁一点，不过也有点缺陷，放在这交流一下。
gen P_skew=.
gen P_kurt=.
tostring Liscd, generate(Liscd1)
encode Liscd1, gen(Liscd2)
forvalues i=1/3156 {
sktest r if Liscd2==`i'
replace P_skew=r(P_skew) if Liscd2==`i'
replace P_kurt=r(P_kurt) if Liscd2==`i'
}
代码里那个3156是在得到Liscd2以后自己查看到的Liscd2的最大值，表示有3156只证券。这个也是问题所在，这一步必须要手工进行。现在看了您的group(Liscd)和forvalues这两行以后相信会有改进空间，多谢。
--------------------------------------
Number of by-groups = 3,156
by-groups with errors = 3,156
by-groups with no data = 0
Observations processed = 76,606
Observations saved = 0
--------------------------------------
这是为什么呢
扫码加好友，拉您进群
全部版块
我的主页
收藏
