tests/testthat/test_update_each_effect.R

In stephenslab/susieR: Sum of Single Effects Linear Regression

context("test_update_each_effect.R")

# test_that("update_each_effect agrees with version 0.3", with(simulate(sparse=T), {
#   original.res = readRDS('vbupdate_original_res.rds')
#   original.res$Xr = as.vector(original.res$Xr)
#   scaledX = set_X_attributes(X)
#   scaledX.sparse = set_X_attributes(X.sparse)
#
#   dense.res = update_each_effect(scaledX,y,s)
#   sparse.res = update_each_effect(scaledX.sparse,y,s)
#
#   expect_equal_susie_update(sparse.res, original.res)
#   expect_equal_susie_update(dense.res, original.res)
# }))