This function calculates the sample size needed to identify transmission links at a predefined false discovery rate, given a final outbreak size \(N\).

translink_samplesize(
  sensitivity,
  specificity,
  N,
  R = NULL,
  tdr,
  min_pairs = 1,
  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

N

scalar or vector giving the final outbreak size

R

scalar or vector giving the effective reproductive number of the pathogen

tdr

scalar or vector giving the desired true discovery rate (1-false discovery rate)

min_pairs

minimum number of linked pairs observed in the sample, defaults to 1 pair (2 samples); this is to ensure reasonable results are obtained

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 sample size needed to meet the given conditions

Author

John Giles, Shirlee Wohl, and Justin Lessler

Examples

translink_samplesize(sensitivity=0.99, specificity=0.995, N=100, R=1, tdr=0.75)
#> [1] 10