myscratchpad

# MASSパッケージ、mvrnorm関数の覚書
sg <- matrix(c(1, 0.5, 0.5,1), nr=2) # 0.5の相関係数をもつ場合の相関行列
sg
vs <- mvrnorm(n=1000, mu=c(0, 0), Sigma=sg, empirical=T)
x <- vs[,1]
y <- vs[,2]

# プロット
plot(x,y, main="r = .50", pch=20)
abline(lm(y~x))


## よく統計の教科書に載ってる相関の表
# データ生成
library(MASS)
crv <- c(0, 0.2, 0.4, 0.6, 0.8, 1.0)
vrlist <- list()
for (i in 1:6) {
  sg <- matrix(c(1, crv[i], crv[i], 1), nr=2)
  vs <- mvrnorm(n=300, mu=c(0, 0), Sigma=sg, empirical=T)
  vrlist[[i]] <- vs
}

# プロット
par(mfrow=c(2,3))
for (i in 1:6) {
 mat <- vrlist[[i]]
 x <- mat[,1]
 y <- mat[,2]
 ttl <- paste("r = ", crv[i], sep="")
 plot(x,y, main=ttl, pch=20)
}

# 回帰直線を追加するプロット
win.graph()
par(mfrow=c(2,3))
for (i in 1:6) {
 mat <- vrlist[[i]]
 x <- mat[,1]
 y <- mat[,2]
 ttl <- paste("r = ", crv[i], sep="")
 plot(x,y, main=ttl, pch=20)
 abline(lm(y~x), col="blue", lwd=2)
}