Mô hình ngẫu nhiên

library(tidyverse)
library(PerformanceAnalytics)
d <- readxl::read_excel("C:/Users/tranh/Downloads/thongke/Mohinhngaunhien.xlsx")
my_data <- d[, c(2,3,4,5,6)]
my_data <- data.frame(sapply(my_data, as.numeric))
chart.Correlation(my_data, histogram=TRUE, pch=19)

corr <- cor(my_data)
round(corr, 2)

##       FDI    TO    IR    ER   INF
## FDI  1.00  0.13 -0.23  0.34 -0.07
## TO   0.13  1.00 -0.06  0.23 -0.10
## IR  -0.23 -0.06  1.00 -0.60  0.48
## ER   0.34  0.23 -0.60  1.00 -0.30
## INF -0.07 -0.10  0.48 -0.30  1.00

library(corrplot)
my_data<- cor(my_data)
corrplot(my_data, type = "upper", order = "hclust", tl.col = "black", tl.srt = 45)

heatmap(x = my_data, col = colorRampPalette(c("blue", "white", "red"))(20), symm = TRUE)