tests/testthat/test-survival-survreg.R

In tidyverse/broom: Convert Statistical Objects into Tidy Tibbles

skip_on_cran()

skip_if_not_installed("modeltests")
library(modeltests)

skip_if_not_installed("survival")
suppressPackageStartupMessages(library(survival))

sr <- survreg(Surv(futime, fustat) ~ ecog.ps + rx, ovarian,
  dist = "exponential"
)

test_that("survreg tidier arguments", {
  check_arguments(tidy.survreg)
  check_arguments(glance.survreg)
  check_arguments(augment.survreg)
})

test_that("tidy.survreg", {
  td <- tidy(sr)
  td2 <- tidy(sr, conf.int = TRUE)

  check_tidy_output(td)
  check_tidy_output(td2)

  check_dims(td, 3, 5)
  check_dims(td2, 3, 7)

  expect_equal(td$term, c("(Intercept)", "ecog.ps", "rx"))
  expect_equal(td2$term, c("(Intercept)", "ecog.ps", "rx"))
})

test_that("glance.survreg", {
  gl <- glance(sr)
  check_glance_outputs(gl)
})

test_that("augment.survreg", {
  check_augment_function(
    aug = augment.survreg,
    model = sr,
    data = ovarian,
    newdata = ovarian
  )
})

test_that("tidy.survreg with robust std err", {
  sr <- survreg(Surv(futime, fustat) ~ ecog.ps + rx, ovarian,
    dist = "exponential", robust = TRUE
  )
  td <- tidy(sr)
  td2 <- tidy(sr, conf.int = TRUE)

  check_tidy_output(td)
  check_tidy_output(td2)

  check_dims(td, 3, 5)
  check_dims(td2, 3, 7)
})