purrrを使ってコインフリップを生成する

Question

今、私はrでpurrrパッケージを使用する方法を学び、1、2、...、99,100コインフリップのそれぞれ5サンプルを生成する方法について考えています。 私のイメージでは、リストを作成することです..

[[1]] 

[[1]] 
[1] 1 0 1 0 0 

[[2]] 

[[1]] 
[1] 1 0 0 0 1 

[[2]] 
[1] 0 1 0 1 1 

[[3]] 

[[1]] 
[1] 0 1 1 1 0 

[[2]] 
[1] 1 0 0 0 1 

[[3]] 
[1] 0 1 1 1 1 

.. 

のようになります。その誰も私がこれを作る手助けができますか？しかし

library(purrr) 
1:100 %>% map(function(x) rerun(x, rbinom(5,1,.5))) 

次のように

Answer 1

あなたはマップ機能を使用して、ベクター1:100の各要素に適用される機能rerunをしたい、複製のためのデフォルトが生成することであるreplicateを、使用することも簡単です列方向アレイ。

lapply(1:100, function(x) replicate(x,rbinom(5,1,0.5))) 

この場合、ベースRの表現がはるかに高速であることに注意してください。あなたは多数表現の法則でちょうど興味を持っている場合

a <- function() 1:100 %>% map(function(x) rerun(x, rbinom(5,1,.5))) 
b <- function() lapply(1:100, function(x) replicate(x,rbinom(5,1,0.5))) 

library(microbenchmark) 
microbenchmark(a(),b()) 

Unit: milliseconds 
expr  min  lq  mean median  uq  max neval cld 
    a() 96.89941 104.83822 117.10245 111.48309 120.28554 391.9411 100 b 
    b() 16.88232 18.47104 23.22976 22.20549 26.31445 49.0042 100 a 

編集、コメント欄にご質問については、次のように行うことができます。

plot(1:100, do.call("c", lapply(b(), mean)), 
type= "l", xlab = "replications", 
ylab = "proportion of heads") 
abline(h = .5) 

Answer 2

私が正しくあなたを理解していれば、これはあなたが後にしているものです：

lapply(1:100, function(x) replicate(x,rbinom(5,1,0.5),simplify = FALSE))