tests/testthat/test-ci.formula.R

In pROC: Display and Analyze ROC Curves

library(pROC)
data(aSAH)

context("ci.formula")

test_that("bootstrap cov works with smooth and !reuse.auc", {
	skip_slow()
	if (getRversion() > "3.6.0") {
		suppressWarnings(RNGkind(sample.kind="Rounding"))
	}
	
	for (pair in list(
			list(ci, list()),
			list(ci.se, list(boot.n = 10)),
			list(ci.sp, list(boot.n = 10)),
			list(ci.thresholds, list(boot.n = 10)),
			list(ci.coords, list(boot.n = 10, x = 0.5)),
			list(ci.auc, list()))) {
		fun <- pair[[1]]

		# First calculate ci with .default
		args.default <- c(
			list(response = aSAH$outcome,
				 predictor = aSAH$s100b),
			pair[[2]])
		set.seed(42) # For reproducible CI
		obs.default <- do.call(fun, args.default)

		# Then with .formula
		args.formula <- c(
			list(formula = outcome ~ s100b,
				 data = aSAH),
			pair[[2]])
		set.seed(42) # For reproducible CI
		obs.formula <- do.call(fun, args.formula)

		# Here we check both returned the same result
		# We ignore attributes, as we have different
		# roc objects, and unfortunately equivalent means
		# we only test near equality
		expect_equivalent(obs.default, obs.formula)
	}
})