apcamargo

bootstraps_mann_whitney <- function(bootstraps_number, sample_1, sample_2) {
  differences <- numeric()
  for (i in sample_1) {
    for (j in sample_2) {
      differences[1 + length(differences)] <- i - j
    }
  }
  bootstraps_statistic <- numeric(bootstraps_number)
  hodges_lehmann_estimator <- median(differences)
  sample_1_delta <- sample_1 - hodges_lehmann_estimator