急问：遇到一个关于计算标准差的难题 - Stata专版 - 经管之家

› 论坛 › 计量经济学与统计论坛五区 › 计量经济学与统计软件 › Stata专版

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

9222

25

收藏 2011-08-31

变量有250个观测值，想从第1个观测值开始，计算第1个至第125个观测值的样本标准差，然后再计算第2个到第126个观测值的样本标准差，以此类推，逐次后退一个观测值计算。想了一大半天，实现不了。老板急着要结果。。
非常非常感谢！！！！！在线等。。。
PS:昨天一着急发到SAS版了，两位老大帮忙，发了代码，可惜用不上，angway，在这里，再一次，非常感谢！！！

二维码

扫码加我拉你入群

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

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

全部回复

2011-8-31 09:27:32

matlab啊，
假设data存的是250个数据。
for i= 1:125
compute = data[i:i+125];
res(i) = std(compute);
end
最后res中的数据就是了。

二维码

扫码加我拉你入群

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

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

2011-8-31 09:29:53

std是matlab有的函数

二维码

扫码加我拉你入群

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

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

2011-8-31 09:46:44

老板就是要我用stata算啊，因为前期都是用的stata，悲催。。。

二维码

扫码加我拉你入群

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

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

2011-8-31 09:48:05

bithuye 发表于 2011-8-31 09:29
std是matlab有的函数

谢谢！！！！因为被要求用stata算，所以。。。用stata的循环能够实现吗？

二维码

扫码加我拉你入群

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

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

2011-8-31 11:57:26

试试下面的程序：
capture drop sd
gen sd=.
local nloop=250-125+1
forvalue i=1(1)`nloop' {
qui sum x if _n>=`i' & _n<=125+(`i'-1)
qui replace sd=r(sd) if _n==`i'
}
list sd if sd!=.

二维码

扫码加我拉你入群

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

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

点击查看更多内容…

2011-8-31 14:00:01

liuxb 发表于 2011-8-31 11:57
试试下面的程序：
capture drop sd
gen sd=.

解决了谢谢大牛

！！！！
ps:想提高自己的stata编程能力，可否推荐一些教材或者资料？？谢谢！！

二维码

扫码加我拉你入群

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

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

2011-8-31 14:05:04

liuxb 发表于 2011-8-31 11:57
试试下面的程序：
capture drop sd
gen sd=.

还想再请教一下，如果有多组，每组都是250个观测值，仍然从第1个观测值开始，计算第1个至第125个观测值的样本标准差，然后再计算第2个到第126个观测值的样本标准差，以此类推，逐次后退一个观测值计算。每一组都这样，如何实现？

二维码

扫码加我拉你入群

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

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

2011-8-31 15:25:30

*1. if group=1…k and data is long form:

capture drop sd
qui gen sd=.
local nloop=250-125+1

qui sum group
local ngroup=r(max)

forvalue k=1(1)`ngroup'{
local istart=(`k'-1)*250+1
local iend=(`k'-1)*250+`nloop'

forvalue i=`istart'(1)`iend' {
      qui sum x if _n>=`i' & _n<=125+(`i'-1)
      qui replace sd=r(sd) if _n==`i'
}
}

bysort group:list sd if sd!=.

*2. if group=1…k and data is wide form:
* For example: k=3

local nloop=250-125+1

foreach var of varlist x1-x3 {
capture drop sd`var'
qui gen sd`var'=.

forvalue i=1(1)`nloop' {
      qui sum `var' if _n>=`i' & _n<=125+(`i'-1)
      qui replace sd`var'=r(sd) if _n==`i'
}
}

list sd* if _n<=126

二维码

扫码加我拉你入群

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

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

2011-8-31 15:47:22

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

capture drop sd

出手好快啊！！先谢谢！！！马上试试

二维码

扫码加我拉你入群

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

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

2011-8-31 16:01:00

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

capture drop sd

capture drop sd
qui gen sd=.
local nloop=250-125+1

qui sum sid
local nsid=r(max)

forvalue k=1(1)`nsid'{
local istart=(`k'-1)*250+1
local iend=(`k'-1)*250+`nloop'

forvalue i=`istart'(1)`iend' {
      qui sum rpercent if _n>=`i' & _n<=125+(`i'-1)
      qui replace sd=r(sd) if _n==`i'
}
}

把变量换成我要处理的，报告语法无效（invalid syntax）

二维码

扫码加我拉你入群

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

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

2011-8-31 16:02:42

xixi221zmw 发表于 2011-8-31 15:47
出手好快啊！！先谢谢！！！马上试试

还是不行哦。。。。。。。。。

二维码

扫码加我拉你入群

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

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

2011-8-31 16:30:54

程序运行没有错误。
把你的数据传上来看看！

二维码

扫码加我拉你入群

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

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

2011-8-31 16:48:14

是不是因为有的数据不是250？

附件列表

大小:3.87 MB

二维码

扫码加我拉你入群

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

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

2011-8-31 17:15:56

1. invalid syntax 是因为你的数据sid不是数值型变量；
2. 我编的程序是根据每组有250观察，需要循环126次，如果一组观察个数小于126，那么在保存的结果中会重迭，造成结果错误。

二维码

扫码加我拉你入群

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

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

2011-8-31 17:17:37

liuxb 发表于 2011-8-31 17:15
1. invalid syntax 是因为你的数据sid不是数值型变量；
2. 我编的程序是每组必须有250观察，否则结果不对。 ...

哦~~那我把sid转成数值型

二维码

扫码加我拉你入群

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

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

2011-8-31 17:25:41

liuxb 发表于 2011-8-31 17:15
1. invalid syntax 是因为你的数据sid不是数值型变量；
2. 我编的程序是每组必须有250观察，否则结果不对。 ...

谢谢！！下午从你的两个程序了学到很多。THANKS!

二维码

扫码加我拉你入群

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

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

2011-9-1 09:55:37

liuxb 发表于 2011-8-31 17:15
1. invalid syntax 是因为你的数据sid不是数值型变量；
2. 我编的程序是根据每组有250观察，需要循环126次 ...

我把程序改成下面的：
capture drop sd

.
. qui gen sd=.

.
. local nloop=242-121+1

.
.
.
. qui sum sid

.
. local nsid=r(max)

.
.
.
. forvalue k=1(1)`nsid'{
  2.
.    local istart=(`k'-1)*242+1
  3.
.    local iend=(`k'-1)*242+`nloop'
  4.
.
.
.    forvalue i=`istart'(1)`iend' {
  5.
.       qui sum rpercent if _n>=`i' & _n<=121+(`i'-1)
  6.
.       qui replace sd=r(sd) if _n==`i'
  7.
.    }
  8.
. }

每组都是242个观测值，sid已经转成了int型，在11.2SE上跑，没显示错误，但是一直运行不出结果。。。。已经快半小时了。  麻烦你再帮我看看吧，必须是250个观测值？不胜感激！！！

二维码

扫码加我拉你入群

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

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

2011-9-1 09:57:26

liuxb 发表于 2011-8-31 17:15
1. invalid syntax 是因为你的数据sid不是数值型变量；
2. 我编的程序是根据每组有250观察，需要循环126次 ...

结果出来了正确。谢谢大牛

二维码

扫码加我拉你入群

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

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

2011-9-1 11:46:09

liuxb 发表于 2011-8-31 17:15
1. invalid syntax 是因为你的数据sid不是数值型变量；
2. 我编的程序是根据每组有250观察，需要循环126次 ...

再请教一个问题，我用函数对r计算了sd后，用replace值付给另一个变量V，再用函数对p计算sd，为何这时V的值全部变成缺省值了呢？

二维码

扫码加我拉你入群

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

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

2011-9-1 23:54:53

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

二维码

扫码加我拉你入群

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

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

2011-9-3 21:22:43

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

呵呵~~高手！！

二维码

扫码加我拉你入群

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

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

2011-9-3 21:23:04

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！！！！

二维码

扫码加我拉你入群

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

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

2011-9-5 22:30:37

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
我想生成一列变量等于收入连续三年的标准差，该如何操作呢？十分感谢啊！

二维码

扫码加我拉你入群

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

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

2011-9-6 09:43:25

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一起用的，是吗？

二维码

扫码加我拉你入群

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

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

2011-10-19 06:52:25

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

二维码

扫码加我拉你入群

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

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

相关推荐

栏目导航

热门文章

推荐文章

扫码加好友，拉您进群

各岗位、行业、专业交流群