tests/testthat/test_single_effect_regression.R

In susieR: Sum of Single Effects Linear Regression

context("test_single_effect_regression.R")

test_that("single_effect_regression agrees with version 0.3", with(simulate(sparse=T), {
  original.res = readRDS('singleReg_original_res.rds')
  V = 0.2

  original.res$alpha = as.vector(original.res$alpha)
  original.res$mu = as.vector(original.res$mu)
  original.res$mu2 = as.vector(original.res$mu2)
  original.res$Xr = as.vector(original.res$Xr)
  original.res$lbf = as.vector(original.res$lbf)

  scaledX = set_X_attributes(X)
  scaledX.sparse = set_X_attributes(X.sparse)

  dense.res = single_effect_regression(y,scaledX,V, optimize_V="none")
  sparse.res = single_effect_regression(y,scaledX.sparse,V, optimize_V="none")

  expect_equal_SER(sparse.res, original.res)
  expect_equal_SER(dense.res, original.res)
}))