Now that by_row() in purrr is going to be (is?) deprecated, what is the new preferred tidyverse implementation of:
somedata = expand.grid(a=1:3,b=3,c=runif(3))
somedata %>%
rowwise() %>% do(binom.test(x=.$a,n=.$b,p=.$c) %>% tidy())
It seems as if you might nest each row into a single column, and then use map(), but I'm not sure how to do that nesting operation...plus it seems like that's a little obscure. Is there a better way?
Here is one way with map
library(tidyverse)
library(broom)
do.call(Map, c(f = binom.test, unname(somedata))) %>%
map_df(tidy)
# estimate statistic p.value parameter conf.low conf.high method alternative
#1 0.3333333 1 1.00000000 3 0.008403759 0.9057007 Exact binomial test two.sided
#2 0.6666667 2 0.25392200 3 0.094299324 0.9915962 Exact binomial test two.sided
#3 1.0000000 3 0.03571472 3 0.292401774 1.0000000 Exact binomial test two.sided
#4 0.3333333 1 0.14190440 3 0.008403759 0.9057007 Exact binomial test two.sided
#5 0.6666667 2 0.55583967 3 0.094299324 0.9915962 Exact binomial test two.sided
#6 1.0000000 3 1.00000000 3 0.292401774 1.0000000 Exact binomial test two.sided
#7 0.3333333 1 0.58810045 3 0.008403759 0.9057007 Exact binomial test two.sided
#8 0.6666667 2 1.00000000 3 0.094299324 0.9915962 Exact binomial test two.sided
#9 1.0000000 3 0.25948735 3 0.292401774 1.0000000 Exact binomial test two.sided
Or with only tidyverse
functions
somedata %>%
unname %>%
pmap(binom.test) %>%
map_df(tidy)
#estimate statistic p.value parameter conf.low conf.high method alternative
#1 0.3333333 1 1.00000000 3 0.008403759 0.9057007 Exact binomial test two.sided
#2 0.6666667 2 0.25392200 3 0.094299324 0.9915962 Exact binomial test two.sided
#3 1.0000000 3 0.03571472 3 0.292401774 1.0000000 Exact binomial test two.sided
#4 0.3333333 1 0.14190440 3 0.008403759 0.9057007 Exact binomial test two.sided
#5 0.6666667 2 0.55583967 3 0.094299324 0.9915962 Exact binomial test two.sided
#6 1.0000000 3 1.00000000 3 0.292401774 1.0000000 Exact binomial test two.sided
#7 0.3333333 1 0.58810045 3 0.008403759 0.9057007 Exact binomial test two.sided
#8 0.6666667 2 1.00000000 3 0.094299324 0.9915962 Exact binomial test two.sided
#9 1.0000000 3 0.25948735 3 0.292401774 1.0000000 Exact binomial test two.sided
If you love us? You can donate to us via Paypal or buy me a coffee so we can maintain and grow! Thank you!
Donate Us With