Standard Error of variance component from the output of lmer .

Cynderella · Jul 29, 2015

I want to extract the "standard error" of variance component from the output of "lmer" .

Code:

    library(lme4)
    model <- lmer(Reaction ~ Days + (1|Subject), sleepstudy)

The following produces estimates of variance component :

Code:

    s2 <- VarCorr(model)$Subject[1]

It is **NOT** the standard error of the variance. And I want the standard error . How can I have it ?

Many Thanks! Regards.