SSブログ

Stan: Reduce Sum をつかった N-mixture model (2) bivariate Poisson 版 [統計]

Reduce Sum をつかった N-mixture modelですが、繰り返し測定が2回なら、bivariate Poisson分布でさらに高速化できます。

続きを読む


タグ:STAn
nice!(1)  コメント(0) 
共通テーマ:日記・雑感

Stan: Reduce Sum をつかった N-mixture model [統計]

きのうのコードは まちがっていたので、ちゃんと うごくようになおしました。

続きを読む


タグ:STAn
nice!(1)  コメント(0) 
共通テーマ:日記・雑感

Stan: Reduce Sum でなぜか結果がかわる [統計]

N-mixtureモデルをStanであてはめるのに個体数ごとの尤度を足し上げる場合、Reduce Sumをつかってマルチスレッドにすると高速化できるのではないかと、やってみました。しかし、マルチスレッドを有効にするとなぜか結果がかわってしまいました。

【2020-10-25 07:43 JST 追記】 一晩ねて落ち着いて考えてみると、log_sum_expでたすべきところを、そのままたしてしまうことになるので、結果がかわるのは当然でした。

続きを読む


タグ:STAn
nice!(1)  コメント(0) 
共通テーマ:日記・雑感

[Stan] 隠れマルコフモデルのテスト2 [統計]

遷移確率と出力確率を推定してみます。

続きを読む


タグ:STAn
nice!(1)  コメント(0) 
共通テーマ:日記・雑感

[Stan] 隠れマルコフモデルのテスト1 [統計]

CmdStan 2.24に隠れマルコフモデルが くみこまれたということでテストです。

続きを読む


タグ:STAn
nice!(1)  コメント(0) 
共通テーマ:日記・雑感

[Stan] beta_proportion 分布 [統計]

Stanにbeta_proportionという分布がはいっていたのでした。平均をつかってベータ分布のあてはめができます。

続きを読む


タグ:STAn
nice!(1)  コメント(0) 
共通テーマ:日記・雑感

[Stan] Modeling wiggly covariate relationships (AHM 10.14) [統計]

Applied Hierarchical Modeling in Ecology 10.14節の「曲がりくねった共変量関係のサイト占有モデル」をStanに移植しました。

続きを読む


タグ:AHM STAn
nice!(1)  コメント(0) 
共通テーマ:日記・雑感

[Stan] ゼロ過剰2項分布で事後予測チェック [統計]

ゼロ過剰2項分布へのあてはめ結果に事後予測チェックを適用した例と、むりやり通常の2項分布をあてはめた場合との比較です。

続きを読む


nice!(1)  コメント(0) 
共通テーマ:日記・雑感

[Stan] ランダム効果のあるモデルで新しい群のデータを加えたとき [統計]

ランダム効果のあるモデルで新しい群のデータを加えたときのパラメータの変化を見てみるということをやってみたメモです。

続きを読む


タグ:STAn
nice!(1)  コメント(0) 
共通テーマ:日記・雑感