Compute raw and pearson residuals for point process models
pp_residual(object, events, start = 0, end = max(events), steps = 1000)
object | point process model containing the parameters |
---|---|
events | vector of event times |
start | start of observation period (default 0) |
end | end of observation period (default final event) |
steps | number of steps for numeric integration (if needed) |
the raw and pearson residuals
Q <- matrix(c(-0.4, 0.4, 0.2, -0.2), ncol = 2, byrow = TRUE) x <- pp_mmhp(Q, delta = c(1 / 3, 2 / 3), lambda0 = 0.9, lambda1 = 1.1, alpha = 0.8, beta = 1.2 ) y <- pp_simulate(x, n = 10) pp_residual(x, events = y$events)#> $raw #> [1] 2.938867 #> #> $pearson #> [1] 3.097838 #>