代码之家  ›  专栏  ›  技术社区  ›  Prometheus

sas-proc混合自由度误差

  •  0
  • Prometheus  · 技术社区  · 6 年前

    我有这样一个数据集:

    enter image description here

    因子A是一个介于受试者之间的因子(有两个水平-高和低)。 因子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。我有预感,这就是错误的征兆。但我一直不明白为什么会这样。任何线索将不胜感激。谢谢。

    1 回复  |  直到 6 年前
        1
  •  1
  •   John Doe    6 年前

    根据SAS文件:

    估计的G矩阵不是正定消息通常指示随机语句上的一个或多个方差分量被估计为零,并且可以/应该从模型中移除。

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

    至于自由度,知道违约 ddfm 对于带有随机语句的模型 contain 需要大样本,我建议你尝试其他方法 KR SAT 是的。

    一旦你的模型产生了结果,你就可以根据你的偏好,根据你拥有的数据修改模型。

    推荐文章