test_that("Bootstrap-Prozess läuft erfolgreich, mehrere Variablen", {
  
  # Erstelle ein Beispiel für ein bsw-Objekt
  set.seed(123)
  
  # Erstellen von 5 unabhängigen Variablen mit einer Normalverteilung
  x1 <- rnorm(1000, 50, 10)
  x2 <- rnorm(1000, 30, 5)
  x3 <- rnorm(1000, 40, 8)
  x4 <- rnorm(1000, 60, 12)
  x5 <- rnorm(1000, 55, 7)
  
  # Berechnung der abhängigen Variablen y mit einer logistischen Transformation
  logit <- -4 + x1 * 0.04 + x2 * 0.05 + x3 * 0.02 + x4 * -0.03 + x5 * 0.01
  p <- 1 / (1 + exp(-logit))  # Logistische Transformation zu Wahrscheinlichkeiten
  y <- rbinom(1000, 1, p)  # Binäre Zielvariable
  
  # Erstellen eines Datensatzes
  data <- data.frame(y, x1, x2, x3, x4, x5)
  
  fit <- aRD(formula = y ~ x1 + x2 + x3, data = data)
  
  # Teste die Ausführung von bootBSW
  result <- bootaRD(fit, ci_level = 0.90, R = 1000L)
  
  # Überprüfe, dass alle Bootstrap-Replikate erfolgreich durchgeführt wurden
  expect_true(is.numeric(result$boot_obj$t0))  # Sicherstellen, dass die Originalschätzungen berechnet wurden
  expect_equal(length(result$boot_obj$t[,1]), 1000)  # Sicherstellen, dass es 1000 Replikate gibt
  expect_gt(sum(complete.cases(result$boot_obj$t)), 950)  # Überprüfen, dass mindestens 95% erfolgreich sind
})