# Plot the estimated intensity functionplot_intensity(fit_nhpp(DataCPSim, tau =826))# Segment a time series using PELTmod <- segment(bogota_pm, method ="pelt")# Plot the estimated intensity function for the NHPP model using the # changepoints found by PELTplot_intensity(fit_nhpp(bogota_pm, tau = changepoints(mod)))