Nothing
test_that("censoring output is valid", {
#The below are a tests that check for valid stratification operations and blq/alq output
obs_data <- as.data.table(tidyvpc::obs_data)
sim_data <- as.data.table(tidyvpc::sim_data)
# Subest MDV = 0
obs_data <- obs_data[MDV == 0]
sim_data <- sim_data[MDV == 0]
# Add LLOQ
obs_data$LLOQ <- obs_data[, ifelse(STUDY == "Study A", 50, 25)]
#Assign observed and simulated data to tidyvpc object
vpc <- observed(obs_data, x = TIME, y = DV )
vpc <- simulated(vpc, sim_data, y = DV)
vpc <- censoring(vpc, blq = DV < LLOQ, lloq = LLOQ)
vpc <- stratify(vpc, ~STUDY)
vpc <- binning(vpc, bin = NTIME)
vpc <- vpcstats(vpc, qpred = c(0.1, 0.5, 0.9))
#Generate logical for blq in observed data
obs_strat_blq <- obs_data[, num_blq := ifelse(DV < LLOQ, TRUE, FALSE), by = STUDY]
#Check sum of TRUE logicals in obs_strat_
expect_equal(sum(vpc$obs$blq), sum(obs_strat_blq$num_blq))
vpc <- observed(obs_data, x = TIME, y = DV )
vpc <- simulated(vpc, sim_data, y = DV)
vpc <- censoring(vpc, alq = DV > LLOQ, uloq = LLOQ)
vpc <- stratify(vpc, ~STUDY)
vpc <- binning(vpc, bin = NTIME)
vpc <- vpcstats(vpc, qpred = c(0.1, 0.5, 0.9))
#Generate logical for blq in observed data
obs_strat_alq <- obs_data[, num_alq := ifelse(DV > LLOQ, TRUE, FALSE), by = STUDY]
#Check sum of TRUE logicals in obs_strat_
expect_equal(sum(vpc$obs$alq), sum(obs_strat_alq$num_alq))
})
test_that("censoring errors are expected", {
obs_data <- obs_data[MDV == 0]
sim_data <- sim_data[MDV == 0]
# Add LLOQ
obs_data$LLOQ <- obs_data[, ifelse(STUDY == "Study A", 50, 25)]
vpc <- observed(obs_data, x = TIME, y = DV )
vpc <- simulated(vpc, sim_data, y = DV)
expect_error(censoring(vpc, blq = DV < LLOQ, lloq = NA))
expect_error(censoring(vpc, blq = DV < LLOQ, lloq = NULL))
expect_error(censoring(vpc, blq = NA, lloq = NA))
expect_error(censoring(vpc, blq = 10))
expect_error(censoring(vpc, blq = TRUE, lloq = -Inf))
expect_error(censoring(vpc, alq = 100))
expect_error(censoring(vpc, alq = NA))
expect_error(censoring(vpc, alq = TRUE, uloq = NULL))
expect_error(censoring(vpc, alq = TRUE, uloq = NA))
expect_error(censoring(vpc, alq = TRUE, uloq = Inf))
})
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.