用mixed模型解决单因素方差分析中方差不齐时的随机模拟及SAS实现

摘要

目的 用mixed模型处理完全随机设计的三组资料方差不齐时的Ⅰ类错误率和对各组总体均值的估计,并与Kruskal-Wallis秩和检验和ANOVA方法比较检验的把握度.方法 利用SAS9.3产生模拟数据,并分别用proc mixed过程、proc nparl way过程和proc glm过程进行mixed模型参数估计、Kruskal-Wallis秩和检验和one-way ANOVA分析.模拟试验重复1000次.结果 方差不齐时,对不均衡资料,mixed法的检验把握度要高于Kruskal-Wallis秩和检验和one-way ANOVA法.对均衡资料,三种方法的检验效能接近,但后两者的Power受方差不齐的影响较大.结论 研究设计时尽量使各组样本均衡或接近则三种方法的Ⅰ类错误率最小,Power达到最大.对不均衡资料方差不齐时,mixed法的检验把握度要高于Kruskal-Wallis秩和检验和one-way ANOVA法,且其Ⅰ类错误率能得到控制.