我正在为随机森林绘制 R 中的一些部分依赖图。我正在使用 pdp 包,它对我在 z 变量中的限制没有任何规范。代码:
pd <- partial(RF, pred.var= c ( "SOC", "Q"))
plotPartial(pd, zlim=c(1000,3000))
我也尝试过:zlim(c(1000,3000)、zmax 和 min 以及类似的
感谢您的建议,这是我的代码
library(randomForest) #fitting RF
library(pdp) #partial dependence analysis
set.seed(100)
RF=randomForest(CorrYield~Year+SOC+pH+ Sorghum+Pbray + FertP+ FertN+ PP_veg +PP_PC, data= RF_sin1 , ntree=500, importance=T, proximity=T)
pd_SOC_P <- partial(RF, pred.var= c("SOC", "Pbray"))
plotPartial(pd_SOC_P,xlab="SOC(%)", ylab="PBray(mg kg^-1)")