为什么SAS不为混合模型中的固定效应提供p值?

时间:2017-08-07 11:38:50

标签: sas mixed-models p-value

我使用混合模型来查看变量的效果。我使用的代码是:

proc mixed data=pb2;
class treat_a treat_b hoknr_ day;
model conc=treat_a|treat_b hoknr_/outp=residuals1 residual;
repeated day/subject=hoknr_(treat_a treat_b)type=vc;
run;

对于treat_a,treat_b或treat_a | treat_b,结果没有p值,但它对hoknr_没有。我排除了重复的陈述,我简化了模型,我改变了类,但我仍然没有得到所有固定效果的p值。我之前使用过这个模型并且它有效,现在我将它安装到这个数据集中,但它并没有完全发挥作用。

EDIT1 类型3固定效果测试表如下所示: Type 3 Table。 治疗可能是不可估计的(treat_a是yes或no,同样对于treat_b)。我已将yes / no更改为0或1,未更改Type 3表。之前我曾经用过治疗方法来表达不会产生这样的表格。

Edit2 当解决方案添加到模型语句时,结果如下:Solution for Fixed effects

这个模型有什么问题,它没有显示所有固定效果的p值?

2 个答案:

答案 0 :(得分:2)

您可能需要使用htype=选项指定所需的测试类型。这听起来像是有人最初没有对函数进行编程的那些过程之一,并且它在开发后期是一种事后的想法(与showpvalues中的proc glmselect选项不同;直到今天我认为这是回归过程中最奇怪的选项。)

https://support.sas.com/documentation/cdl/en/statug/63033/HTML/default/viewer.htm#statug_mixed_sect015.htm#statug.mixed.mixedmodelhtype

  

类型3固定效果测试

     

您可以使用MODEL语句中的HTYPE =选项来获取表格   另外,还有1类(顺序)测试和2类(调整)测试   或者代替第3类(部分)测试表。

     

ODS表名称对于Type 1测试是“Tests1”,对于Type 2测试是“Tests2”,对于Type 3测试是“Tests3”。

或者,可能是某些固定效果不可估量。

答案 1 :(得分:1)

导致结果为0的原因是您的Treat_a和Treat_b是类别变量。 Treat_a = 1和Treat_b = 1是参考水平。因此,您在解决方案表中缺少p值。在您的互动条件中,如果条件包括treat_a = 1或treat_b = 1,则不会有p值。