The researcher used the boottwo() function from lecture. Output and the code she ran is given below:
mean(new_drug)
[1] 136
sd(new_drug)
[1] 2.4
mean(placebo)
[1] 166
sd(placebo)
[1] 2.2
BP.boot.t=boottwo(new_drug, placebo, nboot=5000)
round(quantile(BP.boot.t, probs=c(0.005, 0.025, 0.05, 0.10, 0.90, 0.95, 0.975, 0.995)),4)
0.5% 2.5% 5% 10% 90% 95% 97.5% 99.5%
-2.825 -2.017 -1.664 -1.300 1.260 1.672 1.948 2.787