tests/testthat/test-calculate_lm.R

In jackmwolf/pcsstools: Tools for Regression Using Pre-Computed Summary Statistics

context("calculate_lm is exact")

test_that("calculate_lm betas are exact", {
  ex_data <- pcsstools_example[c("g1", "x1", "y1")]
  means <- colMeans(ex_data)
  covs <- cov(ex_data)
  n <- nrow(ex_data)

  betas <- coef(calculate_lm(means = means, covs = covs, n = n, 
                             add_intercept = TRUE))[, 1]
  expected_betas <- coef(lm(y1 ~ g1 + x1 + 1, data = ex_data))

  expect_equal(betas, expected_betas)
})

test_that("calculate_lm standard errors are exact", {
  ex_data <- pcsstools_example[c("g1", "x1", "y1")]
  means <- colMeans(ex_data)
  covs <- cov(ex_data)
  n <- nrow(ex_data)

  ses <- coef(calculate_lm(means = means, covs = covs, n = n, 
                           add_intercept = TRUE))[, 2]
  expected_ses <- coef(summary(lm(y1 ~ g1 + x1 + 1, data = ex_data)))[, 2]

  expect_equal(ses, expected_ses)
})