# NOT RUN {
fit <- ddhazard(
Surv(time, status == 2) ~ log(bili), pbc, id = pbc$id, max_T = 3600,
Q_0 = diag(1, 2), Q = diag(1e-4, 2), by = 50,
control = list(method = "GMA"))
predict(fit, type = "response", new_data =
data.frame(time = 0, status = 2, bili = 3))
predict(fit, type = "term", new_data =
data.frame(time = 0, status = 2, bili = 3))
# }
Run the code above in your browser using DataLab