直接上代码
#RF code1
library(randomForest)
library(readxl)
data<- read_excel('***.xlsx',,sheet=1)
# Random sampling
samplesize = 0.7 * nrow(data)
# set.seed(80)
index = sample( seq_len ( nrow ( data ) ), size = samplesize )
# # Create training and test set
datatrain = data[ index, ]
datatest = data[ -index, ]
#rf<-randomForest(ER~Tair+Tsoil+RH+VPD,data=datatrain ,mtry=4,ntree=1000)
rf<-randomForest(ER~.,data=datatrain ,mtry=4,ntree=1000)
#predict
prediction <- predict(rf, datatest )
#calculate R2 and RMSE
#RF code2
待续...关于如何选择mtry和ntree