#Run the PCA using the rda() function (NB: rda() is used for both PCA and RDA)
spe.h.pca <- rda(spe.hel)

#Extract the results
summary(spe.h.pca) #overall results