0

How to calculate the standard error for variance in the multilevel model?

m5 <- lmer(peakflow ~ pm25 + gender +age + (1| participant), data = mydata, REML = F)
summary(m5)

Random effects:
 Groups      Name        Variance Std.Dev.
 participant (Intercept) 9410.6   97.01   
 Residual                 632.8   25.16   
Number of obs: 12256, groups:  participant, 97

There are just standard deviations of the variances. How can I know the standard error?

Thanks!

Crystal Lee
  • 1
  • The [VCA package](https://cran.r-project.org/web/packages/VCA/index.html) provides the variance-covariance matrix of the variance components. – Stéphane Laurent Feb 03 '22 at 12:18
  • and https://stats.stackexchange.com/questions/161225/estimates-of-the-variance-of-the-variance-component-of-a-mixed-effects-model# – user20650 Feb 03 '22 at 12:19

0 Answers0