albert-R 非常感谢
现在代码写到这了
rm(list=ls(all=T))
library(data.table)
evi <- fread("J:/x.txt",header = F)
soil <- fread("J:/y.txt",header = F)
library(randomForest)
coe_rf <- matrix(0,nrow=10321465,ncol=1)
evi_rf <- matrix(0,nrow=10321465,ncol=17)
column <- c(1:10321465)
for (i in column){
evi0 <- evi[,i,with=FALSE]
soil0 <- soil[,i,with=FALSE]
grid <-data.frame(evi0,soil0)
names(grid) <- c("evi","soil")
train <- grid[1:16,]
test <- grid[17,]
set.seed(1234)
fit <- randomForest(evisoil,data=train,na.action=na.roughfix)
fit_pred <- predict(fit,newdata=test)
evi_rf[1:16,i] <- fit$predicted
evi_rf[17,i] <- fit_pred
coe_rf[1,i] <- fit$rsq[500]
}
write.table(coe_rf, file ="J:/random_forest_coef.txt", sep ="", row.names =FALSE, col.names =FALSE, quote =FALSE)
write.table(evi_rf, file ="J:/random_forest_evi_predict.txt", sep ="", row.names =FALSE, col.names =FALSE, quote =FALSE)
循环里用小数据测试的时候没问题,但是用全部数据循环的时候显示我导出随机森林的解释率的时候出错,想问问 有没有导出过随机森林结果参数的呀?