Functional-space biplot with KDE density shading and phylopic image overlays.
PCA_plot_funspace <- function(PCoAPlot, PCoACorPlot, multAx1, multAx, title, legend = NULL, probs = c(0.25, 0.5, 0.99),