【R每日一贴】R检验样本均值（t检验）？

文/陆勤

问题描述：从总体中获得一个样本，根据这个样本，推知总体的均值是否相当于一个特殊值m？

解决方案：使用R中的t.test函数，对样本进行t检验。即t.test(x, mu=m)，结果输出包括一个p值。按照惯例，若p<0.05，则总体均值不太可能为m；反之，则不能提供均值不等于m的证据。【注意：若是样本量很小，为了保证t检验的结果有意义，潜在总体必须是正态分布。多大样本量为“小”呢？一个经验法则是指样本量<30】

举例说明：从均值为100的正态总体来模拟抽样，然后采用t检验来检测总体均值是否可能为95。

> x <- rnorm(50, mean=100, sd=15)

> t.test(x, mu=95)

      One Sample t-test

data:  x

t = 4.6915, df = 49, p-value = 2.208e-05

alternative hypothesis: true mean is not equal to 95

95 percent confidence interval:

 100.7539 109.3766

sample estimates:

mean of x

 105.0652

【结论】p值很小，所以95不可能是总体均值。

【想一想】

1 t检验的如何定义？有什么用途？

【做一做】

1 查阅t.test函数的R帮助文档和应用实例？

【温馨话语】

“坚忍有恒”四字，最为办事要着，勿畏难中辍，勿滋生弊窦，勿遗人口实。照此进行，何患不达目的？

我是陆勤，微信号：wangluqing360，你们在R路上的朋友，一起欣赏R路上的风景。PPV课R语言群：342699237，欢迎朋友们加入，共建最好的R社区。

赠人玫瑰，手有余香。若是觉得此文有用，欢迎分享给更多的人，让更多的人受用。