#Preparing the data prior to RDA env.z <- subset(env.z, select = -das) # remove the "distance from the source" variable #Running the RDA ?rda spe.rda <- rda(spe.hel~., data=env.z) ### Extract the results summary(spe.rda, display=NULL) #The results are called using summary: summary(spe.rda, display=NULL) #display = NULL optional