霞と側杖を食らう

ほしいものです。なにかいただけるとしあわせです。[https://www.amazon.jp/hz/wishlist/ls/2EIEFTV4IKSIJ?ref_=wl_share]

学習記録

説明記事

■ 自己の説明 ● 学んできたこと ・大学、大学院:経済学。統計学。 ・プログラミング言語:RとSAS。ほんのすこしpython。 ● 関心事 ・生物統計学(生存時間解析や縦断データ解析)、因果推論、疫学 ● 資格 ・統計検定1級、日本酒検定3級、さんま焼き師 ■ 記事…

美味しさの表現に関する学習記録

【学習動機】 料理やお酒に関心があり、この2,3年間、どのように表現すればいいのかを考えながら、関係する本などを読んでいた。日本酒については、日本酒検定3級を取ったときに記録を書いた。 moratoriamuo.hatenablog.com 日本酒の美味しさの表現について…

カプランマイヤー曲線からデータを復元する方法の学習記録(後編)

【学習動機】 前編からの続き. moratoriamuo.hatenablog.com カプランマイヤー曲線から, 疑似的な個々のtime-to-eventデータを復元したいとする. 復元の手順としては, Kaplan-Meier curve(カプランマイヤー曲線)からGraph Digitization Software(グラフ数値…

カプランマイヤー曲線からデータを復元する方法の学習記録(前編)

【学習動機】 論文に描かれているKaplan-Meier curve(カプランマイヤー曲線)からtime-to-eventデータを再現したいということがある. たとえば, meta-analysis(メタアナリシス)やcost-effectiveness analysis(費用効果分析), その他の二次的な分析をしたいと…

日本酒アプリのさけのわデータの学習記録

code{white-space: pre-wrap;} span.smallcaps{font-variant: small-caps;} span.underline{text-decoration: underline;} div.column{display: inline-block; vertical-align: top; width: 50%;} div.hanging-indent{margin-left: 1.5em; text-indent: -1.5…

『公衆衛生の倫理学 ─国家は健康にどこまで介入すべきか』の学習記録

【学習動機】 昔、経済学を学んでいた頃、メカニズムデザインや行動経済学で、人の行動をより良い方向へ導いていく仕組みが考えられるのは、素晴らしいものだと思っていた。所与の損失関数や利得関数を最小化したり最大化したりする方法を考えることが重要だ…

Rでメモリに乗らないデータを分割して読み込んで処理する方法の学習記録

【学習動機】 Rでメモリに乗らないような大規模データを読み込んで処理するとき, どうすべきかで悩んだ. 昔, 福島『Rによるハイパフォーマンスコンピューティング』 www.amazon.co.jp は読んだが, 今の状況と違うかもしれない. 最新のプラクティス版が出たら…

Matching Adjusted Indirect Comparison (MAIC)の学習記録

【学習動機】 Matching Adjusted Indirect Comparison (MAIC)という比較の手法があることを知った. 直接比較されていない薬剤の効果を間接的に比較したい場合, よく知られているのはネットワークメタアナリシス(Network Meta-analysis : NMA)だと思われる. N…

『因果推論入門~ミックステープ:基礎から現代的アプローチまで』の学習記録

【学習動機】 Scott Cunninghamの"Causal Inference: The Mixtape"の和訳本の『因果推論入門〜ミックステープ:基礎から現代的アプローチまで』をご恵投いただきました。ざっと目を最後まで通したので、書評と言えるほど、たいそうなことは書けないが、読後…

日本酒検定3級の学習記録

【学習動機】 日本酒を飲むことは昔から好きだった。たくさん飲んで、美味しいかどうかは感じるが、それ以上の解像度は持てずにいた。1年を通して、何種類か印象的に美味しいと感じた銘柄を覚えているかどうかといった程度だった。この1年間は4合瓶を酒屋で…

ガンマ分布の乱数からディリクレ分布の乱数を生成する方法に関する学習記録

【学習動機】 ディリクレ分布からの乱数を生成したいが, そんなものはないというケースがある. たとえば, Excel周りで実装しないといけない場合など. そういう場合にどうすれば良いのか調べてみると, ディリクレ乱数はガンマ分布からの乱数で生成できること…

Immortal time biasに対する統計解析手法のシミュレーションの学習記録その2

【学習動機】 その1のつづき. moratoriamuo.hatenablog.com Immortal time bias(不死時間バイアス)に関して, Mi et al. 2016 “Use of the landmark method to address immortal person-time bias in comparative effectiveness research: a simulation study…

Immortal time biasに対する統計解析手法のシミュレーションの学習記録その1

【学習動機】 Immortal time bias (不死時間バイアス)というバイアスがある. フォローアップ開始から治療開始までに時間がある場合, その分の時間はイベントが起きないので (もしもイベントが起きていたら, 未治療群のイベントとして扱われる), その分の時間…

パーセンタイル(四分位点など)から確率分布のパラメータを推定する方法の学習記録

h1 {font-size: 34px;} h1.title {font-size: 38px;} h2 {font-size: 30px;} h3 {font-size: 24px;} h4 {font-size: 18px;} h5 {font-size: 16px;} h6 {font-size: 12px;} code {color: inherit; background-color: rgba(0, 0, 0, 0.04);} pre:not([class]) …

生存時間解析の一般化ガンマ分布のモデルの推定の学習記録

【学習動機】 一般化ガンマ分布について、前回色々調べたが、 一般化ガンマ分布の数式の学習記録 - 霞と側杖を食らう 実際に、生存時間解析の、一般化ガンマ分布のパラメトリックモデルで推定すると、上手く推定できないケースがちょこちょこある。一般化ガ…

一般化ガンマ分布の数式の学習記録

【学習動機】 生存時間解析のパラメトリックモデルを見ていると、一般化ガンマ分布というものが現れるのだが、生存時間解析の和書や日本語資料を読んでもよく分からない。式の形も分かりにくいし、ソフトウェアを見てると2種類の確率密度関数があって、何が…

生存時間解析(生存関数とハザード関数)の学習記録

【学習動機】 つい先日, 生まれてから1万日が経過した. そんなわけなので, 記念に生存時間解析について書く(生存時間解析に関する日本語の良い感じの記事少ないから書く). 生存時間解析をやるとなると, モデルの種類が, ノンパラメトリック(カプランマイヤー…

メタアナリシスの変量効果モデルのDerSimonian-Lairdの方法の学習記録

【学習動機】 メタアナリシスの変量効果モデル(ランダム効果モデル, random effects model)で異質性つまり, 研究間の分散 を推定する方法としてよく使われるDerSimonian-Lairdの方法があるのだが, これが一目だとよく分からない. その周辺について, 少し調べ…

メタアナリシスのリスク比統合のための人工データ生成の学習記録

【学習動機】 メタアナリシスの本や記事を読んでも, 実践例は全て過去のメタアナリシスの例をやるだけで, 人工データで手法が上手くいっていることの確認をしていることがほとんどない. 自分の頭の整理も兼ねて, これから気になったことを少しずつ試したり, …

『効果検証入門』の学習記録

【学習動機】 先日、TokyoRで、じゃんけんを勝ち抜いた結果、安井&ホクソエム『効果検証入門 ~正しい比較のための因果推論/計量経済学の基礎』(以下、『効果検証入門』)を頂いた。書評と言えるほど、たいそうなことは書けないが、読後感を書いていく。 ww…

LaLonde(1986)とその周辺の学習記録

【学習動機】 Rで傾向スコアマッチングをしようとしたら、パッケージのサンプルデータでだいたいLaLonde(1986) "Evaluating the Econometric Evaluations of Training Programs with Experiment Data"(https://www.researchgate.net/publication/4900843_Eva…

折れ線回帰の学習記録

【学習動機】 折れ線回帰や区分回帰(piecewise regression, segmented regressionといったところか)をこれまで触れたことが無かったのだが, 少し触れる機会があった. どうやって推定するものかと頭で考えていたのだけれど, よく分からなかったので, 調べてみ…

{wakefield}パッケージの学習記録

【学習動機】 人工データを作るとき, いつもテキトーに思いついた乱数で作っているが, 楽にできないかと思っていたら, 偶然, Rの{wakefield}パッケージを発見した. データセットを乱数で発生させるパッケージのようだが, 日本語で説明をしているものは見当た…

リリーフ運用指標MRTRの作成の学習記録その3

【(続き)学習内容】 moratoriamuo.hatenablog.com この記事の続き. 追加分析でパリーグも. 作ったものを並べる. 左上から右に順番に西武, ソフトバンク, 楽天, ロッテ, 日ハム, バッファローズ. gridExtra::grid.arrange(p2Lion + theme(legend.position = '…

リリーフ運用指標MRTRの作成の学習記録その2

【(続き)学習内容】 moratoriamuo.hatenablog.com この記事の続き. 作ったものを並べる. 左上から右に順番に巨人, 横浜, 阪神, 広島, 中日, ヤクルト. gridExtra::grid.arrange(p2G + theme(legend.position = 'none'), p2DeNA + theme(legend.position = 'n…

リリーフ運用指標MRTRの作成の学習記録その1

【学習動機】 プロ野球のペナントレースで優勝するためには, どうすればいいのかを考えた. 『メジャー』,『八月のシンデレラナイン』, 『大正野球娘。』, 『ONE OUTS』あたりを見た経験などから色々考えた結果, シーズン中のペース配分が戦略として重要だと…

変量効果モデルとAICと目で見る尤度関数の学習記録

【学習動機】 前回, AICについて少し書いた. moratoriamuo.hatenablog.com その続きで, 線形混合効果モデル(LMM)で変量効果(ランダム効果)をモデルに入れたときって, AICが使えるのか, つまりは, モデルの尤度関数が正規分布で近似できそうなのか否かという…

『AICをぶっ壊す』とAICの復習の学習記録

【学習動機】 分析に用いているモデルの評価にAICが使えるのか不安になったので, AICについて学び直すことにした. 小西・北川『情報量規準』でAICの導出を追い直しつつ, Rで数値実験しているものを探していたら, ノンパラでよく見かける竹澤先生の書いた『AI…

Bradley-Terryモデルの基本の学習記録

【学習動機】 とあるスポーツをやっていて, 過去に試合に出るメンバーを選ばないといけない役割を背負ったことがあった. その際に, 各メンバーの強さを数値化できればいいのにと思ったことがあった. 他にも連携による強さが数値化できたらいいなと思ったこと…

{mlogit}パッケージ(混合ロジットモデル)の学習記録

【学習動機】 前回の続きで少し複雑なロジットモデルをRで推定するために, {mlogit}パッケージの使い方を覚えたい. moratoriamuo.hatenablog.com 調べていたら, Train先生が{mlogit}の使い方のエクササイズを書いているのを発見した. Kenneth Train’s exerci…