1

Rのいくつかの昆虫計数データに対してゼロ膨張負の二項GLMを実行しています。私の問題は、ゼロ膨張を維持するために、Rに種データを1つの積み重ねられた列として読み取らせる方法です。小計して、Abundanceというタイトルの単一行としてRにインポートすると、ゼロが失われ、モデルが機能しなくなります。すでに、私は次のことを試みました：

  1. データを自分でスタックする（80列* 47行あります）ので、手動でスタックした後の3760行では、pscl zeroinfl()コマンドを使用したときにRがどれほど遅くなるか想像できます（私のコンピューターでは20分かかります！、それでも機能します）

  2. 次の問題は空間相関に関するものです。独立性を侵害するために同じ媒体からサンプリングされた特定のサンプラー。モデルの要素としてミディアムを入れることはできますか？

4

1 に答える 1

0

3760 行は PSCL で 20 分かかりますか? なんてこった、私は30.000行の戦いをしています:)それが私のpscl計算が終了しなかった理由です...

ただし、次に、入れ子になったランダム効果 (lme/gamm) を含む GLMM と、シータを低い値に設定する負の二項分布を使用して、分布がゼロ インフレーションを処理できるようにしました。これはゼロの次数によると思います。私の場合は 44% で、残差はかなり良さそうでした。

于 2012-09-01T19:00:24.220 に答える