So-net無料ブログ作成
  • ブログをはじめる
  • ログイン

Stan: 階層距離標本抽出モデル [統計]

Kéry and Royle (2016) “Advanced Hierarchical Modeling in Ecology”の8.5.2節、階層距離標本抽出モデルのうち、距離がビンに分割されているデータの方のモデルをStanに移植してみました。

続きを読む


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

R: geom_sf()をためしてみた [統計]

おくればせながら、ggplot2 3.0のgeom_sf()で地図をかいてみました。

続きを読む


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

ポアソン分布のデータに負の二項分布をあてはめてみるテスト [統計]

ポアソン分布のデータに負の二項分布をむりやりあてはめてみたテストです。

続きを読む


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

R: KFASでオフセット項をあつかう [統計]

KFASではオフセット項もあつかえるとマニュアルにあったので、確認してみました。

続きを読む


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

Stan: 悪い例2題 [統計]

Stanで(に限らず)うまく推定できない、というより、モデルがまちがっているという例です。

続きを読む


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

R: NIMBLEをためす [統計]

ちょっとNIMBLEをためしてみました。

続きを読む


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

経県マップをRでプロットする [統計]

経県マップはFlashを使用していて、自分はもうFlashをインストールしていないので、QGISで地図を作成していたりしましたが、もっとお手軽にできるようにRスクリプトをかきました。谷村さんのNipponMapパッケージを使用しています。



Rplot01.png

経県値は184です。

オリジナルとは、色が違っていたりするのですが、そのあたりは調整してみてください。

【2018-06-23 追記】
データが間違っていたのを修正しました。なお、本家経県値のサイトもFlashからSVGに変更されています。→生涯経県値
タグ:R
nice!(1)  コメント(0) 
共通テーマ:日記・雑感

不等分散のときのKruskal-Wallis検定 [統計]

不等分散のときに、Kruskal-Wallis検定が正しい結果をかえすか、ためしてみました。検索してみると、不等分散のときはKruskal-Wallis、みたいなことが かかれているところも みつかるもので。

続きを読む


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

Dockerでjupyter-notebookでMXNet [統計]

DockerでMXNetをさわってみようとおもったが、純正のmxnet/pythonではターミナルから使うようになっているみたい。

ということで、continuumio/anaconda3をインストールして、jupyter-notebookからTerminalでmxnetをインストールした。

とりあえず、これでつかえるみたい。

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