[R-sig-ME] Testing for differences in random effects between two groups

Stephen T stwebvanuatu at yahoo.com.au
Mon May 27 14:13:33 CEST 2013


Hello, 
 
I am trying to apply linear models to some observational data. 
The data are not balanced. 

Basically, I am wish to describe sexual and individual signatures in calls of a bird species. I have many call recordings and measured several acoustic variables from each recording.

Here's an analysis for one variable: 
 
> lmer(MH11~SEX+(1|BIRD)+(1|NIGHT), data=calls) 
Linear mixed model fit by REML  
Formula: MH11 ~ SEX + (1 | BIRD) + (1 | NIGHT)  
   Data: calls  
  AIC  BIC logLik deviance REMLdev 
 3662 3681  -1826     3663    3652 
Random effects: 
 Groups   Name        Variance Std.Dev. 
 NIGHT    (Intercept)  652.77  25.549   
 BIRD     (Intercept) 1083.67  32.919   
 Residual              966.50  31.089   
Number of obs: 356, groups: NIGHT, 138; BIRD, 57 
 
Fixed effects: 
            Estimate Std. Error t value 
(Intercept)  445.059      8.212   54.19 
SEXM          94.779     10.588    8.95 
 
Correlation of Fixed Effects: 
     (Intr) 
SEXM -0.776 
 
SEX is a fixed effect (Female/Male). BIRD (individuals) and NIGHT are random effects. 
The nesting is SEX/BIRD/NIGHT/CALL.

There is a clear difference between-SEXes (fixed effect), i.e. SEXF and SEX
M have different voices.

Secondly, there is some interesting variation between-BIRDs. Individual BIRDs have different voices.


What if the BIRD random effect was not the same in both SEXF and SEXM groups? That would confound the model. Since there is a fixed effect (sexual signature), perhaps the random effect (individual signature) is different between groups SEXF and SEXM as well (unexpected, but should be tested).

I'm not sure how to test for differences between groups in the random effects. One possibility is to split SEXF and SEXM and run separate lmer models. How would I the compare the results? Another is to compare the distributions of random effects between SEXF and SEXM (a two-sample test?). 
 
I would appreciate some advice on how to proceed?

Stephen from Australia.



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