truediscoveryrate: Calculate true discovery rate of a sample

In HopkinsIDD/phylosamp: Sample Size Calculations for Molecular and Phylogenetic Studies

Calculate true discovery rate of a sample

Description

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

Usage

truediscoveryrate(eta, chi, rho, M, R = NULL, assumption = "mtml")

Arguments

eta	scalar or vector giving the sensitivity of the linkage criteria
chi	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: 'stsl' for the single-transmission single-linkage assumption (prob_trans_stsl()). 'mtsl' for the multiple-transmission single-linkage assumption (prob_trans_mtsl()). 'mtml' for the multiple-transmission multiple-linkage assumption (prob_trans_mtml()).

Value

scalar or vector giving the true discovery rate

Author(s)

John Giles, Shirlee Wohl, and Justin Lessler

See Also

Other discovery_rate: falsediscoveryrate()

Examples

# The simplest case: single-transmission, single-linkage, and perfect sensitivity
truediscoveryrate(eta=1, chi=0.9, rho=0.5, M=100, assumption='stsl')

# Multiple-transmission and imperfect sensitivity
truediscoveryrate(eta=0.99, chi=0.9, rho=1, M=50, R=1, assumption='mtsl')

# Small outbreak, larger sampling proportion
truediscoveryrate(eta=0.99, chi=0.95, rho=1, M=50, R=1, assumption='mtml')

# Large outbreak, small sampling proportion
truediscoveryrate(eta=0.99, chi=0.95, rho=0.5, M=1000, R=1, assumption='mtml')

HopkinsIDD/phylosamp documentation built on May 28, 2023, 3:21 a.m.