## ----include = FALSE----------------------------------------------------------
knitr::opts_chunk$set(collapse = TRUE, comment = "#>", eval = FALSE)

## -----------------------------------------------------------------------------
# library(bqmm)
# 
# fit <- bqmm(distance ~ age + (1 | Subject),
#             data = nlme::Orthodont,
#             tau  = 0.5)
# 
# summary(fit)

## -----------------------------------------------------------------------------
# fit3 <- bqmm(distance ~ age + (1 | Subject),
#              data = nlme::Orthodont,
#              tau  = c(0.1, 0.5, 0.9))
# 
# coef(fit3)      # tau-by-coefficient matrix
# plot(fit3)      # coefficient-versus-tau paths

## -----------------------------------------------------------------------------
# vcov(fit, adjusted = TRUE)
# vcov(fit, adjusted = FALSE)   # naive posterior covariance