PCoA (classical multidimensional scaling) works from any distance matrix. PCA needs numeric data; missing values are replaced by the column mean, and traits (not markers) are scaled to unit variance.
Usage
ordinate(dist = NULL, d = NULL, method = c("pcoa", "pca"), n_axes = 5)