R语言单一样本t检验案例实现,练习随机数、条形图制作

今天学习R语言实现单一样本t检验过程。

某个医生检测40名从事铅作业工人的血红蛋白含量,其均数为130.83g/L,标准差为25.74g/L,试分析从事铅作业的工人血红蛋白含量是否不同于正常成年人平均值140g/L?

没有40个样本的原始监测数据怎么办?用R生成随机数。

set.seed(1234)
x <- rnorm(n=40,mean=130.83,sd=25.74)
hist(x,col="light blue")

转换成密度直方图,添加正态曲线,

hist(x,prob=T,col="light blue")
lines(density(x), col="red", lwd=3)

正态曲线形状呈现倒扣钟形,粗略来看样本来自正态分布总体(这是废话,但由于t检验数据需满足正态性条件,所以还得唠叨一句)。

单一的一组样本,已知总体均值,所以采用单一样本t检验即可考察40例工人和正常人血红蛋白含量是否有差异。

t.test(x,mu=140)

来看统计检验的结果:

t=-5.3357,p=4.321e-06,以a=0.05水准,p<0.05,有统计意义,结合专业认为从事铅作业的工人平均血红蛋白含量低于正常人。(130.83<140)

条形图可视化效果:

(0)

相关推荐