plot.SMNCens function

Plot a smn.clmm object

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

## 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.

Returns

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)