如何将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

产生

this link

但是，我想弄清楚如何在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