如何解决如何将p值绘制到标准法线曲线上?
我试图从视觉上说明一个基本的两尾假设检验。我有学生使用ttest来找到p值。我希望能够绘制标准正态分布并为大于t统计量的区域着色。我使用以下方法生成标准法线:
#delimit ;
graph twoway (function y=normalden(x,1),range(-4 4) lw(medthick)),title("Standard Normal-Distribution (Null)")
xtitle("Normal",size(medlarge)) ytitle("")
xlabel(-4(1)4)
xscale(lw(medthick)) yscale(lw(medthick))
legend(off)
graphregion(fcolor(white));
#delimit cr
产生
但是,我想弄清楚如何在ttest之后填写p值的区域。理想情况下,我希望此示例通常可应用于任何两面ttest。我正在尝试获取如下所示的内容:
我是通过在Paint的基本区域中填充来实现的,但是我正在尝试获得更加自动化和精确的内容。
解决方法
感谢Nick Cox指出这个确切的示例在手册条目twoway function
中。
以下内容提供了我所需的信息:
ttest highestPrev,by(match) // can be switched out for any ttest of interest
local positive_tstat = abs(r(t))
local negative_tstat = -1*abs(r(t))
#delimit ;
twoway
function y=normalden(x),range(-4 `negative_tstat') color(gs12) recast(area)
|| function y=normalden(x),range(`positive_tstat' 4) color(gs12) recast(area)
|| function y=normalden(x),range(-4 4) lstyle(foreground)
||,plotregion(style(none))
ysca(off) xsca(noline)
legend(off)
xlabel(-4 "-4 sd" -3 "-3 sd" -2 "-2 sd" -1 "-1 sd" 0 "0"
1 "1 sd" 2 "2 sd" 3 "3 sd" 4 "4 sd",grid gmin gmax)
xtitle(" " "Standard Normal Distribution (Null)")
;
#delimit cr