# Relative Risk calculation
# sources:
# @see https://www.fda.gov/media/151733/download
# @see https://www.medrxiv.org/content/10.1101/2021.07.28.21261159v1.supplementary-material
# @see https://alexberenson.substack.com/p/more-people-died-in-the-key-clinical
pfizer_vaccine_deaths_november_march = 19
pfizer_placebo_deaths_november_march = 13
pfizer_vaccine_sample_size = 21926
pfizer_placebo_sample_size = 21921

# Given death rate of ~15/22,000, what sample size is needed to get statistically different
# ratio of outcomes
size_multiplier = 3
testTable <- matrix(c(size_multiplier*(pfizer_placebo_sample_size - pfizer_placebo_deaths_november_march),
                  size_multiplier*(pfizer_vaccine_sample_size - pfizer_vaccine_deaths_november_march),
                  size_multiplier*pfizer_placebo_deaths_november_march,
                  size_multiplier*pfizer_vaccine_deaths_november_march),nrow = 2, ncol = 2)


oddsratio.wald(testTable)

$data
A matrix: 3 × 3 of type dbl
Disease1    Disease2    Total
Exposed1    65724   39  65763
Exposed2    65721   57  65778
Total   131445  96  131541
$measure
A matrix: 2 × 3 of type dbl
estimate    lower   upper
Exposed1    1.000000    NA  NA
Exposed2    1.461605    0.972478    2.196749
$p.value
A matrix: 2 × 3 of type dbl
midp.exact  fisher.exact    chi.square
Exposed1    NA  NA  NA
Exposed2    0.06711917  0.08208172  0.06625857
$correction
FALSE