Blog – Kentaro Kamada

Rで連続変数のカテゴリ化

いま1から100までの連続変数があるとする

library(tidyverse)

data <- 
  tibble(
    x = 1:100
  )

data |> print(n = 100)

# A tibble: 100 × 1
        x
    <int>
  1     1
  2     2
  3 …

手動でSuperlearner

Naimi, Ashley I. & Laura B. Balzer, 2018, “Stacked Generalization: An Introduction to Super Learning,” European Journal of Epidemiology, 33(5): 459–64, (https://doi.org/10.10…

学校外教育の参加率

使用するデータは以下の3つ

信頼区間を理解したい

信頼区間の意味がわかりやすいシミュレーション

PIAAC、PISAの標準誤差

\[ \mathrm{SE}_\theta = \sqrt{h\sum_{r=1}^{R} (\hat{\theta}_{(r)} - \hat{\theta})^2} \]

Borusyak et al. (2024)

Causal Inference

Borusyak, Kirill, Xavier Jaravel & Jann Spiess, 2024, “Revisiting Event-Study Designs: Robust and Efficient Estimation,” The Review of Economic Studies, (https://doi.org/10.1…

複利計算

\[\begin{align*} & 1年目:ar^1 \\ & 2年目:ar^1 + ar^2 \\ & 3年目:ar^1 + ar^2 + ar^3 \\ & n年目:\sum_{i=1}^{n} ar^i = \frac{ar(1 - r^n)}{1 - r} \end{align*}\]

マルチアンサーの処理

googleフォームでマルチアンサー形式の設問を作ると、以下のようにカンマで区切られた形式でデータが入力される。

バージョンを指定したRのインストール

イマドキRのインストール事情？～rig、pak、p3m～

BLASの変更によるRの高速化

R言語の計算速度はそれほど速くはないが、いくつか工夫することで高速化が可能。

tensorflowをRstudioで

pythonをRStudioで使うのはreticulateパッケージによりだいぶ楽になった。

Targeted Maximum Likelihood Estimation / Targeted Minimum Loss-based Estimation (TMLE)

Causal Inference

参考ページ：https://www.khstats.com/blog/tmle/tutorial-pt2

rlang学習

マルチレベルモデルは点推定値に影響を及ぼすか？

「マルチレベルモデルでやんなきゃ係数にバイアスが…」←ホント？

Sudharsanan & Bijlsma (2021)

Causal Inference

Sudharsanan, Nikkil & Maarten J. Bijlsma, 2021, “Educational Note: Causal Decomposition of Population Health Differences Using Monte Carlo Integration and the G-Formula,” Int…

観測度数・期待度数から対数尤度を計算する方法

尤度関数を\(L\)とすると

Causal Decomposition Analysis

Causal Inference

library(tidyverse)
library(cfdecomp)
library(gapclosing)
library(causal.decomp)


d <- 
  sMIDUS |> 
  transmute(Y = health |> as.numeric(),  # outcome
            T = edu …

複数時点での処置の因果効果の推定

Causal Inference

Hernán and Robins (2020) のPart3では、複数時点での処置の因果推論について議論されている。ここでは、IPW、g-formulaを用いた因果効果の推定を実際にやってみる。

Chap12 Distributional Assumptions, Set Identification, and Sensitivity Analysis

Causal Inference

Morgan, Stephen L. & Christopher Winship, 2014, Counterfactuals and Causal Inference: Methods and Principles for Social Research, 2nd ed., Cambridge: Cambridge University…

学生生活実態調査の文字コード

文字コードがカオスで調べるのが大変だったのでメモ

新しいパソコンを買ったらやること

特に問題なくできた

Zoteroの設定備忘録

先日ようやく文献管理をはじめた．
社会学周りだと，ここらへんのメソッドが全く確立されていない&共有されないので，試行錯誤が大変だった．
以下では自分なりに見つけた方法をメモっていく．（メモらないと絶対忘れる）

Rでロジットモデル

これは簡単で，2値のfactorを従属変数，リンク関数をglm(family = binomial('logit'))のように指定してやればOK

移動平均のかゆいところに，partial = TRUE

『前処理大全』の4章で移動平均の計算が出てきたのでメモ。

RとPython書き比べ（ユークリッドの互除法）

ユークリッドの互除法は