Nothing
test_that("MDV == 0 are handled properly", {
mod <- exmodel(301, FALSE)
dat <- exdata(301) %>%
filter(time < 72) %>%
mutate(BW = 75, SEX = 0) #base data, 4 conc to fit
dat <- bind_rows(
dat,
dat %>%
mutate(ID = 2, # ID2 =
mdv = c(1,1,0,0,0) # 3 conc to fit instead of 4
),
dat %>%
mutate(ID = 3, # ID 3 =
mdv = c(1,1,0,0,0), # Still one conc ignored, but new info on CL:
BW = c(75, 100, 75, 75, 75) # the line must be read by mrgsim, but not for OFV computation
),
dat %>%
mutate(ID = 4, #ID 4
mdv = c(1,1,0,0,0), # the line is just dropped
BW = c(75, 100, 75, 75, 75)
) %>%
filter(BW != 100)
)
est <- mapbayest(mod, dat)
#Test number of observation fitted
expect_equal(length(est$arg.ofv.id[[1]]$idDV), 4)
expect_equal(length(est$arg.ofv.id[[2]]$idDV), 3)
expect_equal(length(est$arg.ofv.id[[3]]$idDV), 3)
expect_equal(length(est$arg.ofv.id[[4]]$idDV), 3)
#Test different eta are found
# 1 != 2
# 1 != 3
# 1 != 4
# 2 != 3
# 2 == 4
# 3 != 4
eta <- get_eta(est, output = 'list')
expect_false(isTRUE(all.equal(eta[[1]], eta[[2]])))
expect_false(isTRUE(all.equal(eta[[1]], eta[[3]])))
expect_false(isTRUE(all.equal(eta[[1]], eta[[4]])))
expect_false(isTRUE(all.equal(eta[[2]], eta[[3]])))
expect_true (isTRUE(all.equal(eta[[2]], eta[[4]])))
expect_false(isTRUE(all.equal(eta[[3]], eta[[4]])))
})
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.