急问：遇到一个关于计算标准差的难题

g s=sum(x)
g s2=sum(x^2)
replace s=sqrt(((s2-s2[_n-125])-(s-s[_n-125])^2/125)/124) in 126/l
replace s=sqrt((s2-s^2/125)/124) in 125
replace s=. in f/124
drop s2

或

set ob 251
g t=_n
tsset t
tssmooth ma s=x,w(125)
tssmooth ma s2=x^2,w(125)
g sd=sqrt((s2[_n+125]-s[_n+125]^2)*125/124)
drop s s2 t
drop in 251

sungmoo 发表于 2011-9-1 23:54
g s=sum(x)
g s2=sum(x^2)
replace s=sqrt(((s2-s2[_n-125])-(s-s[_n-125])^2/125)/124) in 126/l

呵呵~~高手！！

sungmoo 发表于 2011-9-1 23:54
g s=sum(x)
g s2=sum(x^2)
replace s=sqrt(((s2-s2[_n-125])-(s-s[_n-125])^2/125)/124) in 126/l

3Q！！！！

liuxb 发表于 2011-8-31 15:25
*1. if group=1…k and data is long form:

capture drop sd

求助高手！在stata中，如何生成一个变量如（主营业务收入）三年的标准差？比如下面例子：
代码       年份 收入
000002  1999     10
000002  2000     5
000002  2001     12
000002  2002     7
000002  2003     6
000002  2003     8
000005  1999     4
000005  2000    5
000005  2001    8
000005  2002  10
000005  2003 12
我想生成一列变量等于收入连续三年的标准差，该如何操作呢？十分感谢啊！

sungmoo 发表于 2011-9-1 23:54
g s=sum(x)
g s2=sum(x^2)
replace s=sqrt(((s2-s2[_n-125])-(s-s[_n-125])^2/125)/124) in 126/l

如果有多组，每组都是250个观测值，仍然从第1个观测值开始，计算第1个至第125个观测值的样本标准差，然后再计算第2个到第126个观测值的样本标准差，以此类推，逐次后退一个观测值计算。每一组都这样，如何实现？ 第一个程序不能用了吧 因为in不能和by一起用的，是吗？

*另一种思路：
https://bbs.pinggu.org/thread-1210664-1-1.html