PAVA order restriction under decreasing failure rate (DFR)
This function imposes the PAVA DFR order restriction by eliminating change-points violating the restriction
pava_dfr(time_die,ttot,deaths)
time_die
: event timesttot
: the total time on test (ttot) corresponding to the event timesdeaths
: the number of deaths at each event timetime2: the event times after PAVA ttot2: the corresponding ttot deaths2 the corresponding number of deaths
data(pava_dfrd) t_d = pava_dfrd[,1] t = pava_dfrd[,2] d = pava_dfrd[,3] pava_dfr(t_d, t, d)
Useful links