################################################### ### Paquid study ################################################### library("SmoothHazard") data("Paq1000", package = "SmoothHazard") head(round(Paq1000, 1)) fit.weib <- idm(formula01 = Hist(time = list(l, r), event = dementia, entry = e) ~ certif + gender, formula02 = Hist(time = t, event = death, entry = e) ~ gender, formula12 = ~ 1, data = Paq1000) fit.weib par(mgp = c(4, 1, 0), mar = c(5, 5, 5, 5)) plot(fit.weib, conf.int = TRUE, lwd = 3, citype = "shadow", xlim = c(65, 100), axis2.las = 2, axis1.at = seq(65, 100, 5), xlab = "Age (years)") fit.splines <- idm(formula01 = Hist(time = list(l, r), event = dementia, entry = e) ~ certif + gender, formula02 = Hist(time = t, event = death, entry = e) ~ gender, formula12 = ~ 1, method = "Splines", data = Paq1000) fit.splines par(mgp = c(4, 1, 0), mar = c(5, 5, 5, 5)) plot(fit.splines, conf.int = TRUE, lwd = 3, citype = "shadow", xlim = c(65, 100), axis2.las = 2, axis1.at = seq(65, 100, 5), xlab = "Age (years)") pred <- predict(fit.weib, s = 70, t = 80, newdata = data.frame(certif = 1, gender = 1)) pred LE.weib <- predict(fit.weib, s = 80, newdata = data.frame(certif = 1, gender = 0), lifeExpect = TRUE) LE.weib LE.splines <- predict(fit.splines, s = 80, newdata = data.frame(certif = 1, gender = 0), lifeExpect = TRUE, CI = FALSE) LE.splines