[R-sig-ME] Correlation of -1.00 in a mixed model

Manuel Spínola mspinola10 at gmail.com
Sun Mar 28 15:28:55 CEST 2010


Dear list members,

A few days ago somebody posted a problem like the want that I have now, 
I got a correlation value of - 1.00.  I would like to have some help in 
the context of my model fitting.
My model is a random intercept random slope model (repeated measure 
design).  Variable "time" has been centered.
What will be an alternative model?

 > modelo7 = lmer(ipa ~ time + (time | localidad), data=ipa)
 > summary(modelo7)
Linear mixed model fit by REML
Formula: ipa ~ time + (time | localidad)
   Data: ipa
  AIC  BIC logLik deviance REMLdev
 1917 1937 -952.4     1913    1905
Random effects:
 Groups      Name           Variance     Std.Dev.    Corr  
 localidad  (Intercept)        2232.31    47.247
                   time                545.97     23.366   -1.000
 Residual                           450.92     21.235
Number of obs: 196, groups: localidad, 49

Fixed effects:
            Estimate Std. Error t value
(Intercept)   42.852      6.918   6.194
time         -19.746      3.603  -5.480

Correlation of Fixed Effects:
     (Intr)
time -0.904

Thank you very much in advance.
Best,

Manuel

-- 
Manuel Spínola, Ph.D.
Instituto Internacional en Conservación y Manejo de Vida Silvestre
Universidad Nacional
Apartado 1350-3000
Heredia
COSTA RICA
mspinola at una.ac.cr
mspinola10 at gmail.com
Teléfono: (506) 2277-3598
Fax: (506) 2237-7036




More information about the R-sig-mixed-models mailing list