So-net無料ブログ作成
検索選択

折れ線へのあてはめ [統計]

下のグラフのような折れ線状のデータへのあてはめをおこなうモデルをかんがえる。変化点の位置も不明とする。
Rplot.png

続きを読む


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

JAGSのadaptのテスト [統計]

JAGSでは、adaptationしてからburn-inするとマニュアルにあるので、そのあたりの動作を確認してみました。

続きを読む


タグ:jags MCMC
nice!(1)  コメント(0)  トラックバック(0) 
共通テーマ:日記・雑感

Stanで、負の二項分布へのあてはめと予測区間の描画 [統計]

Rのglm.nb()では予測区間を求めるのがけっこう面倒そうなので、Stanでやってみたメモです。

続きを読む


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

過分散データに負の二項分布のモデルをあてはめ [統計]

過分散データに負の二項分布のモデルをあてはめて、p値をみてみました。MASSパッケージのglm.nb関数を使用。

続きを読む


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

過分散データへのGLM適用時のAICによるモデル選択 [統計]

こんどは、過分散データにGLMを適用して、AICでモデル選択するとどうなるかをためしました。前提として、「AICは正しいモデルを選ぶためのものではない」ことに注意。

続きを読む


タグ:AIC GLM
nice!(3)  コメント(0)  トラックバック(0) 
共通テーマ:日記・雑感

過分散データへのGLM適用時のp値 [統計]

つづいて、過分散データにGLM(誤差構造はポアソン分布、リンク関数はlog)を適用したときの検定はどうなるか、をやってみます。

続きを読む


タグ:GLM
nice!(1)  コメント(0)  トラックバック(0) 
共通テーマ:日記・雑感

過分散データへのGLM適用時のパラメータ推定値 [統計]

過分散データをむりやりGLMにあてはめたとき、パラメータの推定値がどうなるか実験してみました。

続きを読む


タグ:glmm GLM
nice!(1)  コメント(0)  トラックバック(0) 
共通テーマ:日記・雑感

[Stan] 分散を階層化したモデル [統計]

分散を階層化したモデルのメモ。

続きを読む


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

Stan example models/BPA 2.13化 [統計]

先日の Bayesian Population Analysis using WinBUGS のStan移植版の更新ですが、Stan 2.13で導入されたcompound declarationもとりいれて完了しました。

stan-dev/example-models/BPA

12.3節のN-mixutureモデルをを二変量ポアソンモデルで実装しなおしたのが、もっとも効率があがったところです。



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

[Stan] 2回測定のN-mixture model, bivariate Poisson distribution (2) [統計]

きのうのモデルですが、Stan組込みの確率質量関数をつかった方が、理解しやすくて、さらに はやいようなので、そのようにしました。

続きを読む


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