[统计学]简单方便!用统计分析R软件计算置信区间


[统计学]简单方便!用统计分析R软件计算置信区间

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行命令。

统计分析R软件

    「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, ...)