
双因素方差分析对数据适用有明确要求,当数据不满足正态性、方差齐次条件,如果因素水平组合下个案数≥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
数据小兵坚持写博客已经12年坚持写微信公号文章6年坚持更新SPSS视频课程2年坚持一对一答疑讨论2年
绝对超值:一对一答疑 欢迎加入SPSS视频课程竭诚服务

◢点【