data(prostate);
# Fit the survival join point model with zero join points, i.e., fit the proportional hazard relative survival model.
fit1 = joinpoint(~Year, data=prostate, numJPoints = 0);
# Get the estimates and standard errors of the annual changes of cumulative relative survival.
haz_ac = aapc(fit1, type="HAZ_AC(CS)", interval=5);
Run the code above in your browser using DataLab