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

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