# \donttest{
library(dynamichazard)
set.seed(56219373)
fit <- ddhazard(
Surv(time, status == 2) ~ log(bili), pbc, id = pbc$id, max_T = 3000,
Q_0 = diag(1, 2), Q = diag(1e-4, 2), by = 100,
control = ddhazard_control(method = "GMA"))
bt <- ddhazard_boot(fit, R = 999)
plot(fit, ddhazard_boot = bt, level = .9)
# }
Run the code above in your browser using DataLab