一日一膳(当社比)

RとJavaと時々数学

統計学

Rで始める傾向スコア分析による効果検証

最近因果推論の勉強を始めた管理人です. 因果効果の推定の基礎となるポテンシャルアウトカムフレームワーク, 選択バイアス, 傾向スコア(特に逆確率重み付き推定)についてRでの実装つきでまとめたので記事にしました. 大まかな内容については、以下の本を参考…

【R】Search Consoleレポートからクリック率を高めるキーワードを推定する

1. 検索キーワードとクリック率 Search Consoleで取得できる検索パフォーマンスレポートでは, Search Consoleを導入済みウェブサイトがどのような検索キーワードからアクセスされているかを確認できます. このデータを元に「どんなキーワードで自分のウェブ…

【R】エーレンフェストモデル

1.エーレンフェストモデル エーレンフェストモデル(Ehrenfest model)は, 外界から閉ざされた容器内での気体の拡散を説明する場面に使われる確率モデルの一つで、 マルコフ連鎖の例としてしばしば取り上げられます. 以下、Rでのシミュレーションと合わせてこ…

【R】WAICを用いた予測モデル評価を混合Poisson分布で試す

こちらの本の6, 7章の汎化損失・WAICの内容をまとめてみました。 社会科学のための ベイズ統計モデリング (統計ライブラリー ) 作者:浜田 宏,石田 淳,清水 裕士 出版社/メーカー: 朝倉書店 発売日: 2019/12/01 メディア: 単行本(ソフトカバー) 理論と数値…

隠れマルコフ連鎖モデルの変分推論(Rで実装)

今読んでいるベイズ推論のテキストベイズ推論による機械学習入門 の内容より、 隠れマルコフモデルの変分推論についてRで実装したので記事化してみます。 1.隠れマルコフモデルにおける変分推論 モデルの確認 個の状態の間を移りあうマルコフ連鎖を考えます…

Poisson過程とGamma分布

Poisson過程 をパラメータとし,が次の条件で与えられるPoisoon過程であるとします. i. ii.独立増分性( independent increment) 任意のに対して は独立である. iii. 任意のに対して が次の事前分布に従うとします, ただしここではガンマ分布でその確率密度…

基本情報技術者試験まであと一週間を切りました

(この行を消して、ここに「迷い」と「決断」について書いてください) #「迷い」と「決断」 Sponsored by イーアイデム 基本情報技術者試験とは 基本情報技術者試験とは,情報処理に関する知識を有することを示す国家資格です.毎年2回(4月と10月)開催され…

幾何ブラウン運動が対数正規分布に従うことの導出

目次 幾何ブラウン運動 導出 幾何ブラウン運動 幾何ブラウン運動は,確率過程の一つで株価の予測などに応用があります.幾何ブラウン運動は,各時刻において対数正規分布(対数とったら正規分布)に従いますが,今回はその導出を行いたいと思います. 幾何ブラ…

変分近似アルゴリズム

// 変分近似アルゴリズムは,ベイズモデルの事後分布を近似的に求める手段の一つです.今回の記事では,変分近似アルゴリズムの導出についてまとめ,また実際に線形回帰モデルに適用してみたいと思います. 目次 1.変分近似アルゴリズム 2.線型回帰モデルの…

生存時間分析のメモ(~Python ライブラリ lifelines~)

$$\newcommand{\Expo}[1]{ \mathrm{exp}\Bigl(#1 \Bigr)}$$ $$\newcommand{\Prob}[1]{\mathbb{P} \lbrack #1 \rbrack}$$ 生存時間分析の基礎事項についてまとめてみた。pythonの生存時間分析ライブラリであるLifelinesを使った分析例も載せています. 目次 1.…

ブラウン運動の最尤推定

keywords :時系列解析 , 最尤法 , Kullback-Leibuler情報量 , ブラウン運動 今回のテーマは、前回に引き続き時系列解析の話題です。 前回の記事ではワイトノイズ、ブラウン運動、AR(1)過程について その平均、自己共分散がどのように変化するかについて書き…

時系列解析入門

key words:時系列解析 , ブラウン運動( Brownian Motion) ,AR過程 今回のテーマは時系列解析です。 時系列解析の目的は、確率的に変化する情報(確率変数)を時刻ごとに 追って分析することです。数学、物理学、経済学などの各分野で実用の 機会も多い題材…