R 배워보기-8.4. Other interesting graphs

이 다음이 스크립트랑 함수파트인데... 
함수 정의하고 스크립트 실행하는게 링크가 없다... 
뭐 어쩌라는겨... 

아 참고로 이번에도 라이브러리 하나 깔아야됩니다 

install.packages("ellipse")

ㅇㅋ ㄱㄱ

---

Correlation matrix

사실 여기는 이거 하나밖에 없음...ㅋㅋㅋㅋㅋㅋ 
근데 내가 이걸 어디서 본 것 같은디... (가물가물)

 

> set.seed(955)
> vvar <- 1:20 + rnorm(20,sd=3)
> wvar <- 1:20 + rnorm(20,sd=5)
> xvar <- 20:1 + rnorm(20,sd=3)
> yvar <- (1:20)/2 + rnorm(20, sd=10)
> zvar <- rnorm(20, sd=6)

난수를 뭐 이렇게 많이 만드냐... 

 

> data <- data.frame(vvar, wvar, xvar, yvar, zvar)
> head(data)
       vvar       wvar     xvar       yvar      zvar
1 -4.252354  5.1219288 16.02193 -15.156368 -4.086904
2  1.702318 -1.3234340 15.83817 -24.063902  3.468423
3  4.323054 -2.1570874 19.85517   2.306770 -3.044931
4  1.780628  0.7880138 17.65079   2.564663  1.449081
5 11.537348 -1.3075994 10.93386   9.600835  2.761963
6  6.672130  2.0135190 15.24350  -3.465695  5.749642

이걸 굳이 데이터프레임까지 만들어야 하냐... 

 

> library(ellipse)

다음의 패키지를 부착합니다: ‘ellipse’

The following object is masked from ‘package:car’:

    ellipse

The following object is masked from ‘package:graphics’:

    pairs

그리고 새기들아 깔아야 되는 라이브러리는 미리 말하라고... 

 

> ctab=cor(data)
> round(ctab,2)
      vvar  wvar  xvar  yvar  zvar
vvar  1.00  0.61 -0.85  0.75 -0.21
wvar  0.61  1.00 -0.81  0.54 -0.31
xvar -0.85 -0.81  1.00 -0.63  0.24
yvar  0.75  0.54 -0.63  1.00 -0.30
zvar -0.21 -0.31  0.24 -0.30  1.00

아무튼 그려봅시다 

 

> plotcorr(ctab,mar=c(0.1,0.1,0.1,0.1))

어 때깔이... 흑백이네? 

 

> colorfun=colorRamp(c("#f7cac9","#5f4b8b","#91a8d1"),space="Lab")
> plotcorr(ctab,col=rgb(colorfun((ctab+1)/2),maxColorValue=255),mar=c(0.1,0.1,0.1,0.1))

내가 색깔을 잘못 잡았나본데...? 

아무튼 그래요... 저거 근데 보통 히트맵으로 그리지 않음?