09-23-2014 11:38 AM
Please help me sort out this output of Negative Binomial regression in PROC GENMOD. The model includes a binary factor, Factor_B. There are a few p-values associated with Factor_B that I expect to be consistent (see the attachment):
1) In “Analysis of Maximum Likelihood estimates” Wald p-value for the corresponding regression coefficient (0.0226)
2) LR test, Type 3, (0.8667)
3) Wald test, Type 3, (0.8666)
4) Contrast “B1 vs B2”, LR (0.8667)
5) Contrast “B1 vs B2”, Wald (0.8666)
As you can see, 1) seems to be a problem. Why is it so different? Note also that the problem showed up only after I fitted the model with Factor_A*Factor_B interaction. In the additive version, everything was consistent.
09-23-2014 12:59 PM
09-23-2014 11:56 AM
I think the problem is related to the interaction term "Factor_A*Factor_B". It is not meaningfull to test away the maineffect of Factor_B effect when the interactionterm between Factor_A and Factor_B is included in the model.
Will the problem persist if you remove the interactionterm?
09-23-2014 12:59 PM
09-23-2014 02:52 PM
Thank you. The way I see it is as follows:
1) The Wald/LR Type 3 tests and Wald/LR Contrasts aim at testing whether the mean response is the same in the corresponding cells, B1 and B2 in this case. LSMEANS statement is of a similar nature.
All these tools are dependent on GLM parameterization so much so that LSMEANS provides no output when reference parameterization is used.
2) Because there is interaction, the mean response depends not only on the main effects of B but also on the A*B interaction effects. On the other hand, the main table provides Wald p-values just for the main effects of B, i.e. it tests a different hypothesis.
3) It may be possible to remove the discrepancy by switching to reference parametrization, but then the interpretation of 1) as a test for the equality of mean responses will be lost.
09-23-2014 04:02 PM
You don't want to switch to reference parameterization, unless you want this for other purposes. Interpretation of the global tests becomes really strained when there is an interaction in the model.
10-03-2014 01:52 PM
I think 'strained' is putting it mildly. In my opinion, there was a reason Goodnight and Sarle developed the non-full rank parameterization: in agricultural field studies, interactions are the rule rather than the exception, and you just can't get good global tests with full rank parameterizations (see Speed and Hocking for examples).
Need further help from the community? Please ask a new question.