context("test_Eloglik.R")
test_that("Eloglik agrees with version 0.3", with(simulate(sparse = TRUE), {
original.res = readRDS('Eloglik_original_res.rds')
scaledX = set_X_attributes(X)
scaledX.sparse = set_X_attributes(X.sparse)
s$Xr = colSums(compute_MXt(s$alpha*s$mu, scaledX))
dense.res = Eloglik(scaledX, y, s)
sparse.res = Eloglik(scaledX.sparse, y, s)
expect_equal(dense.res, original.res)
expect_equal(sparse.res, original.res)
}))
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.