tests/testthat/test_single_effect_regression.R

In stephenslab/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)
# }))