# NOT RUN {
# Fit msm model:
Q <- rbind(c(0,0.01,0.01), c(0,0,0.01), c(0,0,0))
model <- msm(state~age, subject = id, data = electData,
center = FALSE, qmatrix = Q, deathexact = TRUE,
covariates = ~age+x)
# Estimate function of life expectancies:
sddata <- electData[electData$bsline == 1,]
LEs <- elect(x = model, b.covariates = list(age = 0, x = 1),
statedistdata = sddata, h = 0.5, age.max = 50, S = 25)
plusmin(LEs, index = c(2,4), func = "plus")
# }
Run the code above in your browser using DataLab