Question

我试图在SAS中进行简单的单向随机效应ANOVA。我想知道人口方差是否显着不同于零。

在UCLA's idre site上，他们声明如下使用PROC MIXED：

proc mixed data = in.hsb12 covtest noclprint;
   class school;
   model mathach = / solution;
   random intercept / subject = school;
run;

根据我之前使用PROC MIXED的经验，这对我来说很有意义。

然而，在Murray Logan撰写的Biostatistical Design and Analysis Using R中，他说，对于单向ANOVA，固定和随机效应没有区分，并且（在R中）进行“标准”单向ANOVA，即使他正在测试方差，而不是手段。我发现在SAS中，他的R程序相当于使用以下任何一种方法：

PROC ANOVA
PROC GLM（与ANOVA相同，但用GLM代替ANOVA）
PROC GLM with RANDOM statement

上述三个模型中的 p - 值是相同的，但与UCLA使用的PROC MIXED模型不同。对于我的数据，它是 p = 0.2508和 p = 0.3138的差异。虽然在这种情况下结论没有改变，但我对这种差异感到不舒服。

任何人都可以就哪一个更合适以及为什么会出现这种差异提出建议？

Answer 1

对于您的模型，PROC ANOVA和PROC MIXED之间的差异仅仅是由于数值噪声（PROC MIXED的REML估计）。但是，问题中提到的p值对应于不同的测试。为了使用PROC MIXED中的COVTEST输出获得F值，您需要重新计算MS_groups，同时考虑到不相等的样本大小（手动，如http://bio.classes.ucsc.edu/bio286/MIcksBookPDFs/QK08.PDF的第223页所述，或者仅使用PROC MIXED与与PROC ANOVA中相同的固定模型规范）。除了SAS手册之外，本文（http://isites.harvard.edu/fs/docs/icb.topic1140782.files/S98.pdf）提供了一些PROC MIXED的使用示例。

SAS中的单向随机效应ANOVA：PROC GLM还是MIXED？

1 个答案: