Minitab的基本统计命令与假设测验-SunfineData.PDF

R 的基本统计命令与假设测验 本文介绍了R 的基本统计命令，包括变数的统计描述、单个平均数的统计推断、两 个平均数的假设测验等，并给出了具体实例的R 程序和计算结果的分析。 一、变量的统计描述 （一）命令 R 提供了用来计算变量的各种统计量的命令summary ，其用法和有关的说明如下： 用法：summary(x) 功能：输出向量x 的平均值MEAN ，中位数MEDIAN ，最大值MAX ，最小值MIN ， 四分之三值Q3 ，四分之一值Q1 。 说明：将变量x 的数据由小到大排列，Q3 为3(N+1)/4 位置的值，Q1 为(N+1)/4 位置 的值。 另外，用length 获得样本个数，用stderr 获得样本平均数的标准误。 （二）实例 〖例1〗测定10 个“金皇后”玉米果穗的重量(克) ，得结果为150，170，160，162， 170，168，171，162，151，160。试计算它的平均数、标准差、变异系数等统计量。 1 R程序 # example 4-1 c1 - c(150,170,160,162,170,168,171,162,151,160) summary(c1) # sample size n - length(c1); n # variance s2 - var(c1); s2 # standard deviation s - sqrt(var(c1)); s # standard error of the sample mean se - stderr(c1); se # estimate the coefficient of variation cv - sqrt(var(c1)) /mean(c1) * 100; k1 # geometric mean c2 - log(c1) .1. gmean - exp(mean(c2)) gmean #harmonic mean c3 - 1/c1 hmean - 1/mean(c3) hmean 2 R计算结果 Min. 1st Qu. Median Mean 3rd Qu. Max. 150.0 160.0 162.0 162.4 169.5 171.0 N 10 样本个数 s 7.5748 标准差 se 2.3954 样本平均数的标准误 cv 4.66429 变异系数 gmean 162.2386 几何平均数 hmean 162.0748 调和平均数 结果表明，“金皇后”玉米果穗的重量平均为162.4 克，标准差为7.57 克，变异系数 为4.7% 。 二、单个平均数的假设测验 （一）命令 在R 中，有一个t.test 命令用来完成未知总体方差的单个平均数的假设测验和置信区 间的估计，简介如下： 用法：t.test(x, y = NULL, alternative = c(two.sided, less, greater), mu = 0, paired = FALSE, var.equal = FALSE, conf.level = 0.95, ...) 功能：对列变量x 执行t 测验。 说明：x 为列变量；alternative 可以指明是两尾测验、或大于或小于的一尾测验；mu 是总体平均数，默认值为0 ；paired 指明是否为成对数据 （但y 赋值时）；var.equal 指明 两个样本的方差是否相同；conf.level 指明置信区间的概率值。 R 没有提供正太分布的z 测验函数，我们在应用时自己编写一个。 （二）实例 〖例2〗某蔗糖自动打包机在正常工作状态时的每包蔗糖重量符合N(100 ，2)分布。 某日抽样调查10 包，得结果为100.5、99.8、101