とりあえず plot だけでもしてみるのススメ #みどりぼん

読書会統計

10/21 に開催された「データ解析のための統計モデリング入門」（以下みどりぼん）読書会の最終回にのこのこ参加。主催の yamakatsu さん、参加者＆発表者のみなさん、会場を提供してくださったドワンゴさん、ありがとうございました。懇親会はちょっと断念…

2011-06-28

PRML 8.2章「head-to-head が観測されたら独立性が失われる」のもっとわかりやすい具体例

PRML 統計確率

PRML 8.2章「条件付き独立性」では、 head-to-head の場合は「観測されると遮断が解かれる(つまり一般に条件付き独立性を持たない)」という現象の例として「車の燃料装置」が紹介されている。この例はこの例で悪くはないと思うが、ちょっと実感しにくい。ち…

2011-06-17

ベイズの公式は地味に難しいので、確率の乗法公式を2回使おう

機械学習統計確率

window.twttr = (function(d, s, id) { var js, fjs = d.getElementsByTagName(s)[0], t = window.twttr || {}; if (d.getElementById(id)) return t; js = d.createElement(s); js.id = id; js.src = "https://platform.twitter.com/widgets.js"; fjs.paren…

2010-05-06

PRML 11章の重点サンプリングと SIR を試す

PRML 機械学習統計 sampling R

PRML 11章の予習中。 p(z) = Gamma(3, 1) について、 E[ln z] を求めよう。まずは「正解」。 PRML Appendix B を見ると、Gam(tau | a,b) に対して、E[ln tau] = ψ(a) - ln b とある(ψはディガンマ関数)。 R なら簡単に計算できる。 > digamma(3) - log(1) […

2009-07-07

Excel でベイス線形回帰

統計 Excel PRML 機械学習

こっちが本番。「パターン認識と機械学習」(PRML) 上巻 p152〜p154 の分布図を書いちゃうところまで含めて Excel でやってみた。細かい説明は後にして、実際に動いているところを。観測する訓練データが増えていくと、分布が徐々にとがっていく PRML p154 …

2009-07-07

エクセルで正規分布乱数の作り方

統計 Excel

「正規分布に従った乱数が欲しいんだけど……。まあなんかスクリプトでも書くか」「それエクセルで出来るよ」「Ｒ使えよ」 PRML 読書会に向けて３章予習中。手を動かさないと理解できない人なので、「σ=0.2 のガウスノイズを加えて人工的な訓練データを」な…