12/02/02 15:19:45.32
知恵を貸してください。
0～6ぐらいまでの整数を要素にもつサイズが50のベクトルを作りたい。
できればポアソン分布に従って。
ここまでなら、
> rpois(50,2.5)
[1] 1 1 1 4 4 3 0 3 3 3 0 1 1 0 4 5 2 2 2 2 4 0 4 4 2 1 3 1 1 3 3 2 1 1 3 4 3 3
[39] 3 3 4 1 1 1 5 5 1 0 3 6
こんな感じでできるんだけど、このベクトルの合計をきっちり150にしたい。
固定したいのは、ベクトルサイズ50とベクトル合計値150だけで、lambdaとかは適当でいい。
当たり前だけど、合計値は変動する。
> sum(rpois(50,2.5))
[1] 132
> sum(rpois(50,2.5))
[1] 131
> sum(rpois(50,2.5))
[1] 118
こういうときはどうすればいいんだろうか。
よろしくお願いします。