双因素方差分析对数据适用有明确要求，当数据不满足正态性、方差齐次条件，如果因素水平组合下个案数≥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/AMOS/R语言/Process/Minitab等，在职场中也会用到统计方法，用起来的时候，第一印象是：好简单点两下就出结果，第二印象是：好像有点难哦，第三印象是： 没有一个结果我认识~~
不要听人家说简单易用就真的以为简单，难的是怎么把统计分析方法应用到你的研究问题上，出来的结果怎么解释，结论怎么描述。
你可能恰好是一个没人管的研究生，科研无人带，数据分析无人带，只听导师、师兄师姐咔咔咔讲给你怎么做，具体一做的时候，可能会处处碰壁。
所以，你需要一个略懂统计方法的人来帮你，帮你回答一些疑问，帮你梳理研究分析的思路，帮你解释统计分析结果，帮你撰写呈现统计分析结论。
欢迎来找小兵。
我推出纯答疑、咨询、交流讨论问题的服务，专门和大家讨论统计分析软件、统计分析方法在应用时遇到的一些问题，给大家答疑，帮忙整理分析思路，对一些统计概念进行通俗讲解，帮你解释分析结果，整理统计表格和描述统计结论。
付费的，大家有需要的可以扫描上方二维码，加我具体聊。