TOP

R语言各种假设检验实例整理(常用)(一)
2017-10-10 12:08:58 】 浏览:8650
Tags:语言 各种 假设 检验 实例 整理 常用

版权声明:本文为博主原创文章,转载请注明出处

一、正态分布参数检验

   例1. 某种原件的寿命X(以小时计)服从正态分布N(μ, σ)其中μ, σ2均未知。现测得16只元件的寿命如下:

                 159 280 101 212 224 379 179 264
                  222 362 168 250 149 260 485 170

    问是否有理由认为元件的平均寿命大于255小时?

    解:按题意,需检验

                     H0: μ ≤ 225     H1: μ >  225

      此问题属于单边检验问题

      可以使用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,y是又数据构成e向量,(如果只提供x,则作单个正态总体的均值检验,如果提供x,y则作两个总体的均值检验),alternative表示被则假设,two.sided(缺省),双边检验(H1:μ≠H0),less表示单边检验(H1:μ<μ0),greater表示单边检验(H1:μ>μ0),mu表示原假设μ0,conf.level置信水平,即1-α,通常是0.95,var.equal是逻辑变量,var.equal=TRUE表示两样品方差相同,var.equal=FALSE(缺省)表示两样本方差不同。

      R代码:

       X<-c(159, 280, 101, 212, 224, 379, 179, 264,
              222, 362, 168, 250, 149, 260, 485, 170)
       t.test(X,alternative = "greater",mu=225)

       结果:

       

         可见P值为0.257 > 0.05 ,不能拒绝原假设,接受H0,即平均寿命不大于225小时。

 

    例2.在平炉上进行的一项试验以确定改变操作方法的建议是否会增加刚的得率,试验时在同一个平炉上进行的,每炼一炉刚时除操作方法外,其它条件都尽可能做到相同,先用标准方法炼一炉,然后用新方法炼一炉,以后交替进行,各炼了10炉,其得率分别为

                           标准方法    78.1 72.4 76.2 74.3 77.4 78.4 76.0 75.5 76.7 77.3
                           新方法        79.1 81.0 77.3 79.1 80.0 79.1 79.1 77.3 80.2 82.1

设这两个样本相互独立,且分别来自正态总体N(μ1, σ2)和N(μ2, σ2),其中μ1,μ2和σ2未知。问新的操作能否提高得率?(取α=0.05)

    解1:根据题意,需要假设

              H0: μ1 ≥  μ2     H1: μ1 <  μ2  

    这里假定σ12=σ22=σ2,因此选择t.test,var.equal=TRUE

    R代码:

    X<-c(78.1,72.4,76.2,74.3,77.4,78.4,76.0,75.5,76.7,77.3)
    Y<-c(79.1,81.0,77.3,79.1,80.0,79.1,79.1,77.3,80.2,82.1)
    t.test(X,Y,var.equal = TRUE,alternative = "less")

   结果:

 

    可见P值<0.05,接受备择假设,即新的操作能够提高得率。

    解2:

    因为数据是成对出现的,所以采用成对数据t检验比上述的双样本均值检验更准确。所谓成对t检验就是Zi=Xi-Yi,再对Z进行单样本均值检验

    R代码:

    X<-c(78.1,72.
R语言各种假设检验实例整理(常用)(一) https://www.cppentry.com/bencandy.php?fid=91&id=130237

首页 上一页 1 2 3 4 5 6 7 下一页 尾页 1/8/8
】【打印繁体】【投稿】【收藏】 【推荐】【举报】【评论】 【关闭】 【返回顶部
上一篇R语言实现对基因组SNV进行注释 下一篇时间序列分析算法【R详解】