用R做两因素多水平重复测量方差分析（球对称检验不通过版）

步骤包括：数据导入——边列表转换——异常值检验——正态性检验——ANOVA+球对称检验——校正——组间对比——输出

library(readxl)
data<-read_excel("对比测试.xlsx")

#数据格式转变
library(tidyverse)
data2<-gather(data,"time","电量",-c("group","id"))

data2$id<-factor(data2$id)
data2$time<-factor(data2$time)
data2$group<-factor(data2$group)

data3<-subset(data2,is.na(data2$电量)==F) #去除空行

#异常值检验
boxplot(data3$电量~data3$time,boxwex = 0.25,
        col=c("blue","red","green")) #颜色数量与group数量对应

#正态性检验
library(tidyverse)
data3 %>% group_by(time, group) %>% summarise(p_value = shapiro.test(电量)$p.value)
#p>0.05表示不能拒绝正态性假设

#球检验和anova：方法一
library(ez)
ezANOVA(data = data3, dv = .(电量), wid = .(id), 
        within = .(time), between = .(group), 
        detailed = TRUE, type = 2)

#球检验和anova：方法二
library(rstatix)
anova_test(data = data3,dv = 电量,wid = id,
           within = time, between=group)

#组间对比
library(PMCMRplus)
summary(lsdTest(电量 ~ group, data = data3))