# Plot the estimated intensity function
plot_intensity(fit_nhpp(DataCPSim, tau = 826))
# Segment a time series using PELT
mod <- segment(bogota_pm, method = "pelt")
# Plot the estimated intensity function for the NHPP model using the
# changepoints found by PELT
plot_intensity(fit_nhpp(bogota_pm, tau = changepoints(mod)))
Run the code above in your browser using DataLab