#Run the PCA
env.pca <- rda(env.z) # or rda(env, scale=TRUE)

#Extract the results
summary(env.pca)
summary(env.pca, scaling=2)