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

No.04708　Rのプログラミング　　【K】　2007/11/16(Fri) 02:51

初めて書き込みをします，卒論に取り組む大学生です。実際に観測されたものと，比較するためにシミュレーションを組みたいのですが，プログラムの知識がなく困っています。
内容
0.1 の確率であたり，残りの0.9ではずれとなるくじがあるとします（本当は数値もっと細かいのですが簡単のため）。引いたくじを戻すので確率は一定。それを何千回順番を考慮して繰り返すとき，あたりとあたりの間にあるはずれの回数を表示させたいのです。例えば，あたりを○，はずれを×とすると， ○××○＝2。○○＝0。×○×××○×○＝3，1。2や0といった数字を278個表示させるRのプログラムを組みたいのです。宜しくお願い致します。

No.04710　Re: Rのプログラミング　　【青木繁伸】　2007/11/16(Fri) 11:36

負の二項分布じゃないの？
rnbinom(278, 1, 0.1) でよいのでは？

理論分布が負の二項分布だとすれば，観察データとシミュレーション結果を比較するのではなく，適合度の検定をやればよいのでは？以下のページに負の二項分布ではなくポアソン分布の適合度の検定があるので，それを参考にする。
http://aoki2.si.gunma-u.ac.jp/lecture/GoodnessOfFitness/poissondist.html

どうしても，回り道をやりたければ，
先頭が「はずれ」のときを考慮していないし，あっているかどうか保証しないが
sapply(strsplit(paste(ifelse(runif(10000)<0.1,"o", "x"), collapse=""), "o"), nchar)