Two examples of multiple-bias sensitivity analysis are given in Smith, Mathur, and VanderWeele 2020. The research questions concern whether HIV infection in utero causes wasting (low weight-for-length), and whether vitamin consumption during pregnancy protects against childhood leukemia.
A study by Omoni et al. found that compared to children who were unexposed to HIV, those who had been infected with HIV in utero were significantly more likely to be below a weight-for-length Z-score of -2 as toddlers.1 The odds ratio comparing the two groups was 6.75 (95% CI 2.79, 16.31) at 2 years. The study sample was of participants in a vitamin supplementation trial, but this was an observational study with respect to HIV infection. We may be interested in the sensitivity of this unadjusted result to confounding and to selection bias, since the choice of whether to participate could have been affected by HIV status as well as other factors.
We can declare these biases. The
"increased risk" means that we are willing to make the assumption that the outcome is more likely in the selected portion of both exposure groups. The
"general" argument is in contrast to
"selected", the latter meaning that we are only interested in inference in the selected population. Since
"general" is the default, we could leave it out.
HIV_biases <- multi_bias(confounding(), selection("general", "increased risk"))
Printing the biases prints out the arguments that are required for the
multi_bound() function, making for easy copying and pasting into that function.
HIV_biases #> The following arguments can be copied and pasted into the multi_bound() #> function: RRAUc = , RRUcY = , RRUsYA1 = , RRSUsA1 =
multi_bound(biases = HIV_biases, RRAUc = 2.3, RRUcY = 2.5, RRUsYA1 = 3, RRSUsA1 = 2) #>  2.269737
To calculate a multi-bias E-value, we must provide the observed effect estimate along with the set of biases. The authors reported an odds ratio but we will assume the outcome of extreme wasting is rare enough for it to approximate a risk ratio.
multi_evalue(biases = HIV_biases, est = OR(6.75, rare = TRUE), lo = 2.79, hi = 16.31) #> This multi-bias e-value refers simultaneously to parameters RR_AUc, #> RR_UcY, RR_UsY|A=1, RR_SUs|A=1 . (See documentation for details.) #> point lower upper #> RR 6.750000 2.790000 16.31 #> Multi-bias E-values 4.635703 2.728474 NA
A study by Ross et al. examined the risk of leukemia in children whos mothers did and did not report consuming vitamin supplements during pregnancy.2 They reported an odds ratio of 0.51 (95% CI 0.30, 0.89) for acute lymphoblastic leukemia, conditional on maternal age, race, and a binary indicator of education. We may worry about additional, unmeasured confounding, as well as possible recall bias, as mothers of children with a cancer diagnosis might be more likely to report not taking vitamins even if they did do so.
When assessing exposure misclassification, we have to specify whether the outcome and exposure are rare (irrespective of the actual scale of the effect estimate).
leuk_biases <- multi_bias(confounding(), misclassification("exposure", rare_outcome = TRUE, rare_exposure = FALSE))
leuk_biases #> The following arguments can be copied and pasted into the multi_bound() #> function: RRAUc = , RRUcY = , ORYAa =
Again we can calculate the bound and multi-bias E-value as in the text.
multi_bound(biases = leuk_biases, RRAUc = 2, RRUcY = 1.22, ORYAa = 1.59) #>  1.747568
multi_evalue(biases = leuk_biases, est = OR(0.51, rare = TRUE), lo = 0.3, hi = 0.89) #> This multi-bias e-value refers simultaneously to parameters RR_AUc, #> RR_UcY, RR_YA*|a . (See documentation for details.) #> point lower upper #> RR 0.510000 0.3 0.890000 #> Multi-bias E-values 1.351985 NA 1.058404
Omoni AO, Ntozini R, Evans C, et al. Child growth according to maternal and child HIV status in Zimbabwe. Pediatr Infect Dis J. 2017;36:869–876.↩︎
Ross JA, Blair CK, Olshan AF, et al. Periconceptional vitamin use and leukemia risk in children with down syndrome: A children’s oncology group study. Ann Ny Acad Sci. 2005;104:405– 410.↩︎