1

こんにちは、助けてくれてありがとう！

しばらく前に、確率分布に従って値をランダムに選択することについて質問しました。これは関連していますが、独自の投稿に値すると思います。

前回の質問で作成したベクトルはバイナリでしたが、今度は加重ベクトル (境界付き整数) を生成したいと考えています。長い裾を持つゼロ膨張分布または準ポアソン分布からサンプリングしているため、ゼロを選択する確率は別の値よりもはるかに高くなりますが、大きな値 (たとえば 63) を選択する確率は有限です。

rpois を使用して、ポアソン分布から値を選択し、特定の長さのベクトルを作成できます。これは私がやりたいことと似ているので、例として使用します。

e=seq(0:63)
vec<-c(0,0,0,1,1,1)
ones <- which(vec == 1L)
temp=rpois((sum(vec)),e)
vec[ones]<-temp

これは、ポアソン分布から選択された特定の数の値をベクトルに割り当てるのに適しています。準ポアソンまたはゼロ膨張にする方法はありますか？

4

1 に答える 1

4

ここにはさまざまなディストリビューションの大きなリストがあります: http://cran.r-project.org/web/views/Distributions.html

ゼロ膨張ポアソンのために...

install.packages("gamlss.dist")
library(gamlss.dist)

rZIP(n, mu, sigma)

quasi-poisson の場合、quasipoissonff を使用した VGAM パッケージ内にいくつかの機能があるように見えますが、それは生成ではなくフィッティング用のようです。Arthur Charpentier はここで何かを考えていたようですが、ディストリビューションを正しく行うには、何を探しているのかを知る必要があります: http://freakonometrics.blog.free.fr/index.php?post/2010/10 /21/準ポアソン分布から変数を生成する方法

于 2011-08-12T01:44:17.330 に答える