?mvpart
 
#Preparer les données: enlever la variable “distance from source”
env <- subset(env, select = -das)

# Créer l'arbre
doubs.mrt<-mvpart(data.matrix(spe.hel)~.,env,
legend=FALSE,margin=0.01,cp=0,xv="pick",xval=nrow(spe.hel),xvmult=100,which=4)