我使用混合模型来查看变量的效果。我使用的代码是:
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值?
答案 0 :(得分:2)
您可能需要使用htype=
选项指定所需的测试类型。这听起来像是有人最初没有对函数进行编程的那些过程之一,并且它在开发后期是一种事后的想法(与showpvalues
中的proc glmselect
选项不同;直到今天我认为这是回归过程中最奇怪的选项。)
类型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值。