前回、最大エントロピー法により正規分布を導出しました。 その際の制約条件を少し変えることで、以下のガンマ分布が出てくるそうです。 今回は、このガンマ分布の導出に挑戦します。

正規分布。一番使われる確率分布ですね。 様々な統計手法や機械学習で、前提としてデータに正規分布を仮定している場合が多いです。 みなさん知っての通り正規分布の確率密度関数は ですね。本などでいきなりでてきますが、この関数がどこからきたのか疑問に…

今回は隠れマルコフモデル(HMM)についてです。 隠れマルコフモデルはRやPythonで実装するのは簡単ですが、理解するのは少々根気が必要です(私は必要でした)。 隠れマルコフモデル自体はシンプルなモデルなので動きは簡単ですが、隠れマルコフモデルで用いら…

授業アンケート、誰もが一度はやったことがあると思います。 ですが、実際の教育現場ではアンケート結果の有効な解析ができていない様です。 今回は、依頼されて制作した授業アンケートの可視化アプリについてです。 アンケート結果は数字でみるだけではあま…

分散には、不偏分散と標本分散がありますよね。 標準偏差には、不偏標準偏差と標本標準偏差がありますよね。 しかし、Rのbaseには不偏分散と不偏標準偏差を求める関数しかないです。 私は標本分散、標本標準偏差の方が知りたいことが多いです。 毎回ぐぐって…

Rのグラフの注釈に、数式表現を入れたいと思ったことはありませんか？ 例えば、なんらかの密度関数をプロットして、数式を書きたいときとか。 もちろんそんなことは後からRの外で書けばいいんですが、 やっぱりRの中だけで完結したい！！ そんなときに役立つ…

私はRをメインで使い、サブでPythonを使います。 そんな私は、「基本はRで処理して、この処理だけPythonでやりたい！」と思うときがたまにあります。 そういったとき、RからPythonスクリプトを実行すれば解決します。 だけど、コマンドライン引数や返り値の…

Rでcsvファイルの高速な書き込みと読み込みの方法。 よく使われているcsvの書き込みと読み込みの関数にwrite.csvとread.csvがある。 しかし、どちらの関数も実行時間が遅い。 私が良く使っている関数は、data.tableパッケージのfwriteとfreadです。 こっちの…