Three hundred ancient Turkish potsherds were analyzed by instrumental neutron activation analysis, and the resulting data analyzed by several techniques of multivariate statistical analysis, some only recently developed. The programs AGCLUS, MASLOC, and SIMCA were sequentially employed to characterize and group the samples by type of pottery and by site of excavation or collection. Comparison of the statistical analyses by each methods provided archaeologica insight into the site/type relationships of the samples and ultimately evidence relevant to commercial relations between the ancient communities and to specialization of pottery production over time. The techniques used for statistical analysis were found to be of significant potential utility in the future analysis of other archaeometric data sets. -- CAL