0

我有一个这种格式的数据集:

在此处输入图像描述

因素 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。我有一种预感,这就是错误的症状。但我一直无法弄清楚为什么会这样。任何线索将不胜感激。谢谢。

4

1 回答 1

1

根据 SAS 文档:

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

我建议你从一个更简单的随机语句开始,也许只是用拦截并逐渐添加效果(在考虑交互之前首先主要)。

至于自由度,知道ddfm带有随机语句的模型的默认值是contain需要大样本量,我建议您尝试其他方法,例如KRor SAT

一旦您的模型产生结果,您就可以继续并根据您拥有的数据根据​​您的偏好修改模型。

于 2018-06-26T08:23:32.897 回答