DOING THINGS RIGHT

彩りのある生産的な日々を。

ガンマ分布のscaleとshape(rate)を直感的にわかりやすく理解する【R】

統計モデリングを学ぶ上での名著、

データ解析のための統計モデリング入門――一般化線形モデル・階層ベイズモデル・MCMC (確率と情報の科学)

(緑本)の6章8節にこのような記述がある。

この例題は、架空植物50個の葉の重量と花の重量の関係を調べます。個体iの葉重量をx _
i、花重量をy _ iとするのですが、x _
iが大きくなるにつれてy _ iも大きくなっているようです。(中略)
ここでは、ある個体の花の重量がガンマ分布に従っていることにします。応答変数y _
iは連続地ですが、重量なので正の値しかとりません。したがって、そのばらつきは正規分布ではなくガンマ分布で説明したほうがよさそうです。

なるほど、応答変数が連続値で正の値しかとらない場合はガンマ分布のほうがふさわしいのか。そのガンマ分布とやらを一度Rで見てみようということで、正規分布のグラフを書くノリでcurve()+dgammaを渡してみたところ、

続きを読む

はてなブログにRのコードを出力結果とともに載せる方法【R Markdown】

はてなブログをR Markdownで書きたい!

先日、当ブログでRに関する記事を書いたのですが、その際ソースコードを載せる方法はわかっていたものの出力結果をどう表示したらいいかわからずに苦戦しましたので、備忘録もかねて解決策を書いておきます。

結論から言うと、「R MarkdownからMarkdownに変換、それをはてなブログ(Markdown記法)に貼り付け」です。なんで思いつかなかったんだろう。。。
なおここではR Markdownについて基本的な事項は既知とします。もし初めて聞いたという方には、以下の記事が参考になります。

kazutan.github.io

続きを読む