# NOT RUN {
library(ed50simulation)
library(boot)
pavaData <- preparePava(groupS)
bootResult <- boot(data = groupS,
statistic = bootIsotonicRegression,
R = 10,
sim = 'parametric',
ran.gen = bootIsotonicResample,
mle = list(baselinePava = pavaData,
firstDose = 2.5,
PROBABILITY.GAMMA = 0.5),
baselinePava = pavaData,
PROBABILITY.GAMMA = 0.5)
bootBC.ci(tObserved = bootResult$t0[3],
tBoot = bootResult$t[, 3],
conf = 0.95)
# }
Run the code above in your browser using DataLab