[统计学]简单方便!用统计分析R软件计算置信区间
在某个螺丝工厂工作的小王,接到了领到指示“A型号螺丝的长度平均值给我调查一下,今天下班之前给我!”
A型号螺丝有10万根,而且都在外地的仓库保管,小王手头只有10根螺丝。
小王内心暗骂上司又给自己找难题,但是又不能不干,总之先把手头10根螺丝的长度进行了测量,并计算了平均值。
但是,把这10根螺丝的平均当成10万根螺丝的平均值报告给上司的话,确实太过分了。于是小王找李师傅商量。李师傅告诉他,只要求总体平均的置信区间,并汇报给上司就可以了。
你能帮一下小王吗?(置信水平95%)
螺丝长度及平均值
No. | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 平均 |
---|---|---|---|---|---|---|---|---|---|---|---|
长度(mm) | 31 | 29 | 28 | 35 | 30 | 30 | 32 | 27 | 33 | 29 | 30.4 |
关于置信区间,置信水平请参考→[统计学]推论统计中的区间估计是什么?如何理解置信水平和置信区间?
关于置信区间的计算方法请参考→[统计学]一文看懂如何计算置信区间
参考上面2篇文章,动笔计算也可以帮助小王解决困难。但是这次我们用统计分析R软件来帮助小王解决问题。在你完全掌握背后原理之后,用适当的工具解决实际问题才能事半功倍。
使用统计分析R软件的话,仅仅只需要2行命令,就可以得到置信区间。
统计分析R软件
R是用于统计分析、绘图的语言和操作环境。R是属于GNU系统的一个自由、免费、源代码开放的软件,它是一个用于统计计算和统计制图的优秀工具。
# <-类似于别的语言中的=,赋值符号
screw <- c(31, 29, 28, 35, 30, 30, 32, 27, 33, 29)
# t分布,置信水平95%的置信区间
t.test(screw, conf.level=0.95)
只要在下列命令行中输入上述2行命令。
「95 percent confidence interval」下方的28.6739,32.1261即为置信水平95%的置信区间。
本来需要笔算半天又要查表,只需要2行代码就解决了!有兴趣的话,也可以再详细研究一下t.test函数。
# Description
# Performs one and two sample t-tests on vectors of data.
# Usage
t.test(x, ...) ## Default S3 method: t.test(x, y = NULL, alternative = c("two.sided", "less", "greater"), mu = 0, paired = FALSE, var.equal = FALSE, conf.level = 0.95, ...) ## S3 method for class 'formula' t.test(formula, data, subset, na.action, ...)