tests/testthat/test-condlSample.R

In markwh/condSample: Generate random samples from the assumed error distribution of a regression model, conditional on predictor variable values.

context("conditional samples from model output")


test_that("quantiles are as expected", {
  x = runif(1000, 0, 100)
  y = x / 4 + rnorm(1000, 0, 10)
  datxy = data.frame(x = x, y = y)
  
  lmxy = lm(y ~ x, datxy)
  q10 = condlSample(lmxy, quantile = 0.1)
  q90 = condlSample(lmxy, quantile = 0.9)
  
  expect_lt(sum(q10 > y), 150)
  expect_gt(sum(q10 > y), 50)
  
  expect_lt(sum(q90 < y), 150)
  expect_gt(sum(q90 < y), 50)
})

# test_that("condlSample works with type = 'terms'", {
#   expect_is()
# })