回答:
加えて、mgcvとそのゼロ膨張ポアソン家族(ziP()
とziplss()
)、あなたはまたに見えるかもしれませんBRMSのポール・クリスチャンBürknerによってパッケージ。これは、分布モデルに適合できます(平均よりも多くの場合、モデルのゼロインフレコンポーネントは、count関数と同じように共変量の関数としてモデル化できます)。
あなたは(平均値/数、ゼロ膨張部、等のために)線形予測子を介してのいずれかで平滑化し含めることができるs()
とt2()
、それぞれ簡単な1次元または等方性2次元スプライン、または異方性テンソル積スプラインのための条件を。これは、ゼロ膨張2項、ポアソン、負の2項、およびベータ分布に加えて、ゼロ1膨張のベータ分布をサポートしています。また、ポアソンおよび負の二項応答のハードルモデルもあります(モデルのカウント部分は、ゼロカウントが生成されないように切り捨てられた分布です)。
brmsはSTANを使用してこれらのモデルに適合するため、完全にベイジアンですが、関連する情報を抽出するには、新しいインターフェイスのセットを学習する必要があります。とは言っても、このタスクだけをサポートする機能を提供するパッケージがいくつかあり、brmsにはこれらのセカンダリパッケージを利用するヘルパー関数が記述されています。あなたはSTANがインストールされて取得する必要がありますととして、あなたは、C ++コンパイラが必要になりますBRMSは evalutationためSTANコードにRを使用して定義としてモデルをコンパイルします。
glmmTMB
ブルックスら:パッケージには、この最近のbioRxivの論文に記載されています。(2017)。モデルゼロ膨張しカウントを持つデータglmmTMB
、bioRxiv、DOI:10.1101 / 132753。
ギャビン・シンプソンはまた、比較の素敵なブログの記事があるglmmTMB
とmgcv
、この目的のために:mgcvで、カウントし、ゼロ膨張したカウントGLMMsフィッティング。
brms
実際に非常に素晴らしく柔軟性のあるものを指摘してくれてありがとう。Niki Umlaufと一緒に、bamlss
さらに柔軟な回帰機能を取得するためのカウントファミリーを作成することも計画しています...しかし、これまでのところ、データ分布をカウントするためのラウンドを行いませんでした。