1、两组杀虫剂分别处理蚜虫,5小时后记录它们的死亡数,如图所示,问两组死亡数是否有差异?由于没有假设样本是正态分布的,采用自助法进行计算
2、 输入数据##输入两组数据x1<-c(35,37,43,41,37,32,45,36,44,31,38,31)x2<-c(32,33,36,35,泌驾台佐36,31,34,32,39,38,39,34,33)
3、 皱诣愚继样本随机有放回的抽样到自助样本,计算差值统计量,重复1000次##sample(x1,10,replace = TRUE) 为抽取十个样本形成自助样本##抽取一次形成自助样本,并做差值得到统计量mean1 = mean(sample(x1,10,replace = TRUE))-mean(sample(x2,10,replace = TRUE))##通过for循环,得到第2到1000次的差值统计量,不断的传入mean1中for(i in 2:1000){ a = mean(sample(x1,10,replace = TRUE))-mean(sample(x2,10,replace = TRUE)) mean1 = c(mean1,a)}
4、 将结果排序sort(mean1)
5、 找出97.5%和2.5%的分位点,即初始位置和最末位置的25个数,它们就限制了95%的置信区间quantile(mean1,0.975)quantile(mean1,0.025)
6、 判断0是否在里面,若在,接受原假设,没有显著性差异;若不在,拒绝原假设,有显著性差异。结果显示,0在置信区间内,接受原假设,没有显著性差异