Calculate expected information matrix — Info_Expected

Calculates the expected information matrix from a fitted linear mixed effects model with AR(1) correlation structure in the level-1 errors.

Info_Expected_lmeAR1(m_fit)

Arguments

m_fit: Fitted model of class lme, with AR(1) correlation structure at level 1.

Value

Expected Information matrix corresponding to variance components of m_fit.

Examples

data(Laski)
Laski_RML <- lme(fixed = outcome ~ treatment, 
                 random = ~ 1 | case, 
                 correlation = corAR1(0, ~ time | case), 
                 data = Laski)
Info_Expected_lmeAR1(Laski_RML)
#>               [,1]          [,2]         [,3]
#> [1,]  1.601686e-03  -0.210945594 5.092243e-06
#> [2,] -2.109456e-01 127.633048102 1.874239e-03
#> [3,]  5.092243e-06   0.001874239 4.954985e-05