## Orthodont data
data(Orthodont)
# Random intercept model
fit <- lqmm(distance ~ age, random = ~ 1, group = Subject,
tau = 0.5, data = Orthodont)
fit.boot <- boot.lqmm(fit)
# extract fixed effects
B <- extractBoot.lqmm(fit.boot, "fixed")
# covariance matrix estimated fixed parameters
cov(B)Run the code above in your browser using DataLab