我在r中使用了mvpart包来生成多元回归树。为了确定负责拆分的指标物种,我使用了labdsv包中的indval命令。
我的问题是,如何将指标物种插入回归树?
这是我使用的脚本,
mvpart(回归树)
mfish<-mvpart(as.matrix(fisheries.hell[,2:34])~Time+Fishing.gear+
Fishing.vessel+Fishing.ground+
Fishing.hours+Fishing.days,fisheries.wide,
method="mrt",xv="pick",which="4",
xvmult=100,bars=F)
indval
mfish.indval$indcls[which(mfish.indval$pval<=0.05)]