双因素方差分析不满足正态性或方差齐次时,如何做非参数检验替代方法?

原创 数据小兵  2020-11-08 10:56  阅读 15,230 次
数据小兵成长记

双因素方差分析不满足正态性或方差齐次时,如何做非参数检验替代方法?

双因素方差分析对数据适用有明确要求,当数据不满足正态性、方差齐次条件,如果因素水平组合下个案数≥2,可用基于秩次的 Scheirer-Ray-Hare 检验实现,Scheirer-Ray-Hare 检验是 Kruskal-Wallis H 检验的扩展,可作为双因素方差分析非参数检验替代解决方案。

可惜啊,咱们常用的SPSS统计软件并没有提供Scheirer-Ray-Hare 检验,所以SPSS是无法直接实现的。好消息是,一些科研工作者们已经开始尝试利用SPSS完成该检验方法。

申希平,祁海萍等学者利用SPSS软件的个案排秩、计算变量、方差分析等菜单功能实现了Scheirer-Ray-Hare 检验,他们的期刊论文标题为《两因素非参数方差分析在 SPSS 中的实现》,于2013年12月发表于《中国卫生统计》。

01 下载该论文

在本公众号后台回复关键词【srh】即可获取论文下载链接,并阅读参考,解决SPSS双因素方差分析非参数方法的实现问题。

下面,小兵给大家演示一下,如何使用R语言实现Scheirer-Ray-Hare 检验,很方便,对R较为熟悉的读者,可以看一看练习一下。

02 R案例演示

用R语言实现Scheirer-Ray-Hare检验需要使用rcompanion包,首先安装这个包。

 install.packages("rcompanion")

案例数据采用申希平,祁海萍等《两因素非参数方差分析在 SPSS 中的实现》论文数据,录入excel或txt文本,用R语言读取。

因变量:研究指标ml
处理因子:group三水平含对照、细胞及膜片组
时间因子:time二水平含3周和6周

本例暂不讨论正态性,方差齐次等问题,主要用来演示Scheirer-Ray-Hare检验。

两因素方差分析非参数检验方法Scheirer-Ray-Hare检验,具体由rcompanion包scheirerRayHare()函数来实现。scheirerRayHare()函数的一些用法与R默认的aov函数类似,可以参考aov()函数了解参数要求。

具体语法为:

scheirerRayHare(y~A*B)

A*B表示考察A和B主效应及交互项AB。

本案例Scheirer-Ray-Hare()函数语法

library(rcompanion)

scheirerRayHare(ml~group*time,data = srh)

03 结果解读

R原汁原味的结果如下:

经整理为三线表结果如下:

不同处理组间指标ml有统计学差异(H=6.088,P=0.048<0.05),不同时间周期下指标ml有统计学差异(H=21.389,P<0.001)。处理和时间周的交互作用对指标ml无影响(H=0.027,P=0.987>0.05)。

Scheirer-Ray-Hare检验,为偏态分布、方差不齐或小样本计量资料的两因素或多因素设计提供了有用的统计分析方法。

申希平,祁海萍等《两因素非参数方差分析在 SPSS 中的实现》这篇论文,介绍了如何利用SPSS实现Scheirer-Ray-Hare检验,以及其后续的两两比较问题,文章中有具体的操作方法和计算过程。如果有读者遇到此类问题,且就此问题对SPSS情有独钟的话,可以按这篇文章介绍的计算方法来实现Scheirer-Ray-Hare检验。

本文完
文/图=数据小兵

 

◢ 为你推荐以下文章 ◣

单因素方差分析的非参数方法
曼惠特尼独立样本非参数检验
交叉设计及SPSS多因素方差分析
SPSS正交试验设计及其方差分析
大家别笑啊,我也在偷偷学Python
《SPSS统计分析:快速入门与实践提高》

视频地址:

https://study.163.com/course/introduction/1003945001.htm?share=1&shareId=1149679450

本文地址:http://www.datasoldier.net/archives/2572
版权声明:本文为原创文章,版权归 数据小兵 所有,欢迎分享本文,转载请保留出处!
视频课程《SPSS统计分析:从入门到实践提高》
欢迎订阅SPSS训练营微信公众号

评论已关闭!