【R言語】分散拡大要因のヒートマップ

分散拡大要因のヒートマップ

メモです。

library(MASS)
library(corrplot)
data("Boston")
cor_Boston <- cor(Boston)

corrplot(cor_Boston,addCoef.col = TRUE)

library(tidyr)
library(ggplot2)
library(tidyverse)
#分散拡大要因
vif1<-
Boston %>% 
  dplyr::select(-medv) %>%
  cor
vif<- 1/(1-vif1^2)


vif_long <- vif %>% 
  as_tibble(rownames=NA) %>% 
  mutate(item1 = rownames(.)) %>% 
  gather(key = item2 , value = vif ,-item1)#-item1を除いて


p <-
  ggplot(data=vif_long,aes(x = item1,y = item2, fill = vif))+
           geom_tile()+geom_text(label = round(vif_long$vif,1))+
  scale_fill_gradient(low="white",high="red")
p 

vif