I am trying to calculate several binomial proportion confidence intervals. My data are in a data frame, and though I can successfully extract the estimate
from the object returned by prop.test
, the conf.int
variable seems to be null when run on the data frame.
library(dplyr)
cases <- c(50000, 1000, 10, 2343242)
population <- c(100000000, 500000000, 100000, 200000000)
df <- as.data.frame(cbind(cases, population))
df %>% mutate(rate = prop.test(cases, population, conf.level=0.95)$estimate)
This appropriately returns
cases population rate
1 50000 1e+08 0.00050000
2 1000 5e+08 0.00000200
3 10 1e+05 0.00010000
4 2343242 2e+08 0.01171621
However, when I run
df %>% mutate(confint.lower= prop.test(cases, pop, conf.level=0.95)$conf.int[1])
I sadly get
Error in mutate_impl(.data, dots) :
Column `confint.lower` is of unsupported type NULL
Any thoughts? I know alternative ways to calculate the binomial proportion confidence interval, but I would really like to learn how to use dplyr
well.
Thank you!
sumcases
? :)?prop.test
in theValue
section the description forconf.int
tells us "a confidence interval for the true proportion if there is one group, or for the difference in proportions if there are 2 groups and p is not given, or NULL otherwise" so you need to test either one or two groups to generate a nonNULL
value for youconf.int
, not the 4 groups that are currently being testeddplyr
would be doing a sort of row-wise call ofprop.test
and thus the rows would each be considered individually, and the first part of your quoted section ("a confidence interval for the true proportion if there is one group") would apply. Am I misunderstandingdplyr
?library(purrr) ; library(dplyr) ; df %>% mutate(confint.lower = map2(.x = cases, .y = population, .f = ~ prop.test(.x, .y, conf.level=0.95)$conf.int[1]))
. I haven't really dug into the htest class to figure out why your version isn't working, but this should. EDIT: Actually, on quick reflection, it's probably not working becauseprop.test
isn't vectorized.