lcmm (version 1.4-1)

plot.predict: Plot of class-specific marginal predictions for the longitudinal outcome

Description

This function provides the class-specific predicted trajectories stemmed from a hlme, lcmm or Jointlcmm object.

Usage

## S3 method for class 'hlme':
plot.predict(x,newdata,var.time,legend.loc="topright",ylim=NULL,\dots)
## S3 method for class 'lcmm':
plot.predict(x,newdata,var.time,legend.loc="topright",ylim=NULL,\dots)
## S3 method for class 'Jointlcmm':
plot.predict(x,newdata,var.time,legend.loc="topright",ylim=NULL,\dots)

Arguments

x

an object inheriting from classes hlme, lcmm or Jointlcmm representing respectively a fitted latent class linear mixed-effects model, a general latent class mixed model or a joint latent class mixed model.

newdata

data frame containing the data from which predictions are computed. Data frame should include all the covariates listed in x$Xnames (in the same order). Names of data frame should be exactly x$Xnames.

var.time

A character string containing the name of the variable that corresponds to time in the data frame (x axis in the plot).

legend.loc

keyword for the position of the legend from the list "bottomright", "bottom", "bottomleft", "left", "topleft","top", "topright", "right" and

ylim

optional numeric vector of length 2, giving the y coordinate range.

...

further arguments to be passed to or from other methods. They are ignored in this function.

Value

• Returns a plot

See Also

hlme, lcmm, Jointlcmm

Examples

################# Prediction from linear latent class model
data(data_hlme)
## fitted model
m<-lcmm(Y~Time+X1+X1_time,mixture=~Time,random=~Time,classmb=~X2+X3,
subject='ID',ng=2,data=data_hlme)
## newdata for predictions plot
newdata<-data.frame(intercept=rep(1,100),Time=seq(0,5,length=100),
X1=rep(0,100),X1_time=rep(0,100),X2=rep(0,100),X3=rep(0,100))
plot.predict(m,newdata,"Time","right")
## data from the first subject for predictions plot
firstdata<-cbind(intercept=rep(1,3),data_hlme[1:3,c("Time","X1",
"X1_time","X2","X3")])
plot.predict(m,firstdata,"Time","right")


################# Prediction from a joint latent class model
data(data_Jointlcmm)
## fitted model - see help of Jointlcmm function for details on the model
m3 <- Jointlcmm(fixed= Ydep1~Time+X1+X1_Time,mixture=~Time,random=~Time,
classmb=~X3,subject='ID',survival = Surv(Tevent,Event)~ mixture(X1)+X2
,hazard="5-equi-splines",hazardtype="PH",ng=3,data=data_Jointlcmm)
# class-specific predicted trajectories (with characteristics of subject ID=193)
data <- data_Jointlcmm[data_Jointlcmm$ID==193,]
plot.predict(m3,var.time="Time",newdata=data)