周广肃《stata统计分析与应用》P68,关于z分布应用于假设检验,H0:u=32.5,alpha=.05,其stata code如下:
. quietly sum
. scalar z=(r(mean)-32.5)/ (1.1/sqrt(6))
. scalar p=(1-normal(abs(z))/2
. scalar crit=invnormal(1-.05/2)
. scalar list z p crit
我觉得这里可能有错误。
由于这是two tail test, p值要乘2,然后将p与alpha比较,即正确的code应该为:scalar p=2*(1-normal(abs(z))。
或者,不乘2,然后将p与alpha/2进行比较,即code为:scalar p=1-normal(abs(z))。
作者将p值的计算code设置为 scalar p=(1-normal(abs(z))/2,就比较奇怪了。
不知,我的理解是否正确,尚望大家指正。
谢谢。


雷达卡


学习啦
京公网安备 11010802022788号







