[Stan] 2回測定のN-mixture model, bivariate Poisson distribution: Taglibro de H

[Stan] 2回測定のN-mixture model, bivariate Poisson distribution [統計] ある生物の個体数を測定するとします。調査箇所数がRで、調査箇所iの個体数Niは、平均λのポアソン分布にしたがうとします。各調査地でT回の観測をおこないます。検出確率をpとすると、調査箇所iのt回目の観測数はbinominal(Ni, p)にしたがいます。すなわち Ni ~ Poisson(λ) Yi,t ~ Binomial(Ni, p) こうしたデータからλとpを推定します。 Stanでは離散パラメーターがあつかえないので、Nを周辺化してモデリングしていましたが、multivariate Poisson分布でモデリングできるというので、そちらで実装してみました。ただし、ここではT=2とします。 参考にしたところ Dennis et al. (2015) Com