songzien 发表于 2011-6-3 17:11
两个总体均值之差的估计和置信区间用stata怎么算啊
*估计两个同方差的正态分布变量x1与x2的期望差(x1与x2的样本量可以不同):
forv i=1/2{
su x`i'
sca n`i'=r(N)
sca m`i'=r(mean)
sca s`i'=r(Var)
}
sca s= invttail(n1+n2-2,.025)*sqrt(((n1-1)*s1+(n2-1)*s2)*(1/n1+1/n2)/(n1+n2-2))
di "点估计:" m1-m2
di "95%置信区间:(" m1-m2-s " , " m1-m2+s ")"