2012-09-23から1日間の記事一覧

rbinomのprob引数は複数評価してくれるっぽい

R

以下の二つの書き方は等価っぽい。 x <- rbinom(900, 10, prob = seq(0.001, 0.9, 0.001)) y <- sapply(seq(0.001, 0.9, 0.001), function(x) rbinom(1, 10, prob = x)) 実際、平均と分散はほぼ一致。 > c(mean(x), sd(x)) [1] 4.524444 2.976449 > c(mean(y…