set.seed(1544)
y <- rnorm(500, 17.4, 9)
z <- data.frame(
var1 = rnorm(500, 15, 4),
var2 = rnorm(500, 23, 7.3)
)
# Optionally create artificial missing values to trigger unshaded plot
missing_indices <- sample(seq(nrow(z)), 250)
z$var1[missing_indices] <- NA
x <- paired_equivalence_test(
z, y, "criterion", scale = "relative",
relative_region_width = 0.25
)
plot(x)
Run the code above in your browser using DataLab