SAS proc混合自由度错误

Question

我有一个采用以下格式的数据集：

因子A在主题因子之间（具有2个级别-高和低）。 因子B在主题因子之内（具有3个级别-高，中和低）。

我想运行一个具有嵌套随机效应因子的混合模型。

我正在使用的代码是：

proc mixed data=data.mydata;
class FactorA FactorB;
model DV = FactorA|FactorB;
random FactorB(FactorA) FactorB*FactorA(FactorA);
lsmeans FactorA|FactorB;
run;

日志指出：估计的G矩阵不是正定的。 我也没有得到任何p值（仅显示“。”）。

此外，在输出表中，我看到DF =0。我预感这是错误的征兆。但是我一直无法弄清楚为什么会这样。任何线索将不胜感激。谢谢。

Answer 1

根据SAS文档：

  

“估计的G矩阵不是正定”消息通常表明RANDOM语句上的一个或多个方差分量被估计为零，并且可以/应该从模型中删除。

我建议您从一个简单的随机语句开始，也许只是从拦截开始，然后逐渐增加效果（在考虑交互作用之前首先要考虑）。

关于自由度，知道带有随机语句的模型的默认ddfm是contain，这需要大量样本，我建议您尝试使用其他方法，例如KR或SAT。

一旦模型产生结果，您就可以继续根据自己的偏好根据所拥有的数据修改模型。