# Untitled

a guest
Nov 18th, 2021
51
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
1. # Relative Risk calculation
2. # sources:
4. # @see https://www.medrxiv.org/content/10.1101/2021.07.28.21261159v1.supplementary-material
5. # @see https://alexberenson.substack.com/p/more-people-died-in-the-key-clinical
6. pfizer_vaccine_deaths_november_march = 19
7. pfizer_placebo_deaths_november_march = 13
8. pfizer_vaccine_sample_size = 21926
9. pfizer_placebo_sample_size = 21921
10.
11. # Given death rate of ~15/22,000, what sample size is needed to get statistically different
12. # ratio of outcomes
13. size_multiplier = 3
14. testTable <- matrix(c(size_multiplier*(pfizer_placebo_sample_size - pfizer_placebo_deaths_november_march),
15.                   size_multiplier*(pfizer_vaccine_sample_size - pfizer_vaccine_deaths_november_march),
16.                   size_multiplier*pfizer_placebo_deaths_november_march,
17.                   size_multiplier*pfizer_vaccine_deaths_november_march),nrow = 2, ncol = 2)
18.
19.
20. oddsratio.wald(testTable)
21.
22. \$data
23. A matrix: 3 × 3 of type dbl
24. Disease1    Disease2    Total
25. Exposed1    65724   39  65763
26. Exposed2    65721   57  65778
27. Total   131445  96  131541
28. \$measure
29. A matrix: 2 × 3 of type dbl
30. estimate    lower   upper
31. Exposed1    1.000000    NA  NA
32. Exposed2    1.461605    0.972478    2.196749
33. \$p.value
34. A matrix: 2 × 3 of type dbl
35. midp.exact  fisher.exact    chi.square
36. Exposed1    NA  NA  NA
37. Exposed2    0.06711917  0.08208172  0.06625857
38. \$correction
39. FALSE