我正在使用bootstraping
尝试dplyr
,我遇到了一行简单的代码。
使用函数bootstrap
,我发现out可以执行
library(dplyr)
library(broom)
mtcars %>% bootstrap(10) %>%
do(tidy(sample(.$cyl, 2)))
获得一个很好的直接输出
replicate x
(int) (dbl)
1 1 6
2 1 8
3 2 6
4 2 8
...
然而,获得更多变量(列)会很好,但我无法弄清楚如何。
我想像
mtcars %>% bootstrap(10) %>%
do(tidy(sample(., 2)))
或
mtcars %>% bootstrap(10) %>%
do(tidy(sample_n(2)))
会起作用,但事实并非如此。
任何线索如何我可以对几个变量进行子集化?
想象一下,我希望获得mpg
,cyl
和disp
(输出)
replicate cyl mpg disp
(int) (dbl)
1 1 6 21 ...
2 1 4 22 ...
3 2 6 ...
4 2 8 ...
...
(我随机选择两个案例sample = 2
并重复此例程(bootstrap
)10次)。