Propagate component test specificity into panel specificity

Usage

uncertain_panel_spec(
  false_pos_controls = NULL,
  n_controls = NULL,
  ...,
  spec = spec_prior(),
  samples = 1000,
  na.rm = FALSE,
  fit_beta = FALSE
)

Arguments

false_pos_controls: the number of positives that appeared in the specificity disease-free control group. These are by definition false positives. This is (1-specificity)*n_controls
n_controls: the number of controls in the specificity disease-free control group.
...: not used
spec: the prior specificity of the test as a beta_dist.
samples: number fo random draws of sensitivity and specificity
na.rm: remove missing values
fit_beta: return the result as a beta_dist object?

Value

a vector of possible specificities for the panel or a fitted beta_dist

Examples

uncertain_panel_spec(c(2,3,4,2,2), c(800,800,800,800,800), fit_beta=TRUE)
#> 98.4% [97.4%—99.1%] (N=831.3)