Trying to calculate a combination of p1 (event in experiment group) and p0 (event in control group) with or (odds ratio) of 1.5. nnt = Number needed to treat (100/(p1-p0))
library(tidyverse)
p1 <- seq(0,1, 0.0001)
p0 <- seq(0,1,0.0001)
or <- 1.5
df <- tibble(p1 = as.numeric(), p0 = as.numeric(), nnt = as.numeric())
for (i in p1) {
for (j in p0) {
or_formula <- round((i/(1-i))/(j/(1-j)),3)
if (or_formula == or & !is.na(or_formula)) {
df <- df %>% add_row(p1 = i, p0 = j, nnt = round(1/(i-j), digits = 0))
}
}
}