Population genetic theory predicts that plant populations will exhibit internal spatial autocorrelation when propagule flow is restricted, but as an empirical reality, spatial structure is rarely consistent across loci or sites, and is generally weak. A lack of sensitivity in the statistical procedures may explain the discrepancy. Most work to date, based on allozymes, has involved pattern analysis for individual alleles, but new PCR-based genetic markers are coming into vogue, with vastly increased numbers of alleles. The field is badly in need of an explicitly multivariate approach to autocorrelation analysis, and our purpose here is to introduce a new approach that is applicable to multiallelic codominant, multilocus arrays. The procedure treats the genetic data set as a whole, strengthening the spatial signal and reducing the stochastic (allele-to-allele, and locus-to-locus) noise. We (i) develop a very general multivariate method, based on genetic distance methods, (ii) illustrate it for multiallelic codominant loci, and (iii) provide nonparametric permutational testing procedures for the full correlogram. We illustrate the new method with an example data set from the orchid Caladenia tentaculata, for which we show (iv) how the multivariate treatment compares with the single-allele treatment, (v) that intermediate frequency alleles from highly polymorphic loci perform well and rare alleles poorly, (vi) that a multilocus treatment provides clearer answers than separate single-locus treatments, and (vii) that weighting alleles differentially improves our resolution minimally. The results, though specific to Caladenia, offer encouragement for wider application.