[R] Very Simple Regression Question

i_finlayson37 ian.finlayson at blueyonder.co.uk
Thu Feb 28 04:57:04 CET 2008


I've just been running my first ever regression and I'm using R for an lmer.
I've created two models, a null one (part ~ 1 + (1 | id) + (1 | word)) and
another with a predictor I want to test (part ~ 1 + conf + (1 | id) + (1 |
word)). I've compared them both and the model with a predictor is
significantly better, but I can't see from the data which direction the
prediction goes (i.e. are conf and part inversely correlated or not). Could
somebody tell me where I can find this information? Here's my data, thanks
in advance:

Null Model:

Generalized linear mixed model fit using Laplace 
Formula: part ~ 1 + (1 | id) + (1 | word) 
   Data: align 
 Family: binomial(logit link)
   AIC BIC logLik deviance
 537.9 551 -265.9    531.9
Random effects:
 Groups Name        Variance Std.Dev.
 word   (Intercept) 2.309122 1.51958 
 id     (Intercept) 0.055252 0.23506 
number of obs: 601, groups: word, 32; id, 20

Estimated scale (compare to  1 )  0.8671394 

Fixed effects:
            Estimate Std. Error z value Pr(>|z|)    
(Intercept)  -1.8920     0.3072  -6.158 7.37e-10 ***
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Model with predictors:

Generalized linear mixed model fit using Laplace 
Formula: part ~ 1 + conf + (1 | id) + (1 | word) 
   Data: align 
 Family: binomial(logit link)
   AIC   BIC logLik deviance
 490.1 507.7 -241.1    482.1
Random effects:
 Groups Name        Variance Std.Dev.
 word   (Intercept) 2.977549 1.72556 
 id     (Intercept) 0.097771 0.31268 
number of obs: 601, groups: word, 32; id, 20

Estimated scale (compare to  1 )  0.9142627 

Fixed effects:
            Estimate Std. Error z value Pr(>|z|)    
(Intercept)  -3.1001     0.4014  -7.724 1.13e-14 ***
conf          1.7941     0.2734   6.563 5.29e-11 ***
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 

Correlation of Fixed Effects:
     (Intr)
conf -0.505

ANOVA:

Data: align
Models:
align.nul: part ~ 1 + (1 | id) + (1 | word)
align.conf: part ~ 1 + conf + (1 | id) + (1 | word)
           Df     AIC     BIC  logLik  Chisq Chi Df Pr(>Chisq)    
align.nul   3  537.85  551.05 -265.93                             
align.conf  4  490.14  507.73 -241.07 49.713      1  1.780e-12 ***
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

-- 
View this message in context: http://www.nabble.com/Very-Simple-Regression-Question-tp15729195p15729195.html
Sent from the R help mailing list archive at Nabble.com.



More information about the R-help mailing list