一个是Boruta, 一个是carat
双重降维(Boruta包和caret包)
#correlation of texture
setwd("C:/Users/Administrator/Desktop/mission/correlation")
library(readxl)
T2WI<- read_excel("Texture features and eGFR.xlsx",
sheet = "T2WI")
SWI<- read_excel("Texture features and eGFR.xlsx",
sheet = "SWI")
BOLD<- read_excel("Texture features and eGFR.xlsx",
sheet = "BOLD")
blood<- read_excel("Texture features and eGFR.xlsx",
sheet = "blood")
##step 1
#T2WI
T2WI_selected<-c()
T2WI_feature<-T2WI[,-c(1,2)]
for (feature in names(T2WI_feature)){
result_temp <- kruskal.test(T2WI[[feature]]~T2WI$Group)
if (result_temp$p.value<0.05) {
T2WI_selected<-c(T2WI_s