一般化線形モデルの勉強の一環として、条件付き期待値と確率分布の性質についてRの乱数シミュレーションで確認してみる 条件付き期待値の性質 相関関係のある乱数y1,y2を生成 set.seed(1) rho <- 0.7 x <- rnorm(n=1000,mean=0,sd=3) e1 <- rnorm(n=1000,mean=0,sd=3) e2 <- rnorm(n=1000,mean=0,sd=3) y1 <- sqrt(rho)*x+sqrt(1-rho)*e1 y2 <- sqrt(rho)*x+sqrt(1-rho)*e2 cor(y1,y2) plot(y1, y2) df<-as…