統計– category –

library(tidyverse) library(gtsummary) #wordpressでの出力でgtsummaryがうまく動かなかったため、モデルを並べる自作関数を作成 compare_models <- function(models, labels){ procamod <- function(amodel, alabel){ broom::tidy(amodel, conf.in...

【はじめに】 先日、Rで行列計算のみで重回帰分析を実装してみましたが、この記事ではその行列計算の導出を行ってみたいと思います。個人的な勉強ノートを公開しているだけなので、もし「ここが間違っている」とかがあれば教えていただけたら非常に嬉しい...

Causal Mediation Analysisの勉強のために自分が理解するためにまとめたノートを公開します。 尚、内容はKRSKさんの素晴らしいブログの記事を思いっきり参考にしています。 因果効果のメカニズムを検討する：媒介分析(Causal Mediation Analysis)入門①～既...

Causal Mediation Analysisの勉強をはじめたところ、$Y_{a=1}$などの記号の意味がだいぶ怪しくなっていたため、改めて自分のための勉強ノートを作成しがてら、スライドにまとめてみました。 参考文献は、 What If データから因果関係をどう導く？：統計的...

library(tidyverse) library(gtsummary) #wordpressでの出力でgtsummaryがうまく動かなかったため、モデルを並べる自作関数を作成 compare_models <- function(models, labels){ procamod <- function(amodel, alabel){ broom::tidy(amodel, conf.in...

library(tidyverse) library(gtsummary) #wordpressでの出力でgtsummaryがうまく動かなかったため、モデルを並べる自作関数を作成 compare_models <- function(models, labels){ procamod <- function(amodel, alabel){ broom::tidy(amodel, conf.in...

【中心極限定理】 中心極限定理は何らかの分布（例えば上の図では「サイコロ」の出目）がある場合に、その分布から出てくる数字の平均値が、正規分布するという話です。この記事では、このことが成り立つということを、私が学んで理解できた範囲で解説を行...

ロジスティック回帰分析を頑張って実装する記事の5つめ、最後の記事です。ここでは、いよいよRを使ってロジスティック回帰分析を実装していきます。 (ブログの引っ越しにともない、summary(model)の表記がうまくアップロードできないため、broom::tidy(mod...

ロジスティック回帰分析を頑張って実装する記事の4つめです。 一つ目の記事では、データの対数尤度、$L$を求めました。 $$L = \log(\ell)$$ $$=\log \left( \prod_{i=1}^{n} p_i^{y_i}(1-p_i)^{1-y_i} \right)$$ $$= \sum_{i=1}^{n} \bigg( y_{i}\log{(p_i...

ロジスティック回帰分析を頑張って実装する記事の3つめです。 一つ目の記事では、データの対数尤度、$L$を求めました。 $$L = \log(\ell)$$ $$=\log \left( \prod_{i=1}^{n} p_i^{y_i}(1-p_i)^{1-y_i} \right)$$ $$= \sum_{i=1}^{n} \bigg( y_{i}\log{(p_i...