powered by
This function imposes the PAVA IFR order restriction by eliminating change-points violating the restriction
pava_ifr(time_die,ttot,deaths)
event times
the total time on test (ttot) corresponding to the event times
the number of deaths at each event time
time2 the event times after PAVA ttot2 the corresponding ttot after PAVA deaths2 the corresponding number of deaths after PAVA
# NOT RUN { data(pava_dfrd) t_d = pava_dfrd[,1] t = pava_dfrd[,2] d = pava_dfrd[,3] pava_ifr(t_d, t, d) # }
Run the code above in your browser using DataLab