tests/testthat/test-surrogate.R

In sure: Surrogate Residuals for Ordinal and General Regression Models

context("Surrogate response values")


test_that("surrogate works for \"polr\" objects", {

  # Skips
  skip_on_cran()
  skip_if_not_installed("MASS")

  # Load data
  data(df1, package = "sure")

  # Fit cumulative link model
  fit <- MASS::polr(y ~ x, data = df1, method = "probit")

  # Compute residuals
  set.seed(101)  # for reproducibility
  s <- surrogate(fit)  # surrogate response values
  set.seed(101)  # for reproducibility
  r <- resids(fit)  # surrogate-based residuals
  mr <- sure:::getMeanResponse.polr(fit)  # mean response

  # Expectations
  expect_equivalent(r, s - mr)

})