Hi Petar,

last I tried this, the estimated variance of the random effects is
(optionally) stored as an attribute.  So, something like this should work
...

rfg <- ranef(my.lmer, postVar=TRUE)

my.se = sqrt(as.numeric(attributes(rfg$group)$postVar))

Cheers

Andrew

On Sat, Nov 20, 2010 at 8:46 PM, Petar Milin <pmilin@ff.uns.ac.rs> wrote:

> Hello all!
> This is my third attempt (first at the general R-list, and one here five
> days or so ago), and I apologize for being tedious, but this started to
> bother me and my students are making my life a misery with this:
> How can one get upper and lower limits of a prediction interval for
> random-effect levels; the exact exact values, numbers? They are shown on
> caterpillar plot using ranef() with
> argument postVar=TRUE, but I would like to know them. A while ago, some
> discussions were opened on "Confidence Intervals for Random Effect
> BLUP's", but the answer was never clear:
> http://www.mail-archive.com/r-help@r-project.org/msg04820.html
>
> Thanks!
> PM
>
> _______________________________________________
> R-sig-mixed-models@r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
>



-- 
Andrew Robinson
Program Manager, ACERA
Senior Lecturer in Applied Statistics                      Tel:
+61-3-8344-6410
Department of Mathematics and Statistics            Fax: +61-3-8344 4599
University of Melbourne, VIC 3010 Australia
Email: a.robinson@ms.unimelb.edu.au    Website: http://www.ms.unimelb.edu.au

	[[alternative HTML version deleted]]


