adjustment_formula: Specify the cycle threshold adjustment model formula

In seabbs/epict: Flexible Cycle Threshold Modelling

Specify the cycle threshold adjustment model formula

Description

Specify the cycle threshold adjustment model formula

Usage

adjustment_formula(formula = ~1, obs, beta_default = c(0, 0.1))

Arguments

formula	A model formula defaulting to ~1.
obs	A data.frame with the following variables: id: An integer vector uniquely identifying eahc infection. test_id: An integer vector uniquely identiying each test ct_value: Numeric cycle threshold value. t: Relative (to a baseline) time of the test yielding a Ct value. t_rel_uncensored: Time of test relative to the first uncensored Ct value for that id. onset_t: Relative (to a baseline) time of onset for each infection onset_t_rel_uncensored: Time of onset relative to the first uncensored Ct value for that id. (optional). NA if unavailable/asymptomatic. censored: Logical, indicating if the Ct has been censored.
beta_default	A vector of length two containing the default mean and standard deviation to use as the prior for all covariate effect sizes.

Value

A named list including the design matrix ("design") and a data.table of priors for covariate effects ("beta").

Author(s)

Sam Abbott

See Also

Functions used to design and setup models epict_formula_as_list(), epict_individual_priors_as_list(), epict_inference_opts(), epict_model_opts(), epict_obs_as_list(), epict_onset_obs_as_list(), epict_population_priors_as_list(), epict_posterior_as_prior(), piecewise_formula(), select_piecewise_parameters()

Examples

obs <- data.frame(
 age = c(1, 2, 3), cats = c(1, 2, 1), status = c("h", "l", "h")
)
adjustment_formula(~ cats + age + status, obs)

seabbs/epict documentation built on July 3, 2023, 6:42 a.m.