So-net無料ブログ作成

R: High Sierraでggplot2でヒラギノフォント [統計]

macOS 10.13 High Sierraで、ggplot2でヒラギノフォントをつかうメモです。

続きを読む


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

R: ggplot2でNoto Sansフォントをつかってみるテスト [統計]

ggplot2でNoto Sansフォントをつかってみるテストです。環境はmacOS 10.13、R 3.4.1です。

続きを読む


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

R: BayesianToolsをためす(2) [統計]

前回のつづきです。今度はZero-Infrated Poissonモデル。

続きを読む


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

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

下のグラフのような折れ線状のデータへのあてはめをおこなうモデルをかんがえる。変化点の位置も不明とする。
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) 
共通テーマ:日記・雑感