【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社区。
赠人玫瑰,手有余香。若是觉得此文有用,欢迎分享给更多的人,让更多的人受用。
请关注“恒诺新知”微信公众号,感谢“R语言“,”数据那些事儿“,”老俊俊的生信笔记“,”冷🈚️思“,“珞珈R”,“生信星球”的支持!