[R-sig-ME] Binary response ordering
John Haart
another83 at me.com
Wed Aug 4 10:54:15 CEST 2010
Dear List,
I have a quick question regarding the setup of my data for analysis with a glmm. I hope this is the appropriate list, i apologise if it is not.
I have a response variable, TRUE or FALSE. I have coded this as 0 = False and 1 = TRUE in excel.
I have 3 categorical factors with C,D and E
I then read in the data frame and run the model as follows-
lmer(trueorfalse~1+(1|A/B) + C + D+ E ,family=binomial)
And this is the output
Generalized linear mixed model fit by the Laplace approximation
Formula: threatornot ~ 1 + (1 | A/B) + C + D+ E ,family=binomial)
AIC BIC logLik deviance
1410 1450 -696.8 1394
Random effects:
Groups Name Variance Std.Dev.
family:order (Intercept) 6.7869e-01 8.2382e-01
order (Intercept) 7.8204e-11 8.8433e-06
Number of obs: 1116, groups: A:B, 43; B, 9
Fixed effects:
Estimate Std. Error z value Pr(>|z|)
(Intercept) 0.11281 0.42232 0.267 0.7894
C1 -0.02414 0.19964 -0.121 0.9038
D2 -0.16482 0.38602 -0.427 0.6694
E2 0.95381 0.54316 1.756 0.0791 .
E3 0.75733 0.87275 0.868 0.3855
E4 0.03044 0.47328 0.064 0.9487
What i am unsure about is the inference, if a term is significant does this relate to TRUE or FALSE?
I.E E2 has a p value of 0.079, does this 0.079 relate to the probability of it resulting in a true or false response? Does it matter how i code the input i.e FALSE = 1, TRUE =2 for instance?
Maybe i am reading the output wrong?
Thanks
John
More information about the R-sig-mixed-models
mailing list