library(randomForest)
library(corrplot)

rawdata <- read.csv("rawdata.csv")
rawdata$CD4 <- as.numeric(rawdata$CD4) 
rawdata1 <- as.data.frame(na.roughfix(rawdata))
rawdata1 <- rawdata1[,1:ncol(rawdata1)-1]
matrix <- cor(rawdata1,method = "spearman")
pdf("cor.heatmap.pdf", width = 8, height = 8)
corrplot(matrix, method = "color", type = "full", tl.col = "black",
         addCoef.col = "black", number.cex = 0.5)
dev.off()