曾经一个英国专家留给我们的命令如下,是用stata7.0来做的标准误,你可以看看把某些变量名改一下运行运行试试。
cap program drop jackn
program def jackn
version 7
tempname simulation
tempfile jknresults
postfile `simulation' valuex psx x using `jknresults', replace
local i=1
while `i'<=15 {
use temp, clear
qui summ mathscorehd [aw=studentsweight]
local x=r(mean)
qui drop if newpsu==`i'
qui gen weight1=(15/14)*studentsweight
qui summ mathscorehd [aw=weight1]
local psx=r(mean)
local valuex=(15/14)*((`psx'-`x')^2)
post `simulation' (`valuex') (`psx') (`x')
local i=`i'+1
}
postclose `simulation'
drop _all
use `jknresults'
egen stderrx=sum(valuex)
qui replace stderrx=sqrt(stderrx)
summ stderrx psx x
end