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