「統計学関連なんでもあり」の過去ログ---026

1019.　教えて下さい　教えて下さい　2003/11/14 (金) 00:48
└1021.　Re: 教えて下さい　青木繁伸　2003/11/14 (金) 11:38

1019.　教えて下さい　教えて下さい　　2003/11/14 (金) 00:48

M種類のボールがそれぞれ1/Mの確率で箱にたくさんある。箱の中を見ずに任意にN個のボールを取り出した時，全ての種類のボールが取り出せる確率は何？
N >= M > 0とする。

　　　　　[このページのトップへ]

1021.　Re: 教えて下さい　青木繁伸　　2003/11/14 (金) 11:38

解析的に求めることを明示されていないので，ひねくれて以下のような解を

> simulation <- function(m=2, n=4, k=1000) + { + x <- matrix(floor(runif(k*n, max=m)), ncol=n) + sum(apply(x, 1, function(arg) length(table(arg))) == m)/k + } > simulation(5,10) [1] 0.535 > > simulation(5,10) [1] 0.529 > simulation(5,10) [1] 0.516 > simulation(5,10) [1] 0.536 > simulation(5,10) [1] 0.502
あってるかな？

　　　　　[このページのトップへ]