ジョンとヨーコのイマジン日記

キョウとアンナのラヴラブダイアイリー改め、ジョンとヨーコのイマジン日記です。

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

ポリア・ガンマ分布を用いた多項ロジスティック回帰のギブスサンプラー

R

ポリア・ガンマ分布 まず Polson et al. (2013, [1205.0310] Bayesian inference for logistic models using Polya-Gamma latent variables) の主な結果を紹介する. ポリア・ガンマ分布の密度関数を とすると, が成り立つ. ポリア・ガンマ分布というのは, あ…

シグモイド型成長曲線のパラメータを閉じた形で最尤推定する

R

モデル 次のような確率モデルを考える.ここで は差分を意味する記号であり, を表す. また は試行回数パラメータ , 成功確率パラメータ の二項分布を表し, 記号 (実は読み方を知らない. 「チルダ」でいいのか?)は左辺の が右辺の二項分布に従う確率変数で…

変な形の尤度関数をプロットする

R

とりあえずプロット 例1: 混合正規分布 例2: 非線形回帰 例3: ガンマ分布 R のコード 参考文献 最尤推定の基礎 とりあえずプロット 例1: 混合正規分布 モデルとして次の分布を考える. (この例は『ベイズ統計の理論と方法』に出てくる.) ここで は平均0, 分…

2変量正規分布の等高線(確率楕円)の描き方

計算ノート 2変量正規分布の密度関数 は次式で表される:ここで は2次元のベクトル.分散共分散行列の逆行列を と置くと, と書ける.密度関数の等高線を描くと(関数を同じ高さのところで輪切りにすると)楕円になる. この投稿では, 2変量正規分布の密度関数を…

カトル(いか;Cuttle)

www.youtube.com Cuttle - card game rules Cuttle - Wikipedia を参考にカトルを紹介します。 マジック・ザ・ギャザリングに似た二人用ゲームだそうですが、マジック・ザ・ギャザリングをよく知らない。マジック・ザ・ギャザリングがこのゲームから影響を受…

テキストが原稿用紙何枚分か数える

原稿用紙1枚を20字×20行としています。 ここに、 なにか、 書く… 枚数: function clickBtn1() { var tx = document.getElementById("text1").value; var tv = tx.split("\n"); var tlen = tv.map(function(word){ return Math.ceil(word.length/20) }) var …

p 値関数のアニメーション

R

前置き この文書は統計的仮説検定とかを一度は学んだことがある人向けに書いている. 小ネタです.オリジナリティのない教科書がそうであるように, 1標本の t 検定を例にとることをお許しください. プロット; 信頼区間 t 検定では, 統計モデルとして正規分布を…

フィールドレコーディング(逆張り)

写真を撮るのが苦手なのでみんなが写真を撮ってるときに音を録ってたことある。 これはおそらくブルージュ(ベルギー)のどこかの時計台の鐘: これはおそらくその近くを歩いていた馬の足音: これはおそらくマラケシュ(モロッコ)のだれかのクルアーンかな…

コレスキー分解された分散共分散行列を持つ正規分布で事後分布を近似する変分ベイズ(in Julia)

アルゴリズムの導出 機械学習スタートアップシリーズ ベイズ推論による機械学習入門作者:須山敦志講談社Amazon変分事後分布に正規分布を仮定して、変分パラメータを勾配法で推定する方法はそれなりによく知られている。例えば須山「ベイズ推論による機械学習…