一日一膳(当社比)

RとJavaと時々数学

R

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

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

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

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

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

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

【R】listにlistを要素として追加しようとして1時間溶かした

Rでlistにlistを追加しようとしただけなのに1時間溶かしてしまったので、始末書がてら記事化します。 1. 行いたかったこと つぎのようなlist l1, l2があるとします。 l1 <- list( list(x=1, y=2), list(x=3, y=4) ) l2 <- list(x=5, y=6) これらを結合してつ…

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

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

Macでrstanの実行環境を用意する

パソコン買い換えた時にrstanの実行環境を用意した時のメモです。 想定するOS/Rのバージョン OS Catalina 10.15.1 R 3.6.1 手順1.C++コンパイラの導入 Stanの実行に必要なのでc++のコンパイラをPCにいれてない人はこの手順1を行ってください。(他にXcodeとい…

「蜜蜂と遠雷」というキーワードと一緒にツイートされる作曲家をRで集計

「蜜蜂と遠雷」に関するネタ分析を行ってみた。 「蜜蜂と遠雷」について 「蜜蜂と遠雷」は恩田陸による長編小説です. ピアノコンクールの出場者を題材とした作品で、 10/4から映画が公開されています. 大まかな内容はというと、ピアノコンクールの出場者を題…

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

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

Rで複数のcsvファイルからデータをまとめて読み込む

このところ暑い日が続いてますが,本気で20度くらい気温下がってほしいですよね. 普段の業務で表形式のデータの集計なんかしている方でしたら, csvを複数一気に読み込まないといけないなんていう状況があると思います. 今日はそんな場面で使えるRコードを…

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

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

もうggらない{dplyr}

R

目次 {dplyr}パッケージ 使用例(データフレーム) 行抽出はfilter() 集計関数はsummarise() グループ化はgroup_by() ソートはarrange() {dplyr}パッケージ Rパッケージ{dplyr}でデータ整形に使えるパッケージです.毎回ggるのが面倒なのでまとめてみます.git…

1次元Ising模型をシミュレーションする.

$$\newcommand{\Expo}[1]{ \mathrm{exp}\Bigl(#1 \Bigr)}$$ $$\newcommand{\optim}[1]{ \underset{#1}{\mathrm{arg \ min}} }$$ 目次 1.1次元Ising模型を解析的に求める 1次元Ising模型 分配関数を求める R, rstanによるシミュレーション stanコード Rコード…