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

統計学備忘録 since2016

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

相関行列と共分散行列

相関行列

x1 <- c(93, 89, 115, 90, 75)
x2 <- c(121, 136, 121, 161, 125)
x3 <- c(101, 115, 118, 122, 106)

#データフレームにすること
dat <- data.frame(x1, x2, x3)       
 
cor(dat)

           x1         x2        x3
x1  1.0000000 -0.2220449 0.3981155
x2 -0.2220449  1.0000000 0.6767614
x3  0.3981155  0.6767614 1.0000000

散布図

pairs(dat)

f:id:yoshida931:20190604112501p:plain:w400

共分散行列

これもデータフレームで    
var(dat)

       x1     x2    x3
x1 207.80 -54.15 49.80
x2 -54.15 286.20 99.35
x3  49.80  99.35 75.30

参考

install.packages("psych")
library(psych)
x <- rnorm(100,5,3)
y <- rnorm(100,15,8)
z <- rnorm(100,2,0.5)
dat2 <- data.frame(x,y,z)
pairs.panels(dat2)

f:id:yoshida931:20190604113119p:plain:w400