Ichijoji Sagarimatsu

統計モデリング 生態学のデータ解析・麻雀の解析など

2021-01-01から1年間の記事一覧

【麻雀】運:実力 = 50 : 1 ? 着順の統計モデリング

麻雀は実力と運が勝敗を左右するゲームですが、運と実力の比率はどのくらいなのでしょうか? 7:3くらい? 8:2くらいでしょうか? あるいは100:1? RとStanを用いて、この問いに答えてみたいと思います。 使用データ データには、オンライン対戦麻雀【天鳳…

【Stan】【農学】平均と標準偏差でガンマ分布を定義する:混合ガンマ分布による除草剤抵抗性雑草の割合推定

複数の確率分布を混合した分布(混合分布)で集団における各クラスタの割合推定をすることができます。 解析には2つのガンマ分布の混合分布を用いますが、その際、ガンマ分布のパラメータとして「平均・標準偏差」を用いるのが便利です。 今回は、農業雑草…

Rでシェープファイルを図示する(Mac)

Rでシェープファイルを図示する方法はいろいろなところで紹介されていますが、 Mac環境のRでシェープファイルが図示できない問題があります。 以前GADM(https://gadm.org/)から入手したシェープファイルを図示しようとしたのですが、うまくいきませんでした…

ggplot2:別々のデータセットから描画した図に凡例を追加する

Rで別々のデータセット、あるいはワイドフォーマットからggplot2で作図する話をします。 (どうってことなくない?)と思う方も多いかもしれませんが、凡例の付け方にかなり悩んでしまったので、記事にします。 Rではロングフォーマットが求められる という…