R语言做Scheirer–Ray–Hare检验

毕设有组重复测量数据非正态分布，分析要用到R语言，下面这个算是小白自学的步骤记录！

R Handbook: Scheirer–Ray–Hare Test

Clear examples in R. Scheirer Ray Hare test

https://rcompanion.org/handbook/F_14.html

这网站里的内容更专业全面（奈何俺不太会操作）

1. 安装Rstudio有很多教程，需要注意的是安装程序包，做Scheirer–Ray–Hare检验需要rcompanion和FSA两个数据包，直接在Rstudio的界面敲进去

if(!require(rcompanion)){install.packages("rcompanion")}
if(!require(FSA)){install.packages("FSA")}

2. 查找默认的工作目录（这里也可以用setwd()设置），等会要把准备好的txt文档放在这个目录里

getwd()

3. 把SPSS里的数据搞成长条，做成txt格式文件，命名最好用英文

不同的文件格式好像都需要下载程序包，导入数据前可以先查查都需要什么。

4. 导入自己做好的txt文件，公式里的file就是你自己设置的文件名，文件第一行为表头的就是TRUE，否则就是FALSE，（简写T和F好像也行）

mydata <- read.table("file.txt",header = TURE)
mydata

5. 这步应该是为了强制把导入的mydata转换成R语言可以识别的string类型？

str(mydata)

成功后会显示

6. 接下来就可以直接进行分析了，我要做的是组间、组内和交互效应，这里的y代表因变量，A和B分别是两种主效应。

library(rcompanion)
scheirerRayHare(y ~ A*B,data = mydata)