data("hagelloch")
fit <- twinSIR(~ household, data = hagelloch) # a simplistic model
## extract the "residual process", i.e., the fitted cumulative intensities
residuals(fit)
## assess goodness of fit based on these residuals
checkResidualProcess(fit) # could be better
Run the code above in your browser using DataLab