tests/testthat/test_aic.R

In chjackson/flexsurv-dev: Flexible Parametric Survival and Multi-State Models

test_that("Information criteria",{
  p <- 3
  n <- nrow(ovarian)
  fitg <- flexsurvreg(formula = Surv(ovarian$futime, ovarian$fustat) ~ 1, dist="gengamma")
  expect_equal(nobs(fitg), n)
  AIC(fitg)
  expect_equal(BIC(fitg), AIC(fitg, k=log(n)))
  nevent <- sum(ovarian$fustat)
  expect_equal(BIC(fitg, cens=FALSE), AIC(fitg, k=log(nevent)))
  expect_equal(BIC(fitg), BIC.flexsurvreg(fitg))
  expect_equal(BIC(fitg,cens=FALSE), BIC.flexsurvreg(fitg,cens=FALSE))
  expect_equal(AICc(fitg), AIC(fitg, k=(2*n) / (n - p - 1)))
  expect_equal(AICc(fitg,cens=FALSE), AIC(fitg, k=(2*nevent) / (nevent - p - 1)))
  expect_equal(AICC(fitg), AICc(fitg))
  expect_equal(AICc.flexsurvreg(fitg), AICc(fitg))
  expect_equal(AICC.flexsurvreg(fitg), AICc(fitg))
})