理学療法士がまとめたノート

統計学備忘録 since2016

Rを使って統計学を勉強するブログです

相関係数のイメージ

パッケージmvtnormを使用して相関係数0.0, 0.2, 0.5, 0.7, 0.8, 0.9のグラフを作成してみます

install.packages("mvtnorm")
library(mvtnorm)

共分散行列.分散を全て1に設定しているので共分散=相関係数となります.

sigma00 <- matrix(c(1,0,0,1), ncol=2)  
sigma02 <- matrix(c(1,0.2,0.2,1), ncol=2)
sigma05 <- matrix(c(1,0.5,0.5,1), ncol=2)
sigma07 <- matrix(c(1,0.7,0.7,1), ncol=2)
sigma08 <- matrix(c(1,0.8,0.8,1), ncol=2)
sigma09 <- matrix(c(1,0.9,0.9,1), ncol=2)

乱数から各相関図に300個プロットします

n00 <- rmvnorm(n=300, mean=c(0,0), sigma=sigma00)  
n02 <- rmvnorm(n=300, mean=c(0,0), sigma=sigma02)  
n05 <- rmvnorm(n=300, mean=c(0,0), sigma=sigma05)  
n07 <- rmvnorm(n=300, mean=c(0,0), sigma=sigma07)  
n08 <- rmvnorm(n=300, mean=c(0,0), sigma=sigma08)  
n09 <- rmvnorm(n=300, mean=c(0,0), sigma=sigma09)  

図を作成

par(mfrow = c(2,3))
plot(n00,xlab = "ρ=0.0",ylab="")  
plot(n02,xlab = "ρ=0.2",ylab="") 
plot(n05,xlab = "ρ=0.5",ylab="")    
plot(n07,xlab = "ρ=0.7",ylab="")
plot(n08,xlab = "ρ=0.8",ylab="")
plot(n09,xlab = "ρ=0.9",ylab="")

f:id:yoshida931:20180514142213p:plain:w800