From de7f8b6b1a234ded1dac55fd97267f5ccbc1c2be Mon Sep 17 00:00:00 2001 From: Vincent Arel-Bundock Date: Fri, 15 Sep 2023 09:19:23 -0400 Subject: [PATCH] Issue #903 --- inst/tinytest/test-pkg-ggeffects.R | 28 ++++++++++++++++++++++++++++ 1 file changed, 28 insertions(+) create mode 100644 inst/tinytest/test-pkg-ggeffects.R diff --git a/inst/tinytest/test-pkg-ggeffects.R b/inst/tinytest/test-pkg-ggeffects.R new file mode 100644 index 000000000..847a3a070 --- /dev/null +++ b/inst/tinytest/test-pkg-ggeffects.R @@ -0,0 +1,28 @@ +source("helpers.R") +using("marginaleffects") + +# Issue #903 +set.seed(123) +dat <- data.frame( + outcome = rbinom(n = 100, size = 1, prob = 0.35), + var_binom = as.factor(rbinom(n = 100, size = 1, prob = 0.2)), + var_cont = rnorm(n = 100, mean = 10, sd = 7), + groups = sample(letters[1:4], size = 100, replace = TRUE) +) +m <- glm(outcome ~ var_binom * var_cont + groups, + data = dat, family = binomial() +) + +d <- structure(list(var_binom = structure(1:2, levels = c("0", "1" +), class = "factor"), var_cont = c(9.24717241397544, 9.24717241397544 +), groups = structure(c(1L, 1L), levels = "b", class = "factor")), class = "data.frame", row.names = c(NA, +-2L)) + +p <- predictions( + m, newdata = d, by = "var_binom", hypothesis = "pairwise", transform = "exp", type = "response" +) +expect_equal(p$estimate, 0.9867827, tolerance = 0.0001) +p <- predictions( + m, newdata = d, by = "var_binom", hypothesis = "pairwise", transform = "exp" +) +expect_equal(p$estimate, 1.61968, tolerance = 0.0001) \ No newline at end of file