predint: Predict confidence interval from an 'lqmm' object

In gforge/lqmm_gforge: Linear Quantile Mixed Models

Description

A confidence interval predictor for an [lqmm] model fit ([lqm] not implemented) at the moment.

Usage

predint(
  object,
  level = 0,
  alpha = 0.05,
  R = 50,
  seed = round(runif(1, 1, 10000)),
  newdata
)

Arguments

object	The model fit
level	an optional integer vector giving the level of grouping to be used in obtaining the predictions.
alpha	1-alpha is the confidence level.
R	number of bootstrap replications.
seed	optional random number generator seed.
newdata	Defaults to the original data used for generating the model. Currently only implemented for 'level = 0' predictions.

Value

A data frame or a list of data frames for predint.lqmm containing predictions, lower and upper bounds of prediction intervals, and standard errors.

Author(s)

Marco Geraci

References

Geraci M and Bottai M (2014). Linear quantile mixed models. Statistics and Computing, 24(3), 461–479.

See Also

lqmm, ranef.lqmm, coef.lqmm

Examples

## Orthodont data
data(Orthodont)

# Random intercept model
fitOi.lqmm <- lqmm(distance ~ age, random = ~ 1, group = Subject,
	tau = c(0.1,0.5,0.9), data = Orthodont)

# 95% confidence intervals
predint(fitOi.lqmm, level = 0, alpha = 0.05)

gforge/lqmm_gforge documentation built on Dec. 20, 2021, 10:42 a.m.