R语言做重复测量的球形检验？

球形检验 (Mauchly's Test of Sphericity) 是前提条件，该检验综合评价组间以及各个时间点测量指标是否满足方差齐性的，P值大于等于0.1表示满足球形检验。

  R语言实例  

既然球形检验那么重要，那么《R语言实战》为何不进行介绍？难道很难实现？且看如下简洁代码。以下实例仅作为方法的讲解，临床价值参考需谨慎。

1. #加载R包，仅为了其中的数据OBrienKaiser

2. library(car)

3. data=OBrienKaiser[,-c(1:2)]

5. #将data从数据框转换为矩阵，data的格式为16行×15列

6. #15列分别是3个组的5个时间点：

7. #A1,A2,A3,A4,A5,B1,B2,B3,B4,B5,...
   data=as.matrix(data)

9. #进行多元线性回归 (多个因变量)

10. mlmfit=lm(data~1)

12. #定义数据的列结构

13. group=factor(rep(c("A","B","C"), c(5, 5, 5)))

14. time=ordered(rep(1:5,3))

15. idata=data.frame(group,time)

17. #球形检验

18. <