set.seed(123)

n <- 500

# --- Generate moderator
z <- sort(rnorm(n))
z1 <- z[z <= 0]
z2 <- z[z > 0]

# --- x -> y when m <= 0
x1 <- rchisq(length(z1), df = 4) - 4
e1 <- rchisq(length(z1), df = 3) - 3
y1 <- 0.5 * x1 + e1

# --- y -> x when m > 0
y2 <- rchisq(length(z2), df = 4) - 4
e2 <- rchisq(length(z2), df = 3) - 3
x2 <- 0.5 * y2 + e2

# Combine data
y <- c(y1, y2)
x <- c(x1, x2)
dat <- data.frame(x, y, z)

# Fit model
m <- lm(y ~ x * z, data = dat)

# Test object creation and summary
test.dda.indep <- dda.indep(formula = m,
                            pred = "x",
                            data = dat,
                            nlfun = 2,
                            B = 200)

test_that("dda.indep works", {
  expect_silent(test.dda.indep)
})

test_that("dda.indep print works", {
  expect_output(print(test.dda.indep))
})