
   Hello, dear list:

   I am trying to extract the random factor estimates from a mixed model I 
implemented. I've made some trials with the function "ranef" but I'm not sure 
I'm doing it alright, so I'd much appreciate some help.

   From the summary output I get this:

   
Random effects:
 Groups   Name        Variance Std.Dev.
 year     (Intercept) 0.060513 0.24599 
 Residual             0.092965 0.30490 
Number of obs: 31, groups: year, 2

   I've seen in some papers that they present the random effect results in the 
form "estimate ± SD", and I´m not sure about where I can take that information 
from. Using "ranef" I only got this:

   
mincoefvar<-lmer(inslarge~tlengthcoef+brivencoef+brithrcoef+condcoef+sex+tlengthcoef:brivencoef+

brivencoef:brithrcoef+brivencoef:condcoef+brithrcoef:condcoef+brithrcoef:sex+brivencoef:brithrcoef:condcoef+(1|year),na.action=na.omit)

> ranef(mincoefvar)
$year
     (Intercept)
2009   0.2436456
2010  -0.2436456

   I'm pretty sure I'm doing something wrong and not getting the results I'm 
looking for, but I can't see where the mistake is. 



   Thank you very much for your help.

   Iker Vaquero-Alba
   Centre for Ecology and Conservation
   Daphne du Maurier Building
   University of Exeter, Cornwall Campus
   Treliever Road
   TR10 9EZ Penryn
   U.K.

   http://biosciences.exeter.ac.uk/cec/staff/postgradresearch/ikervaquero-alba/

	[[alternative HTML version deleted]]


