Calculate false discovery rate of identifying transmission pairs in a sample

This function calculates the false discovery rate (proportion of linked pairs that are false positives) in a sample given the sensitivity and specificity of the linkage criteria, and sample size \(M\). Assumptions about transmission and linkage (single or multiple) can be specified.

translink_fdr(sensitivity, specificity, rho, M, R = NULL, assumption = "mtml")

Arguments

sensitivity

scalar or vector giving the sensitivity of the linkage criteria

specificity

scalar or vector giving the specificity of the linkage criteria

rho

scalar or vector giving the proportion of the final outbreak size that is sampled

M

scalar or vector giving the number of cases sampled

R

scalar or vector giving the effective reproductive number of the pathogen (default=NULL)

assumption

a character vector indicating which assumptions about transmission and linkage criteria. Default = 'mtml'. Accepted arguments are:

1. 'stsl' for the single-transmission single-linkage assumption.

2. 'mtsl' for the multiple-transmission single-linkage assumption.

3. 'mtml' for the multiple-transmission multiple-linkage assumption.

Value

scalar or vector giving the true discovery rate

See also

Other transmission linkage functions: translink_expected_links_obs_mtml(), translink_expected_links_obs_mtsl(), translink_expected_links_obs_stsl(), translink_expected_links_obs(), translink_expected_links_true_mtml(), translink_expected_links_true_mtsl(), translink_expected_links_true_stsl(), translink_expected_links_true(), translink_prob_transmit_mtml(), translink_prob_transmit_mtsl(), translink_prob_transmit_stsl(), translink_prob_transmit(), translink_samplesize(), translink_tdr()

Author

John Giles, Shirlee Wohl, and Justin Lessler

Examples

# The simplest case: single-transmission, single-linkage, and perfect sensitivity
translink_fdr(sensitivity=1, specificity=0.9, rho=0.5, M=100, assumption='stsl')
#> [1] 0.4999926

# Multiple-transmission and imperfect sensitivity
translink_fdr(sensitivity=0.99, specificity=0.9, rho=1, M=50, R=1, assumption='mtsl')
#> [1] 0.1373857

# Small outbreak, larger sampling proportion
translink_fdr(sensitivity=0.99, specificity=0.95, rho=1, M=50, R=1, assumption='mtml')
#> [1] 0.5427252

# Large outbreak, small sampling proportion
translink_fdr(sensitivity=0.99, specificity=0.95, rho=0.5, M=1000, R=1, assumption='mtml')
#> [1] 0.9805463