lqmm {BeQut}

lqmm fits linear quantile mixed model

Description

Function using 'JAGS' software to estimate the linear quantile mixed model assuming asymmetric Laplace distribution for residual error.

Usage

lqmm(
  formFixed,
  formRandom,
  formGroup,
  data,
  tau,
  RE_ind = FALSE,
  n.chains = 3,
  n.iter = 10000,
  n.burnin = 5000,
  n.thin = 1,
  n.adapt = NULL,
  precision = 10,
  save_jagsUI = TRUE,
  parallel = FALSE
)

Arguments

formFixed

formula for fixed part of longitudinal submodel with response variable

formRandom

formula for random part of longitudinal submodel without response variable

formGroup

formula specifying the cluster variable (e.g. = ~ subject)

data

dataset of observed variables

tau

the quantile(s) to be estimated. This must be a number between 0 and 1, otherwise the execution is stopped. If more than one quantile is specified, rounding off to the 4th decimal must give non–duplicated values of tau, otherwise the execution is stopped.

RE_ind

Boolean denoting if the random effects are assumed independent ; default is FALSE

n.chains

the number of parallel chains for the model; default is 1.

n.iter

integer specifying the total number of iterations; default is 10000

n.burnin

integer specifying how many of n.iter to discard as burn-in ; default is 5000

n.thin

integer specifying the thinning of the chains; default is 1

n.adapt

integer specifying the number of iterations to use for adaptation; default is NULL

precision

variance by default for vague prior distribution

save_jagsUI

If TRUE (by default), the output of jagsUI package is returned by the function. Warning, if TRUE, the output can be large.

parallel

see jagsUI::jags() function

Value

A Blqmm object is a list with the following elements:

mean

list of posterior mean for each parameter

median

list of posterior median for each parameter

modes

list of posterior mode for each parameter

StErr

list of standard error for each parameter

StDev

list of standard deviation for each parameter

ICs

list of the credibility interval at 0.95 for each parameters excepted for covariance parameters in covariance matrix of random effects. Otherwise, use save_jagsUI=TRUE to have the associated quantiles.

data

data included in argument

sims.list

list of the MCMC chains of the parameters and random effects

control

list of arguments giving details about the estimation

random_effect

list for each quantile including both posterior mean and posterior standard deviation of subject-specific random effects

out_jagsUI

only if save_jagsUI=TRUE in argument: list including posterior mean, median, quantiles (2.5%, 25%, 50%, 75%, 97.5%), standart deviation for each parameter and each random effect. Moreover, this list also returns the MCMC draws, the Gelman and Rubin diagnostics (see output of jagsUI objects)

Author(s)

Antoine Barbieri

References

Marco Geraci and Matteo Bottai (2014). Linear quantile mixed models. Statistics and Computing, 24(3):461-479. doi: 10.1007/s11222-013-9381-9.

Examples

#---- Use dataLong dataset
data(dataLong)

#---- Fit regression model for the first quartile
lqmm_075 <- lqmm(formFixed = y ~ visit,
                 formRandom = ~ visit,
                 formGroup = ~ ID,
                 data = dataLong,
                 tau = 0.75,
                 n.iter = 10000,
                 n.burnin = 1000)

#---- Get the posterior means
lqmm_075$mean

#---- Visualize the trace for beta parameters
jagsUI::traceplot(lqmm_075$out_jagsUI, parameters = "beta")

#---- Summary of output
summary(lqmm_075)