lcmm (version 1.5.7)

plot.baselinerisk: Plot of the class-specific baseline risk functions estimated from a joint latent class mixed model

Description

Displays the estimated baseline risk functions for the time-to-event of interest in each latent class

Usage

## S3 method for class 'Jointlcmm':
plot.baselinerisk(x,legend.loc="topleft",ylim=NULL,...)

Arguments

x

an object inheriting from class Jointlcmm, representing a fitted joint latent class mixed model

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.

...

other parameters to be passed through to plotting functions

See Also

Jointlcmm

Examples

data(data_Jointlcmm)
#### estimation with 3 latent classes (ng=3) - see Jointlcmm 
#### help for details on the model
m3 <- Jointlcmm(fixed= Ydep1~Time*X1,mixture=~Time,random=~Time,
classmb=~X3,subject='ID',survival = Surv(Tevent,Event)~ X1+mixture(X2),
hazard="3-quant-splines",hazardtype="PH",ng=3,data=data_Jointlcmm,
B=c( 0.7667 ,  0.4020,  -0.8243,  -0.2726,   0.0000 ,  0.0000 ,  0.0000 ,
 0.3020, -0.6212,   2.6247 ,  5.3139,  -0.0255,   1.3595,   0.8172, -11.6867,
  10.1668,10.2355 , 11.5137,  -2.6209,  -0.4328,  -0.6062 ,  1.4718 ,
 -0.0378  , 0.8505, 0.0366,   0.2634  , 1.4981))


# Class-specific predicted baseline risk & survival functions in the 
# 3-class model retained (for the reference value of the covariates) 
plot.baselinerisk(m3,bty="l")
plot.baselinerisk(m3,ylim=c(0,5),bty="l")
plot.survival(m3,bty="l")