基于R语言的区间估计自编函数，部分

作业1：为比较两个小麦品种的产量，选择18块条件相似的试验田，采用相同的耕作方法做实验，结果播种甲品种的8块实验田的单位面积产量和播种乙品种的10块实验田的单位面积产量分别为：

甲品种 628 583 510 554 612 523 530 615

乙品种 535 433 398 470 567 480 498 560 503 426

假定每个品种的单位面积产量均服从正态分布，试求这两个品种平均面积产量差的置信区间（取α=0.05）

要求使用两种方式进行求解，自行编写函数 & t.test()。

作业2：甲、乙两台机床分别加工某种轴承，轴承的直径都服从正态分布，从各自加工的轴承中分别抽取若干个轴承测其直径，结果如下表所示。试求两台机床加工的轴承直径的方差比σ12/σ22的0.95置信区间。

甲(20.5,19.8,19.7,20.4,20.1,20.0,19.0,19.9)
乙(20.7,19.8,19.5,20.8,20.4,19.6,20.2)
要求使用两种方式进行求解，自行编写函数 & var.test()。

程序

two.sample.ci<-function(x,y,conf.level=0.95)

{

  options(digits=4)

  m = length(x); n =length(y)

  xbar=mean(x)-mean(y); alpha=1-conf.level;

  vp=(var(x)/m+var(y)/n)^2 *(m-1)*(n-1)/((var(x)/m)^2*(n-1)+(var(y)/n)^2*(m-1))

  tstar= qt(alpha/2,vp)*(var(x)/m+var(y)/n)^(1/2)

  xbar +c(-tstar, +tstar)

}

x1 = c(628,583,510,554,612,523,530,615)

x2 = c(535,433,398,470,567,480,498,560,503,426)

 two.sample.ci(x1,x2)

t.test(x1,x2,var.equal = F)

作业2：

two.sample.ci<-function(x,y,conf.level=0.95){

  options(digits=4)

  m = length(x); n = length(y)

   alpha = 1 - conf.level

 F1<-qf(1-alpha/2,m-1,n-1)

 F2<-qf(alpha/2,m-1,n-1)

di=var(x)/var(y)

di/c(F2,F1)

}

two.sample.ci(y1,y2)

var.test(y1,y2)