设为首页 加入收藏

TOP

R语言各种假设检验实例整理(常用)(六)
2017-10-10 12:08:58 】 浏览:8057
Tags:语言 各种 假设 检验 实例 整理 常用
试分析两种饲料养猪有无显著差异。

                   

     R语言代码:

     x<-scan()
     25 30 28 23 27 35 30 28 32 29 30 30 31 16

     y<-scan()
     19 32 21 19 25 31 31 26 30 25 28 31 25 25

     binom.test(sum(x<y), length(x))

   

sum(x < y)表示样品X小于样品Y的个数。计算出P值大于0.05,无法拒绝原假设,可以认为两种饲料养猪无显著差异。

 

例15.某饮料店为了解顾客对饮料的爱好情感,进一步改进他们的工作,对顾客喜欢咖啡还是喜欢奶茶,或者两者同样爱好进行了调查,顾客喜欢咖啡超过奶茶用正号表示,喜欢奶茶超过咖啡用负号表示,两者同样爱好用0表示。现将调查的结果列在下表中。试分析顾客是喜欢咖啡还是喜欢奶茶。

                               

解:根据题意可检验如下假设:

    H0:顾客喜欢咖啡等于喜欢奶茶;      H1:顾客喜欢咖啡超过奶茶。

    以上资料中有以人(即6号顾客)表示对咖啡和奶茶有同样爱好,用0表示,因此在样本容量中不加计算,所以实际上N=12.如果H0假设为真,那么符合p为1/2的二项分布,如果H1为真,那么顾客喜欢奶茶的人数小于理论值,al="l",因此用R软件进行计算,显著性水平取α = 0.10,

    R语言代码:

    binom.test(3,12,p=1/2, al="l", conf.level = 0.90)

 

 

    可见P值 < 0.1 ,置信区间也不包括0.5,因此拒绝原假设人口喜欢咖啡的人超过喜欢奶茶的人。

    在符号检验法中,只计算符号的个数,而不考虑每个符号差所包含的绝对值的大小,因此常常使用弥补了这个缺点的wilcoxon符号秩检验。

3.3.4.符号秩检验

    例16.假定某电池厂宣称该厂生产的某种型号电池寿命的中位数为140安培小时。为了检验改厂生产的电池是否符合其规定的标准,现从新近生产的一批电池中抽取了随即样本,并对这20个电池的寿命进行了测试,其结果如下(单位:安培小时):

                 137.0 140.0 138.3 139.0 144.3 139.1 141.7 137.3 133.5 138.2

                 141.1 139.2 136.5 136.5 135.6 138.0 140.9 140.6 136.3 134.1

试用Wilcoxon符号秩检验分析该厂生产的电池是否符合其标准。

    解:根据题意假设:

                 H0:电池中位数M≥ 140安培小时;

                 H1:电池中位数<140安培小时。

          在R语言中进行符号秩检验可以使用wilcox.test( )

          wilcox.test(x, y = NULL,
              alternative = c("two.sided", "less", "greater"),
              mu = 0, paired = FALSE, exact = NULL, correct = TRUE,
              conf.int = FALSE, conf.level = 0.95, ...)

           其中x,y是观察数据构成的数据向量。alternative是备择假设,有单侧检验和双侧检验,mu待检参数,如中位数M0.paired是逻辑变量,说明变量x,y是否为成对数据。exact是逻辑变量,说明是否精确计算P值,当样本量较小时,此参数起作用,当样本两较大时,软件采用正态分布近似计算P值。correct是逻辑变量,说明是否对P值的计算采用连续性修正,相同秩次较多时,统计量要校正。conf.int是逻辑变量,说明是否给出相应的置信区间。

         R语言代码:

         X<-scan()
    

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

最新文章

热门文章

Hot 文章

Python

C 语言

C++基础

大数据基础

linux编程基础

C/C++面试题目