skip_on_cran()
skip_if_offline()
skip_if_not_installed("rstanarm")
data("pbcLong", package = "rstanarm")
m1 <- download_model("stanmvreg_1")
skip_if(is.null(m1))
test_that("clean_names", {
expect_identical(
clean_names(m1),
c("logBili", "albumin", "year", "id", "sex")
)
})
test_that("find_predictors", {
expect_identical(
find_predictors(m1),
list(
y1 = list(conditional = "year"),
y2 = list(conditional = c("sex", "year"))
)
)
expect_identical(find_predictors(m1, flatten = TRUE), c("year", "sex"))
expect_identical(
find_predictors(m1, effects = "all", component = "all"),
list(
y1 = list(conditional = "year", random = "id"),
y2 = list(
conditional = c("sex", "year"),
random = "id"
)
)
)
expect_identical(
find_predictors(
m1,
effects = "all",
component = "all",
flatten = TRUE
),
c("year", "id", "sex")
)
})
test_that("find_response", {
expect_equal(
find_response(m1, combine = TRUE),
c(y1 = "logBili", y2 = "albumin")
)
expect_equal(
find_response(m1, combine = FALSE),
c(y1 = "logBili", y2 = "albumin")
)
})
test_that("get_response", {
expect_equal(nrow(get_response(m1)), 304)
expect_equal(colnames(get_response(m1)), c("logBili", "albumin"))
})
test_that("find_statistic", {
expect_null(find_statistic(m1))
})
test_that("find_variables", {
expect_identical(
find_variables(m1),
list(
response = c(y1 = "logBili", y2 = "albumin"),
y1 = list(conditional = "year", random = "id"),
y2 = list(
conditional = c("sex", "year"),
random = "id"
)
)
)
expect_identical(
find_variables(m1, flatten = TRUE),
c("logBili", "albumin", "year", "id", "sex")
)
expect_identical(
find_variables(m1, effects = "random"),
list(
response = c(y1 = "logBili", y2 = "albumin"),
y1 = list(random = "id"),
y2 = list(random = "id")
)
)
})
test_that("find_terms", {
expect_identical(
find_terms(m1),
list(
y1 = list(
response = "logBili",
conditional = "year",
random = "id"
),
y2 = list(
response = "albumin",
conditional = c("sex", "year"),
random = c("year", "id")
)
)
)
expect_identical(
find_terms(m1, flatten = TRUE),
c("logBili", "year", "id", "albumin", "sex")
)
})
test_that("n_obs", {
expect_equal(n_obs(m1), 304)
})
test_that("find_paramaters", {
expect_equal(
find_parameters(m1, component = "all"),
structure(
list(
y1 = list(
conditional = c("(Intercept)", "year"),
random = sprintf("b[(Intercept) id:%i]", 1:40),
sigma = "sigma"
),
y2 = list(
conditional = c("(Intercept)", "sexf", "year"),
random = sprintf(
c("b[(Intercept) id:%i]", "b[year id:%i]"),
rep(1:40, each = 2)
),
sigma = "sigma"
)
),
is_mv = "1"
)
)
expect_equal(
find_parameters(m1),
structure(
list(
y1 = list(
conditional = c("(Intercept)", "year"),
random = sprintf("b[(Intercept) id:%i]", 1:40)
),
y2 = list(
conditional = c("(Intercept)", "sexf", "year"),
random = sprintf(
c("b[(Intercept) id:%i]", "b[year id:%i]"),
rep(1:40, each = 2)
)
)
),
is_mv = "1"
)
)
expect_equal(
find_parameters(m1, effects = "fixed", component = "all"),
structure(
list(
y1 = list(
conditional = c("(Intercept)", "year"),
sigma = "sigma"
),
y2 = list(
conditional = c("(Intercept)", "sexf", "year"),
sigma = "sigma"
)
),
is_mv = "1"
)
)
expect_equal(
find_parameters(m1, effects = "fixed"),
structure(
list(
y1 = list(conditional = c("(Intercept)", "year")),
y2 = list(conditional = c("(Intercept)", "sexf", "year"))
),
is_mv = "1"
)
)
expect_equal(
find_parameters(m1, effects = "random", component = "all"),
structure(
list(
y1 = list(random = sprintf("b[(Intercept) id:%i]", 1:40)),
y2 = list(random = sprintf(
c("b[(Intercept) id:%i]", "b[year id:%i]"),
rep(1:40, each = 2)
))
),
is_mv = "1"
)
)
expect_equal(
find_parameters(m1, effects = "random"),
structure(
list(
y1 = list(random = sprintf("b[(Intercept) id:%i]", 1:40)),
y2 = list(random = sprintf(
c("b[(Intercept) id:%i]", "b[year id:%i]"),
rep(1:40, each = 2)
))
),
is_mv = "1"
)
)
})
test_that("get_parameters", {
expect_equal(
colnames(get_parameters(m1)),
c(
"y1|(Intercept)",
"y1|year",
"y2|(Intercept)",
"y2|sexf",
"y2|year"
)
)
expect_equal(
colnames(get_parameters(m1, effects = "all")),
c(
"y1|(Intercept)",
"y1|year",
sprintf("b[y1|(Intercept) id:%i]", 1:40),
"y2|(Intercept)",
"y2|sexf",
"y2|year",
sprintf(
c("b[y2|(Intercept) id:%i]", "b[y2|year id:%i]"),
rep(1:40, each = 2)
)
)
)
})
test_that("linkfun", {
expect_false(is.null(link_function(m1)))
expect_length(link_function(m1), 2)
})
test_that("linkinv", {
expect_false(is.null(link_inverse(m1)))
expect_length(link_inverse(m1), 2)
})
test_that("is_multivariate", {
expect_true(is_multivariate(m1))
})
test_that("clean_parameters", {
expect_identical(
clean_parameters(m1),
structure(
list(
Parameter = c(
"(Intercept)",
"year",
"(Intercept)",
"sexf",
"year",
"b[(Intercept) id:1]",
"b[(Intercept) id:2]",
"b[(Intercept) id:3]",
"b[(Intercept) id:4]",
"b[(Intercept) id:5]",
"b[(Intercept) id:6]",
"b[(Intercept) id:7]",
"b[(Intercept) id:8]",
"b[(Intercept) id:9]",
"b[(Intercept) id:10]",
"b[(Intercept) id:11]",
"b[(Intercept) id:12]",
"b[(Intercept) id:13]",
"b[(Intercept) id:14]",
"b[(Intercept) id:15]",
"b[(Intercept) id:16]",
"b[(Intercept) id:17]",
"b[(Intercept) id:18]",
"b[(Intercept) id:19]",
"b[(Intercept) id:20]",
"b[(Intercept) id:21]",
"b[(Intercept) id:22]",
"b[(Intercept) id:23]",
"b[(Intercept) id:24]",
"b[(Intercept) id:25]",
"b[(Intercept) id:26]",
"b[(Intercept) id:27]",
"b[(Intercept) id:28]",
"b[(Intercept) id:29]",
"b[(Intercept) id:30]",
"b[(Intercept) id:31]",
"b[(Intercept) id:32]",
"b[(Intercept) id:33]",
"b[(Intercept) id:34]",
"b[(Intercept) id:35]",
"b[(Intercept) id:36]",
"b[(Intercept) id:37]",
"b[(Intercept) id:38]",
"b[(Intercept) id:39]",
"b[(Intercept) id:40]",
"b[(Intercept) id:1]",
"b[year id:1]",
"b[(Intercept) id:2]",
"b[year id:2]",
"b[(Intercept) id:3]",
"b[year id:3]",
"b[(Intercept) id:4]",
"b[year id:4]",
"b[(Intercept) id:5]",
"b[year id:5]",
"b[(Intercept) id:6]",
"b[year id:6]",
"b[(Intercept) id:7]",
"b[year id:7]",
"b[(Intercept) id:8]",
"b[year id:8]",
"b[(Intercept) id:9]",
"b[year id:9]",
"b[(Intercept) id:10]",
"b[year id:10]",
"b[(Intercept) id:11]",
"b[year id:11]",
"b[(Intercept) id:12]",
"b[year id:12]",
"b[(Intercept) id:13]",
"b[year id:13]",
"b[(Intercept) id:14]",
"b[year id:14]",
"b[(Intercept) id:15]",
"b[year id:15]",
"b[(Intercept) id:16]",
"b[year id:16]",
"b[(Intercept) id:17]",
"b[year id:17]",
"b[(Intercept) id:18]",
"b[year id:18]",
"b[(Intercept) id:19]",
"b[year id:19]",
"b[(Intercept) id:20]",
"b[year id:20]",
"b[(Intercept) id:21]",
"b[year id:21]",
"b[(Intercept) id:22]",
"b[year id:22]",
"b[(Intercept) id:23]",
"b[year id:23]",
"b[(Intercept) id:24]",
"b[year id:24]",
"b[(Intercept) id:25]",
"b[year id:25]",
"b[(Intercept) id:26]",
"b[year id:26]",
"b[(Intercept) id:27]",
"b[year id:27]",
"b[(Intercept) id:28]",
"b[year id:28]",
"b[(Intercept) id:29]",
"b[year id:29]",
"b[(Intercept) id:30]",
"b[year id:30]",
"b[(Intercept) id:31]",
"b[year id:31]",
"b[(Intercept) id:32]",
"b[year id:32]",
"b[(Intercept) id:33]",
"b[year id:33]",
"b[(Intercept) id:34]",
"b[year id:34]",
"b[(Intercept) id:35]",
"b[year id:35]",
"b[(Intercept) id:36]",
"b[year id:36]",
"b[(Intercept) id:37]",
"b[year id:37]",
"b[(Intercept) id:38]",
"b[year id:38]",
"b[(Intercept) id:39]",
"b[year id:39]",
"b[(Intercept) id:40]",
"b[year id:40]",
"sigma",
"sigma"
),
Effects = c(
"fixed",
"fixed",
"fixed",
"fixed",
"fixed",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"random",
"fixed",
"fixed"
),
Component = c(
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"conditional",
"sigma",
"sigma"
),
Group = c(
"",
"",
"",
"",
"",
"Intercept: id",
"Intercept: id",
"Intercept: id",
"Intercept: id",
"Intercept: id",
"Intercept: id",
"Intercept: id",
"Intercept: id",
"Intercept: id",
"Intercept: id",
"Intercept: id",
"Intercept: id",
"Intercept: id",
"Intercept: id",
"Intercept: id",
"Intercept: id",
"Intercept: id",
"Intercept: id",
"Intercept: id",
"Intercept: id",
"Intercept: id",
"Intercept: id",
"Intercept: id",
"Intercept: id",
"Intercept: id",
"Intercept: id",
"Intercept: id",
"Intercept: id",
"Intercept: id",
"Intercept: id",
"Intercept: id",
"Intercept: id",
"Intercept: id",
"Intercept: id",
"Intercept: id",
"Intercept: id",
"Intercept: id",
"Intercept: id",
"Intercept: id",
"Intercept: id",
"Intercept: id",
"year: id",
"Intercept: id",
"year: id",
"Intercept: id",
"year: id",
"Intercept: id",
"year: id",
"Intercept: id",
"year: id",
"Intercept: id",
"year: id",
"Intercept: id",
"year: id",
"Intercept: id",
"year: id",
"Intercept: id",
"year: id",
"Intercept: id",
"year: id",
"Intercept: id",
"year: id",
"Intercept: id",
"year: id",
"Intercept: id",
"year: id",
"Intercept: id",
"year: id",
"Intercept: id",
"year: id",
"Intercept: id",
"year: id",
"Intercept: id",
"year: id",
"Intercept: id",
"year: id",
"Intercept: id",
"year: id",
"Intercept: id",
"year: id",
"Intercept: id",
"year: id",
"Intercept: id",
"year: id",
"Intercept: id",
"year: id",
"Intercept: id",
"year: id",
"Intercept: id",
"year: id",
"Intercept: id",
"year: id",
"Intercept: id",
"year: id",
"Intercept: id",
"year: id",
"Intercept: id",
"year: id",
"Intercept: id",
"year: id",
"Intercept: id",
"year: id",
"Intercept: id",
"year: id",
"Intercept: id",
"year: id",
"Intercept: id",
"year: id",
"Intercept: id",
"year: id",
"Intercept: id",
"year: id",
"Intercept: id",
"year: id",
"Intercept: id",
"year: id",
"Intercept: id",
"year: id",
"Intercept: id",
"year: id",
"",
""
),
Response = c(
"y1",
"y1",
"y2",
"y2",
"y2",
"y1",
"y1",
"y1",
"y1",
"y1",
"y1",
"y1",
"y1",
"y1",
"y1",
"y1",
"y1",
"y1",
"y1",
"y1",
"y1",
"y1",
"y1",
"y1",
"y1",
"y1",
"y1",
"y1",
"y1",
"y1",
"y1",
"y1",
"y1",
"y1",
"y1",
"y1",
"y1",
"y1",
"y1",
"y1",
"y1",
"y1",
"y1",
"y1",
"y1",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y2",
"y1",
"y2"
),
Cleaned_Parameter = c(
"(Intercept)",
"year",
"(Intercept)",
"sexf",
"year",
"id:1",
"id:2",
"id:3",
"id:4",
"id:5",
"id:6",
"id:7",
"id:8",
"id:9",
"id:10",
"id:11",
"id:12",
"id:13",
"id:14",
"id:15",
"id:16",
"id:17",
"id:18",
"id:19",
"id:20",
"id:21",
"id:22",
"id:23",
"id:24",
"id:25",
"id:26",
"id:27",
"id:28",
"id:29",
"id:30",
"id:31",
"id:32",
"id:33",
"id:34",
"id:35",
"id:36",
"id:37",
"id:38",
"id:39",
"id:40",
"id:1",
"id:1",
"id:2",
"id:2",
"id:3",
"id:3",
"id:4",
"id:4",
"id:5",
"id:5",
"id:6",
"id:6",
"id:7",
"id:7",
"id:8",
"id:8",
"id:9",
"id:9",
"id:10",
"id:10",
"id:11",
"id:11",
"id:12",
"id:12",
"id:13",
"id:13",
"id:14",
"id:14",
"id:15",
"id:15",
"id:16",
"id:16",
"id:17",
"id:17",
"id:18",
"id:18",
"id:19",
"id:19",
"id:20",
"id:20",
"id:21",
"id:21",
"id:22",
"id:22",
"id:23",
"id:23",
"id:24",
"id:24",
"id:25",
"id:25",
"id:26",
"id:26",
"id:27",
"id:27",
"id:28",
"id:28",
"id:29",
"id:29",
"id:30",
"id:30",
"id:31",
"id:31",
"id:32",
"id:32",
"id:33",
"id:33",
"id:34",
"id:34",
"id:35",
"id:35",
"id:36",
"id:36",
"id:37",
"id:37",
"id:38",
"id:38",
"id:39",
"id:39",
"id:40",
"id:40",
"sigma",
"sigma"
)
),
class = c("clean_parameters", "data.frame"),
row.names = c(NA, -127L)
)
)
})
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.