实验训练2

实验训练2

一、实验内容

数据：设甲、乙两种零件彼此可以代替，但乙种零件比甲种零件制造简单，造价也低，经过试验获得它们的抗拉强度分别为（单位：kg/cm2）

甲：88，87，92，90，91; 乙：89，89，90，84，88.

假设两种零件的抗拉强度都服从正态分布。

1. 求甲种零件的抗拉强度的期望和乙种零件的抗拉强度的方差的置信区间；

2. 求两种零件的抗拉强度的均值差的置信区间（假设两个总体方差相等）；

3. 求两种零件的抗拉强度的方差比的置信区间；

4. 检验甲种零件的抗拉强度的均值是否为$90kg/c{m}^2$；

5. 检验甲种零件的抗拉强度的均值与乙种零件的抗拉强度的均值有无显著差异。

（提示：两个正态总体样本容量一样，检验均值的差异性，可以两总体取差转化为一个正态总体，方差未知，均值是否为0的假设检验问题。）

置信度都取0.95，显著性水平都取为0.05.

二、问题求解

导入数据：

> a = c(88, 87, 92, 90, 91)
> b = c(89, 89, 90, 84, 88)

（ 一 ）求甲种零件的抗拉强度的期望和乙种零件的抗拉强度的方差的置信区间

1. 求甲种零件的抗拉强度的期望的置信区间

总体均值的区间估计（方差未知）采用t. test()函数来计算，调用格式如下：

t.test(x, y = NULL, alternative=c("two sided", "less", "greater"), mu = 0,paired = TRUE, var.equal = FALSE,
       conf.level = 0.95,...)

其中，x是样本数据，当y=NULL（系统默认值）时，表示对单个正态总体进行区间估计；alternative表示所求置信区间的类型，默认为双尾检验；mu表示均值，均值未知时不需要赋值；paired表示是否是成对检验；var.equal表示双样本的方差是否相等

其中， conf.level用来确定要计算的置信区间的置信度.

求甲种零件的抗拉强度的期望的置信区间的代码如下：

> t.test(a, y=NULL, mu = 0, paired = FALSE, var.equal = FALSE, conf.level = 0.95)

输出结果如图：

由输出结果可知，在置信度为0.95时，甲种零件的抗拉强度的期望的置信区间为 [ 87.02523，92.17477 ] 。

2. 求乙种零件的抗拉强度的方差的置信区间

根据均值已知/未知情况，用卡方分布估计总体方差置信区间，实际情况中均值多为未知。自行编写函数，可用以下代码实现：

var.conf.int=function(x,mu=Inf,alpha){
     n=length(x)
     if(mu<Inf){
        s2=sum((x-mu)^2)/n
        df=n}
     else{
        s2=var(x)
        df=n-1}
     c(df*s2/qchisq(1-alpha/2,df),df*s2/qchisq(alpha/2,df))}

求乙种零件的抗拉强度的方差的置信区间的代码如下：

> var.conf.int = function(x, mu = Inf, alpha){
+  n = length(x)
+  if(mu < Inf){
+    s2 = sum((x - mu)^2)/n
+    df = n
+  }
+  else{
+    s2 = var(x)
+    df = n - 1
+  }
+  c(df*s2/qchisq(1-alpha/2,df),df*s2/qchisq(alpha/2,df))
+ }
> var.conf.int(b, mu = Inf, alpha = 0.05)

输出结果如图：

由输出结果可知，在置信度为0.95时，乙种零件的抗拉强度的方差的置信区间为 [ 1.974283，45.415271 ] 。

（ 二 ）求两种零件的抗拉强度的均值差的置信区间

当两总体方差未知但相等或两总体方差未知且不等，求均值差${\mu }_1-{\mu }_2$的置信区间时采用t.test()函数来计算，调用格式如下：

t.test(x, y , alternative=c("two sided","less","greater"), mu = 0,paired = FALSE, var.equal = FALSE, 
       conf.level = 0.95,...)

其中，当y的取值不为NULL（系统默认值）时，t.test()函数是对两组样数据x和y进行分析，paired=FALSE(系统默认值)表示两组样本不是配对数据，当取值为TRUE时，表示两组样本为配对数据. Var.equal是逻辑变量，当取值为FALSE(系统默认值)时，表示两样本对应总体方差不等，取值为TRUE时，表示总体方差相等.

求两种零件的抗拉强度的均值差的置信区间的代码如下：

> t.test(a, b, mu = 0, paired = FALSE, var.equal = TRUE, conf.level = 0.95)

输出结果如下：

由输出结果可知，当置信度为0.95时，两种零件的抗拉强度的均值差的置信区间为 [ -1.628406，4.828406 ] 。

（ 三 ）求两种零件的抗拉强度的方差比的置信区间

求两总体方差比的置信区间时，采用var.test()函数计算，调用格式如下：

var.test(x, y, ratio = 1, alternative = c("two.sided","less","greater"), conf.level = 0.95,...)

其中，x,y为样本数据；ratio为原假设的方差比值；alternative设置检验类型为双尾或是单尾；conf.level为置信水平.

求两种零件的抗拉强度的方差比的置信区间的代码如下：

> var.test(a, b, ratio = 1, conf.level = 0.95)

输出结果如下：

由输出结果可知，当置信度为0.95时，两种零件的抗拉强度的方差比的置信区间为 [ 0.08140098，7.50899609 ] 。

（ 四 ）检验甲种零件的抗拉强度的均值是否为$90kg/c{m}^2$

我们采用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是样本数据，当y=NULL（系统默认值）时，表示对单个正态总体进行假设检验；

alternative=c(“two.sided”, “less”, “greater”)规定了假设检验备择假设，two.sided对应着备择假设$H_1:\mu \ne {\mu }_0$，less对应着备择假设$H_1:\mu <{\mu }_0$，greater对应着备择假设$H_1:\mu >{\mu }_0$，mu为比较的均值，系统默认为0.

检验甲种零件的抗拉强度的均值是否为$90kg/c{m}^2$时，备择假设为$H_1:\mu \ne 90$，代码如下：

> t.test(a, alternative="two.sided", mu=90)

输出结果如下：

由输出结果可知，显著概率$p=0.6885$，在显著性水平$\alpha =0.05$时，$p>\alpha$，接受原假设，则甲种零件的抗拉强度的均值为$90kg/c{m}^2$。

（ 五 ）检验甲种零件的抗拉强度的均值与乙种零件的抗拉强度的均值有无显著差异

根据提示，两个正态总体样本容量一样，检验均值的差异性，可以两总体取差转化为一个正态总体，方差未知，均值是否为0的假设检验问题，则处理数据将两个数据向量相减：

> c = a - b

接下来检验c的均值是否为0，原假设$H_0:\mu =0$，备择假设$H_1:\mu \ne 0$，检验代码如下：

> t.test(c, alternative="two.sided", mu=0)

输出结果如下：

由输出结果可知，显著概率$p=0.3274$，在显著性水平$\alpha =0.05$时，$p>\alpha$，接受原假设，则甲种零件的抗拉强度的均值与乙种零件的抗拉强度的均值无显著差异。

三、问题总结

1. 甲种零件的抗拉强度的期望的置信区间为 [ 87.02523，92.17477 ] ，乙种零件的抗拉强度的方差的置信区间为 [ 1.974283，45.415271 ] ；
2. 两种零件的抗拉强度的均值差的置信区间为 [ -1.628406，4.828406 ] ；
3. 两种零件的抗拉强度的方差比的置信区间为 [ 0.08140098，7.50899609 ] ；
4. 甲种零件的抗拉强度的均值为$90kg/c{m}^2$；
   87.02523，92.17477 ] ，乙种零件的抗拉强度的方差的置信区间为 [ 1.974283，45.415271 ] ；
5. 两种零件的抗拉强度的均值差的置信区间为 [ -1.628406，4.828406 ] ；
6. 两种零件的抗拉强度的方差比的置信区间为 [ 0.08140098，7.50899609 ] ；
7. 甲种零件的抗拉强度的均值为$90kg/c{m}^2$；
8. 甲种零件的抗拉强度的均值与乙种零件的抗拉强度的均值无显著差异。