В следующем коде я использую начальную загрузку для вычисления C.I. и значение p при нулевой гипотезе о том, что два разных удобрения, применяемые к растениям томата, не влияют на урожайность растений (и альтернативой является то, что «улучшенное» удобрение лучше). Первый случайный образец (x) поступает от растений, на которых использовалось стандартное удобрение, тогда как «улучшенный» образец использовался на растениях, откуда поступил второй образец (y).

x <- c(11.4,25.3,29.9,16.5,21.1)
y <- c(23.7,26.6,28.5,14.2,17.9,24.3)
total <- c(x,y)
library(boot)
diff <- function(x,i) mean(x[i[6:11]]) - mean(x[i[1:5]])
b <- boot(total, diff, R = 10000)

ci <- boot.ci(b)
p.value <- sum(b$t>=b$t0)/b$R

Что мне не нравится в приведенном выше коде, так это то, что повторная выборка выполняется так, как если бы была только одна выборка из 11 значений (отделяя первые 5 как принадлежащие образцу x, оставляя оставшуюся часть выборке y). Не могли бы вы показать мне, как этот код должен быть изменен, чтобы рисовать повторные выборки размера 5 с заменой из первого образца и отдельные повторные выборки размера 6 из второго образца, чтобы при повторной выборке при начальной загрузке имитировался дизайн «отдельных образцов», в результате которого исходные данные?