[R-sig-ME] Fwd: Could the random effect at the level of each observation be a trap?

Billy billy.requena at gmail.com
Fri Dec 17 22:59:51 CET 2010


---------- Forwarded message ----------
From: Billy <billy.requena at gmail.com>
Date: Fri, Dec 17, 2010 at 3:57 PM
Subject: Re: [R-sig-ME] Could the random effect at the level of each
observation be a trap?
To: Drew Tyre <atyre2 at unl.edu>


Hi Drew!

Well, below I'm sending the summary of two corresponding models (just
for simplicity, because I have at least 5 models including onlye the
'MaleID' random effect and other 5 including 'MaleID' + 'resid'
randiom effects):

m1 <- glmer (y~ x + (1|MaleID), family=binomial)
m1.1 <- glmer (y~ x + (1|MaleID) + (1|resid), family=binomial)

> summary(m1)
Generalized linear mixed model fit by the Laplace approximation
Formula: y ~ x + (1 | MaleID)
 AIC  BIC logLik deviance
 1012 1027 -502.9     1006
Random effects:
 Groups Name        Variance   Std.Dev.
 MaleID (Intercept) 6.9038e-11 8.309e-06
Number of obs: 1045, groups: MaleID, 464

Fixed effects:
           Estimate Std. Error z value Pr(>|z|)
(Intercept)   -8.273      2.036  -4.063 4.84e-05 ***
x                   1.192      0.354   3.366 0.000762 ***
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Correlation of Fixed Effects:
  (Intr)
x  -0.999

> summary(m1.1)
Generalized linear mixed model fit by the Laplace approximation
Formula: y~ x + (1 | MaleID) + (1 | resid)
  AIC   BIC logLik deviance
 619.4 639.2 -305.7    611.4
Random effects:
 Groups Name        Variance Std.Dev.
 resid  (Intercept) 2089.1   45.707
 MaleID (Intercept)    0.0    0.000
Number of obs: 1045, groups: resid, 1045; MaleID, 464

Fixed effects:
           Estimate Std. Error z value Pr(>|z|)
(Intercept)  -19.049     85.552  -0.223    0.824
x                    1.281     14.844   0.086    0.931

Correlation of Fixed Effects:
  (Intr)
x  -0.999


We can see that the variance associated to the 'resid' random variable
is really huge, but I don't really know if that's a problem or it's
just expected due the nature of this individual observation-based
variable.
Furthermore, the correlation between fixed effects seems pretty
strong, but again I don't know if it's a problem.
Does anyone have an insight to help me?

Thanks again

Billy


--
Gustavo Requena
PhD student - Laboratory of Arthropod Behavior and Evolution
Universidade de São Paulo
Correspondence adress:
a/c Glauco Machado
Departamento de Ecologia - IBUSP
Rua do Matão - Travessa 14 no 321 Cidade Universitária, São Paulo - SP, Brasil
CEP 05508-900
Phone number: 55 11 3091-7488

http://ecologia.ib.usp.br/opilio/gustavo.html



-- 
Gustavo Requena
PhD student - Laboratory of Arthropod Behavior and Evolution
Universidade de São Paulo
Correspondence adress:
a/c Glauco Machado
Departamento de Ecologia - IBUSP
Rua do Matão - Travessa 14 no 321 Cidade Universitária, São Paulo - SP, Brasil
CEP 05508-900
Phone number: 55 11 3091-7488

http://ecologia.ib.usp.br/opilio/gustavo.html




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