test_that("modify_standata() creates the expected list", {
# single left-censoring limit per variable:
data <- data.frame(x = rep(NA, 11), cens_x = c(rep("none", 10), "left"))
out <- modify_standata(sdata = list(), data, 999, "x", "cens_x")
expect_equal(list(Ncens_x = 1, Jcens_x = as.array(11), U_x = 999), out)
# multiple left-censoring limits per variable:
data <- data.frame(x = rep(NA, 11), cens_x = c(rep("none", 9), "left", "left"))
out <- modify_standata(sdata = list(), data, list(c(998, 999)), "x", "cens_x")
expect_equal(list(Ncens_x = 2, Jcens_x = as.array(10:11), U_x = c(998, 999)), out)
})
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.