傾向スコア分析 傾向スコアについてＲでやっているサイトがあったので、それをほぼもろパクリして自分でやってみた。 install.packages("Matching") library(Matching) data(lalonde) logi <- glm(treat~age+educ+black+hisp+married+nodegr+re74+re75+u74+u…

検出力について 数理統計の教科書を読んでいたら、検定における「検出力」に関する内容があって、すぐにはピンと来なかったので、Ｒで検出力を求める関数をつくってみた。 （ここへ来てより道しすぎな感は否めないが…。）正規母集団 N(mu,1)からn個の無作為…

plot3dについて 先生のアドバイスを頂いて、plot3dについてようやく分かった。あるサイトに library(rgl) z <- matrix(1:20, ncol = 4, nrow = 5) x <- c(1, 2, 3, 4) # z 行列の各行の座標 y <- c(1, 2, 3, 4, 5) # z 行列の各列の座標 plot3d(x, y, z) こ…

回帰モデルによる交絡の制御 2 ただの入力ミスだった。 k0 <- 0.2 k1 <- 0.3 k2 <- 2 N <- 10000 n <- 100 out. <- function(x){ k0+k1*x[1]+k2*x[2]+rnorm(1) } a1.d <- sample(20:80,N,replace=T) a2.d <- sample(0:1,N,replace=T) a3.d <- apply(cbind(a1…

ロジスティック回帰モデル ３Ｄプロットは僕の理解の範疇を超えているようである。意味が分からなさすぎて発狂しそうになったので、諦めてperspでグラフをつくってみた。 age <- c(18,21,22,25,26,28,33,34,35,37,42,47,55,56,58,61,65,68,75,77) out <- c(0…