r - R 可以可视化 t.test 或其他假设检验结果吗？

Question

我需要在 R 中使用许多假设检验并呈现结果。下面是一个例子:

> library(MASS)
> h=na.omit(survey$Height)
> 
> pop.mean=mean(h)
> h.sample = sample(h,30)
> 
> t.test(h.sample,mu=pop.mean)

    One Sample t-test

data:  h.sample
t = -0.0083069, df = 29, p-value = 0.9934
alternative hypothesis: true mean is not equal to 172.3809
95 percent confidence interval:
 168.8718 175.8615
sample estimates:
mean of x 
 172.3667 

有什么方法可以使 t.test 或其他假设检验结果可视化？

以下是我正在寻找的示例:

Answer 1

你可以做很多事情。这只是我从标准正态分布中抽取随机样本，然后进行 t 检验，绘制观察到的 t 和拒绝均值等于 0 的零假设所需的图。

N=20 #just chosen arbitrarily
samp=rnorm(N)
myTest=t.test(samp)
tcrit=qt(0.025, df=(N-1))

dum=seq(-3.5, 3.5, length=10^4)#For the plot

plot(dum, dt(dum, df=(N-1)), type='l', xlab='t', ylab='f(t)')
abline(v=myTest$statistic, lty=2)
abline(v=tcrit, col='red', lty=2)
abline(v=-tcrit, col='red', lty=2)

当然，每次重新运行此代码时，您观察到的 t 看起来都会有所不同，如果重复运行，这可能是一个很好的说明。