plot.SMNclmm {skewlmm}

Plot a smn.clmm object

Description

Fitted values versus residuals plot. Censored values are imputed using their conditional expectation from the fitting algorithm.

Usage

  ## S3 method for class 'SMNclmm'
plot(x, level = "conditional", useweight = TRUE,
                         alpha = 0.3, ...)

Arguments

x

An object inheriting from class SMNclmm, representing a fitted scale mixture of normal censored linear mixed model.

level

Level of residual that should be used. For details see residuals.SMNclmm. Default is "conditional".

useweight

A logical value indicating if the estimated weights should be used as color in the plot.

alpha

Transparency parameter to be used (0<alpha<1). Meaningless if useweight = TRUE.

...

Additional arguments.

Value

A ggplot object.

Author(s)

Fernanda L. Schumacher, Larissa A. Matos, Victor H. Lachos and Katherine L. Valeriano

See Also

ggplot, smn.clmm, fitted.SMNclmm, residuals.SMNclmm

Examples

nj1 = 5; m = 30
time = rep(1:nj1, times=m)
groups = as.factor(rep(1:m, each=nj1))
dat1 = rsmsn.clmm(time, groups, cbind(1,time), rep(1,m*nj1), sigma2=0.7,
                  D=0.5*diag(1), beta=c(1,2), depStruct="UNC", distr="t", nu=4)
# Estimation
fm1 = smn.clmm(dat1, formFixed=y~x, groupVar="ind", depStruct="UNC", ci="ci",
               lcl="lcl", ucl="ucl", distr="t", control=lmmControl(max.iter=30))
plot(fm1)