tests/testthat/test_predict.R

In logitr: Logit Models w/Preference & WTP Space Utility Parameterizations

context("predict method")
library(logitr)

model <- logitr(
    data    = yogurt,
    outcome = 'choice',
    obsID   = 'obsID',
    pars    = c('price', 'feat', 'brand')
)

data <- subset(
  yogurt, obsID %in% c(42, 13),
  select = c('obsID', 'alt', 'price', 'feat', 'brand'))

test_that("predict() uses model data if newdata == NULL", {
  p <- predict(model)
  expect_equal(nrow(p), nrow(yogurt))
})

test_that("predict() uses newdata if provided", {
  p <- predict(model, newdata = data, obsID = "obsID")
  expect_equal(nrow(p), nrow(data))
})

test_that("predict() returns the correct user-specified prediction types", {
  x <- predict(model, newdata = data, obsID = "obsID")
  expect_true(
    (! "predicted_outcome" %in% names(x)) &
    ("predicted_prob" %in% names(x))
  )
  x <- predict(model, newdata = data, obsID = "obsID", type = "prob")
  expect_true(
    (! "predicted_outcome" %in% names(x)) &
    ("predicted_prob" %in% names(x))
  )
  x <- predict(model, newdata = data, obsID = "obsID", type = "outcome")
  expect_true(
    ("predicted_outcome" %in% names(x)) &
    (! "predicted_prob" %in% names(x))
  )
  x <- predict(
    model, newdata = data, obsID = "obsID",
    type = c("outcome", "prob")
  )
  expect_true(
    ("predicted_outcome" %in% names(x)) &
    ("predicted_prob" %in% names(x))
  )
  x <- predict(
    model, newdata = data, obsID = "obsID",
    type = c("prob", "outcome")
  )
  expect_true(
    ("predicted_outcome" %in% names(x)) &
    ("predicted_prob" %in% names(x))
  )
})