multiple_test_diagnostic.R

two_test_parallel_AND <- function(sens, spec) {
  
  list(
    sensitivity = sens*sens,
    specificity = spec + spec - spec*spec)
}

three_test_parallel_AND_2 <- function(sens, spec) {
  
  list(
    sensitivity = 3*sens*sens - 2*sens*sens*sens,
    specificity = 3*spec*spec - 2*spec*spec*spec)
}


two_test_parallel_AND(0.9, 0.9)
three_test_parallel_AND_2(0.9, 0.9)



diseased:
  test1_pos:
    p: sens
    test2_pos:
      p: sens
    test2_neg:
      p: 1 - sens
  test1_neg:
    p: 1 - sens
      test2_pos:
      p: sens
      test2_neg:
      p: 1 - sens
not_diseased:
    test1_pos:
      p: sens
      test2_pos:
        p: sens
      test2_neg:
        p: 1 - sens
    test1_neg:
      p: 1 - sens
      test2_pos:
        p: sens
      test2_neg:
        p: 1 - sens