首 页IT知识库翔宇问吧收藏内容
当前位置:算搜IT乐园IT知识库R语言

R语言中进行多重均值检验的方法

减小字体增大字体作者:suansuo.cn  来源:本站原创  发布时间:2019-03-28 19:58:01

本站在《使用R语言进行单因素方差分析》一文中介绍了在R语言中使用aov函数进行方差分析的方法,并给出具体的例子讲解了该函数的使用方法及其正态性检验和方差齐性检验的方法。在实际中,在进行单因素方差分析时,可能检验结果拒绝了零假设,即因素A的r个水平效应具有显著的差异,也就是说r个均值之间有显著区别。但这并不能说明r个均值之间都会存在这个显著差异,只能说明至少2个总体的均值有显著差异。那么,到底是哪些水平之间存在差异呢?我们需要对每两个水平间的均值做一一比较,即多重比较。

同时比较任意两个水平均值间有无明显差异的问题称为多重比较。通常采用多重T检验方法进行多重比较,这种方法本质上就是针对每组数据进行T检验,只不过估计方差时利用的是全体数据,因而自由度变大。具体地说,要比较第i个总体与第j个总体的均值是否相同,即检验:(薛毅,陈立萍.R语言在统计中的应用)

H0: μi = μj,H1:μi ≠ μj,i ≠ j,i,j=1,2,...,r

在R中,使用pairwise.t.test()函数可以帮助我们完成多重T检验,该函数的原型如下:

pairwise.t.test(x, g, p.adjust.method = p.adjust.methods, pool.sd = !paired,
paired = FALSE, alternative = c("two.sided", "less", "greater"), ...)

其中,主要参数的含义如下:

x 是一个向量,代表响应变量(因变量);

g 代表响应向量的因子或向量;

p.adjust.method 调整p值所使用的方法;

paired 逻辑变量,表示是否做成对的t检验,默认为FALSE,此变量和pool.sd不能同时为TRUE;

alternative 字符串,表示备择假设。取值为two.sided(双侧检验),less(单侧小于),greater(单侧大于),默认为two.sided.

下面用一个例子来说明本函数的具体使用方法:

为了减小失业保险支出,促进就业,政府试图为失业者提供再就业奖励:如果失业者可以在限定的时间内重新就业,他将可以获得一定数额的奖金,这个政策是否有效呢?下图给出了不同奖励下再就业的天数。

R语言多重T检验的例子

(本例子来源于网络)

编写程序如下:

#定义数据
x<-c(92,100,85,88,89,90,94,80,78,
86,108,93,88,89,75,78,72,79,
96,92,90,77,79,71,82,75,81,
78,75,76,87,73,83,82,68,72)
A<-factor(rep(1:4,c(9,9,9,9)))

#定义数据框
dat<-data.frame(x,A)

#方差分析
aovResult<-aov(x~A,data=dat)

#提取方差分析表
summary(aovResult)

#多重比较
pairwise.t.test(x,A)

方差分析结果如下图所示:

 多重比表T检验

由上图方差分析的结果来看,p值=0.0433<0.05,则可以认为有充分理由认为水平间的效用是不同的。

那么,不同产生在哪些水平间呢,下面是多重比较的结果:

多重比较的结果

 从上图可以看出1和4之间的p值为0.04,则可以认为有充足的理由认为水平1和水平4之间差异显著,即无奖金组和高奖金组的再就业时间有显著性差异。而其它组之间没有差异。

上图还给出了p值调整所用的方法:holm修正的方法。

因为在多重T检验时,当因素的水平较多,而检验又是同时进行的,则多次重复使用T检验会增大犯第一类错误的概率,得到的“有显著差异”的结论不可靠。为了克服这一缺陷,统计学家提出了多种方法来调整P值,这里使用的是holm调整方法。

-------------------------------------------

本文为本站原创,如需转载请注明出处:算搜IT乐园(www.suansuo.cn)-R语言中进行多重均值检验的方法(http://www.suansuo.cn/article/R/2019/6402.html)

微信搜索“优雅的代码”关注本站的公众号,或直接使用微信扫描下面二维码关注本站公众号,以获取最新内容。

个人成长离不开各位的关注,你的关注就是我继续前行的动力。

知识评论评论内容只代表网友观点,与本站立场无关!

   评论摘要(共 0 条,得分 0 分,平均 0 分) 查看完整评论
愿您的精彩评论引起共鸣,带来思考和价值。
用户名:查看更多评论
分 值:100分 90分 80分 70分 60分 40分 20分
内 容:
验证码:
关于本站 | 网站帮助 | 广告合作 | 网站声明 | 友情连接 | 网站地图
本站部分内容来自互联网,如有侵权,请来信告之,谢谢!
Copyright © 2007-2021 suansuo.cn. All Rights Reserved.
名词解释产业结构,中国目前的产业支柱是哪些行业?可预见的未来会实现哪些产业升级与转型?_策鸾股票网公司最佳负债比率,新公司资产负债表怎么填写?_厕锫股票网大盘指数股票行情走势预测,大盘走势5月21日回顾及5月24日预测_侧殆股票网烟的利润是多少比例,一盒香烟的成本大约是多少?_参雄股票网中国居民负债率是多少,中国地方政府的巨额债务是怎么形成的?_册篡股票网金融股包括哪些板块,牛市发展路径及板块轮动规律_侧戳股票网原始股东减持要求,低税率股权转让、自然人减持、税收优化、财政扶持_策焯股票网底部打压吸筹盘中表现,庄家吸筹分布全程曝光,死记“三大吸筹手法”,及时出手,随时暴涨_策纂股票网h股什么意思与a股的区别,a股和h股同时上市,原股东的股票在哪儿交易?_册谆股票网怎样判断股票趋势的拐点,判断拐点的简单方法如何判断基金走势,如何判断某只基金相对底部长期均线金叉,KDJ指标,30分钟金钗,60分钟死叉120分钟金叉,日线死叉,以哪个周期为准?月数据的趋势图怎么做,数据的可视化——学习制作图表和初步分析第四换手率0.55怎么样,阐述a股和港股的三大区别低位高换手率放量上涨说明什么?,涨停板换手高的股票预示什么?次新股换手率50以上要出吗,历史上有没有日换手率达到或者接近100%的股票?换手率到多少易开板,股票换手率高说明什么?换手率百分之50,在股市里,一只股票换手率高达45%,意味着什么?值得所有人反省简述市盈率法计算股票价格的方法,「市盈率(P/E)」是如何定义与计算的,有何参考价值?查询指数基金市盈率的网站,股票、基金、财报、行业数据查阅网站市盈率很低的股票可以买吗,那些4倍市盈率的股票靠谱吗?市盈率动静多少才合理,中国股市:献给每一个决心支持自己的股东PB是什么意思?,有没有老粉科普一下兴家的pb呀,从哪些方面可以了解到,十分感谢!?通达信筹码分布指标,通达信的筹码。它是怎么知道当前位置的筹码,是从哪个部位的筹码换手来的?筹码峰什么形态是买点,“筹码集中度”寻找牛股诀窍,学会赚钱一步到位,至少8个涨停筹码双峰形态是什么意思,国内“筹码分布”第一人,玩转股市的奥秘在于这两种筹码形态筹码分布指标,428家公司透露最新股东户数 306家筹码趋向集中筹码是什么,股票筹码90%集中度是怎么理解的?铜陵有色金属集团股份有限公司官网,请问鹤山天山金属材料制品有限公司是国企吗?什么叫流通股票,股票的底层逻辑汇总(一)