A PROVENANCE STUDY OF FRENCH LIMESTONE BASED ON VARIABLE SELECTION FROM COMPOSITIONAL PROFILES

The present study shows how multivariate analysis and variable selection techniques can be used in archaeological provenances studies to improve classification performances. Neutron activation analysis (NAA), in combination with stepwise linear discrimination analysis (SLDA) (capable of simultaneously performing variable selection and classification), was applied to differentiate among limestone samples from different quarries across the north of France based on their compositional fingerprint. A hierarchical classification approach was followed, aimed at progressively assigning limestone samples to more specific sources of origin, from the broadest classification units (French regions) to the narrowest ones (individual quarries). The application of the stepwise variable selection procedure to extract the most discriminating compositional variables prior to each classification development allowed us to obtain a perfect separation between the limestone categories considered at every classification stage (all samples were correctly classified and predicted in all cases). The high‐quality results obtained were even more remarkable considering the relatively small number of significant variables selected in each case using the SLDA method. An illustrative example was provided in order to demonstrate that the classification strategy proposed can actually allocate an unknown sculpture to a particular quarry of origin.     

THE ESTABLISHMENT OF A DATA BASE OF NEUTRON ACTIVATION ANALYSES OF WHITE MARBLE*

The establishment of a data base of trace element analysis results obtained by neutron activation analysis of 183 marble samples, each weighing about 200 mg, from eight different quarrying localities for use in provenance studies is described and a summary of the results presented. The problems of sampling museum artefacts are discussed, and details of tests of the reproducibility of the technique are given together with those obtained from multiple-samples from single museum objects. Replicate analyses of a marble standard are also discussed. Thus these data are a measure of the accuracy and precision of the technique, as well as allowing an assessment of the variation in trace element composition within quarry sources and single blocks of marble, as used for artefacts.     

MULTIVARIATE ANALYSIS OF COMPOSITIONAL DATA: APPLIED COMPARISONS FAVOUR STANDARD PRINCIPAL COMPONENTS ANALYSIS OVER AITCHISON'S LOGLINEAR CONTRAST METHOD

There has been debate about whether standard principal components analysis is appropriate for the multivariate analysis of compositional data (e.g. oxide composition of glass), Loglinear transformation has been recommended by Aitchison as a prerequisite. This paper argues that previous comparisons of methodological merits have tended to circularity of argument by making assumptions about the form of a good multivariate result. To break the circularity of argument the authors have introduced randomized variables into five data sets. A good result must recognize these randomized variables as noise and place them near the centroid of the principal components scattergram of variable loadings. Standard principal components analysis is found to perform better than loglinear transformation in its ability to recognize the randomized variables. It is concluded that loglinear transformation tends to introduce spurious structure into a table of compositional data. This paper is followed by a comment by M. J. Baxter.     

MODIFIED MAHALANOBIS FILTERS FOR GROUPING POTTERY BY CHEMICAL COMPOSITION*

A statistical procedure for grouping pottery in provenance studies by chemical data is presented, which now is routinely in use in our laboratory. It is based on the Mahalanobis filter method and X² -statistics, and can be used for both establishing groups and assigning single sherds to already known groups, thus replacing principal components analysis or cluster analysis and avoiding their problems in grouping pottery. The new method is able to consider correlations, uncertainties of measurement and constant shifts of the data in case of dilution. In particular, considering dilution effects results in both a better assignment to and separation of reference groups and is also equivalent to the compositional data approach, if log-transformed data are used. Other distortions of data (e.g., mixing of clays) can also be considered.     

ROMAN AND MEDIEVAL GLASS FROM THE ITALIAN AREA: BULK CHARACTERIZATION AND RELATIONSHIPS WITH PRODUCTION TECHNOLOGIES*

Compositional investigations were performed on 81 Roman and medieval glass fragments (first to 14th centuries ad ) from four Italian archaeological sites. The samples were soda–lime–silica in composition, with natron as flux for the Roman and early medieval glass samples, and with plant ash as flux for the late medieval ones. The varying colours are due to the differing FeO, Fe ₂ O ₃ , MnO and Sb ₂ O ₃ contents. Hierarchical cluster analysis identified six compositional groups related to age, which were compared with those found in the literature. In this way, technological continuity from the Roman to the early medieval period and the appearance of plant ash technology in the ninth century, 200 years in advance of the period previously believed, are demonstrated.     

MATHEMATICS,STATISTICS AND ARCHAEOMETRY: THE PAST 50 YEARS OR SO*

This review of developments in the use of mathematics and statistics in archaeometry over the past 50 years is partial, personal and ‘broad‐brush’. The view is expressed that it is in the past 30 years or so that the major developments have taken place. The view is also expressed that, with the exception of methods for analysing radiocarbon dates and increased computational power, mathematical and statistical methods that are currently used, and found to be useful in widespread areas of application such as provenance studies, don't differ fundamentally from what was being done 30 years ago.     

RQ-MODE PRINCIPAL COMPONENTS ANALYSIS OF CERAMIC COMPOSITIONAL DATA*

RQ-mode principal components analysis (PCA) is a means for calculating variable and object loadings on the same axes, so that elements can be displayed along with data points on a single diagram. The biplots resulting from RQ-mode PCA preserve both Euclidean relations among the objects and variance-covariance structure. When used with data on the chemical composition of archaeological pottery, such biplots facilitate recognizing compositional subgroups and determining the chemical basis of group separation. RQ-mode PCA is illustrated in this paper with neutron activation data on Mesoamerican Plumbate pottery.     

MEDIEVAL GLASS FROM ROCCA DI ASOLO (NORTHERN ITALY): AN ARCHAEOMETRIC STUDY

An archaeometric study was performed on 33 medieval glass samples from Rocca di Asolo (northern Italy), in order to study the raw materials employed in their production, identify analogies with medieval glass from the Mediterranean area and possible relationships between chemical composition and type and/or production technique, contextualize the various phases of the site and extend data on Italian medieval glass. The samples are soda–lime–silica in composition, with natron as flux for early medieval glasses and soda ash for the high and late medieval ones. Compositional groups were identified, consistent with the major compositional groups identified in the western Mediterranean during the first millennium AD . In particular, Asolo natron glass is consistent with the HIMT group and recycled Roman glass; soda ash glass was produced with the same type of flux (Levantine ash) but a different silica source (siliceous pebbles, and more or less pure sand). Cobalt was the colouring agent used to obtain blue glass; analytical data indicate that at least two different sources of Co were exploited during the late medieval period. Some data, analytical and historical, suggest a Venetian provenance for the high/late medieval glass and a relationship between type of object (beaker or bottle) and chemical composition.     

CHEMISTRY VERSUS DATA DISPERSION: IS THERE A BETTER WAY TO ASSESS AND INTERPRET ARCHAEOMETRIC DATA?

Given the common use of chemical concentration data to define ceramic groups that aid in the exploration of ancient technology, trade and provenance, it is important to reflect on how we collectively establish and define both chemical groups and outliers. In this paper, we argue that commonly used data analysis procedures, such as principal component analysis and centred log‐ratio principal component analysis favoured in the examination of ceramic chemical data, although rapid and easy, may overlook existing chemical groups and outliers, especially when the ratio of non‐diagnostic to diagnostic elements is high. To evaluate whether geochemistry is more important than data dispersion in data assessment, we re‐examine chemical concentration data from previously published ceramic, clay and daub samples from the lower Ohio River Valley. We begin by briefly discussing steps we took to ensure that the data set reflects geochemical differences, rather than analytical or data transfer errors. Next, we use bivariate plots, as well as PCA and CLR–PCA, to examine different versions of our altered data, using varying numbers of element combinations. We propose that the careful examination of bivariate plots is critical in establishing the elements that should be included in PCA and other multivariate analyses.     

A MULTIDISCIPLINARY APPROACH TO THE CHARACTERIZATION OF ROMAN GRAVESTONES FROM AQUILEIA (UDINE,ITALY)*

The results of an archaeometric study on gravestone production in Aquileia (northeastern Italy) during Roman Republican times are reported here. The artefacts are made of limestone, quarried from the sedimentary sequence of the Trieste Karst (Carnican Alps, northeastern Italy), as testified by historical sources which refer to the ‘Roman Quarry’ at Aurisina. The rock types used for the gravestones are petrographically, geochemically and geochronologically homogeneous, and display similar faunal assemblages. In addition, their age is restricted to the Santonian – Lower Campanian. Comparisons between archaeological finds and rock specimens from various active and inactive quarries in this area suggest that quarrying during Roman Republican times was limited to a specific stratigraphic horizon of the sedimentary sequence, corresponding to the higher part of the ‘Roman Quarry’.