An R Library for Compositional Data Analysis in Archaeometry

Compositional data naturally arises from the scientific analysis of the chemical composition of archaeological material such as ceramic and glass artefacts. Data of this type can be explored using a variety of techniques, from standard multivariate methods such as principal components analysis and cluster analysis, to methods based upon the use of log-ratios. The general aim is to identify groups of chemically similar artefacts that could potentially be used to answer questions of provenance. This paper will demonstrate work in progress on the development of a documented library of methods, implemented using the statistical package R, for the analysis of compositional data. R is an open source package that makes available very powerful statistical facilities at no cost. We aim to show how, with the aid of statistical software such as R, traditional exploratory multivariate analysis can easily be used alongside, or in combination with, specialist techniques of compositional data analysis. The library has been developed from a core of basic R functionality, together with purpose-written routines arising from our own research (for example that reported at CoDaWork'03). In addition, we have included other appropriate publicly available techniques and libraries that have been implemented in R by other authors. Available functions range from standard multivariate techniques through to various approaches to log-ratio analysis and zero replacement. We also discuss and demonstrate a small selection of relatively new techniques that have hitherto been little-used in archaeometric applications involving compositional data. The application of the library to the analysis of data arising in archaeometry will be demonstrated; results from different analyses will be compared; and the utility of the various methods discussed

Compositional analysis of archaeological glasses

At CoDaWork'03 we presented work on the analysis of archaeological glass composi- tional data. Such data typically consist of geochemical compositions involving 10-12 variables and approximates completely compositional data if the main component, sil- ica, is included. We suggested that what has been termed `crude' principal component analysis (PCA) of standardized data often identi ed interpretable pattern in the data more readily than analyses based on log-ratio transformed data (LRA). The funda- mental problem is that, in LRA, minor oxides with high relative variation, that may not be structure carrying, can dominate an analysis and obscure pattern associated with variables present at higher absolute levels. We investigate this further using sub- compositional data relating to archaeological glasses found on Israeli sites. A simple model for glass-making is that it is based on a `recipe' consisting of two `ingredients', sand and a source of soda. Our analysis focuses on the sub-composition of components associated with the sand source. A `crude' PCA of standardized data shows two clear compositional groups that can be interpreted in terms of di erent recipes being used at di erent periods, re ected in absolute di erences in the composition. LRA analysis can be undertaken either by normalizing the data or de ning a `residual'. In either case, after some `tuning', these groups are recovered. The results from the normalized LRA are di erently interpreted as showing that the source of sand used to make the glass di ered. These results are complementary. One relates to the recipe used. The other relates to the composition (and presumed sources) of one of the ingredients. It seems to be axiomatic in some expositions of LRA that statistical analysis of compositional data should focus on relative variation via the use of ratios. Our analysis suggests that absolute di erences can also be informative

A compositional approach to stable isotope data analysis

Isotopic data are currently becoming an important source of information regarding sources, evolution and mixing processes of water in hydrogeologic systems. However, it is not clear how to treat with statistics the geochemical data and the isotopic data together. We propose to introduce the isotopic information as new parts, and apply compositional data analysis with the resulting increased composition. Results are equivalent to downscale the classical isotopic delta variables, because they are already relative (as needed in the compositional framework) and isotopic variations are almost always very small. This methodology is illustrated and tested with the study of the Llobregat River Basin (Barcelona, NE Spain), where it is shown that, though very small, isotopic variations comp lement geochemical principal components, and help in the better identification of pollution sources

Elemental analysis of New Zealand Nephrite Jade (pounamu) for forensic purposes

Pounamu (NZ jade), or nephrite, is a protected mineral in its natural form following the transfer of ownership back to Ngai Tahu under the Ngai Tahu (Pounamu Vesting) Act 1997. Any theft of nephrite is prosecutable under the Crimes Act 1961. Scientific evidence is essential in cases where origin is disputed. A robust method for discrimination of this material through the use of elemental analysis and compositional data analysis is required. Initial studies have characterised the variability within a given nephrite source. This has included investigation of both in situ outcrops and alluvial material. Methods for the discrimination of two geographically close nephrite sources are being developed. Key Words: forensic, jade, nephrite, laser ablation, inductively coupled plasma mass spectrometry, multivariate analysis, elemental analysis, compositional data analysis

A hybrid method for sag source location in power network

The work presented in this paper belongs to the power quality knowledge area and deals with the voltage sags in power transmission and distribution systems. Propagating throughout the power network, voltage sags can cause plenty of problems for domestic and industrial loads that can financially cost a lot. To impose penalties to responsible party and to improve monitoring and mitigation strategies, sags must be located in the power network. With such a worthwhile objective, this paper comes up with a new method for associating a sag waveform with its origin in transmission and distribution networks. It solves this problem through developing hybrid methods which hire multiway principal component analysis (MPCA) as a dimension reduction tool. MPCA reexpresses sag waveforms in a new subspace just in a few scores. We train some well-known classifiers with these scores and exploit them for classification of future sags. The capabilities of the proposed method for dimension reduction and classification are examined using the real data gathered from three substations in Catalonia, Spain. The obtained classification rates certify the goodness and powerfulness of the developed hybrid methods as brand-new tools for sag classification

Response of the benthic macroinvertebrate community to a point source in La Tordera stream (Catalonia, NE Spain)

Vam monitoritzar paràmetres físics i químics, macroinvertebrats bentònics, clorofil·la a, productors primaris i matèria orgànica durant un any (2001-2002) per examinar els efectes d'una font puntual sobre la composició taxonòmica, la estructura de la comunitat, l'organització funcional, la utilització de l'habitat i la estoquiometria al riu la Tordera (Catalunya). Aigües avall de la font puntual, concentració de nutrients, cabal i conductivitat eren majors que al tram d'aigües amunt, mentre que oxigen dissolt era menor. La densitat de macroinvertebrats era més elevada al tram d'aigües avall però la biomassa era similar als dos trams. La riquesa taxonòmica al tram de dalt era un 20% més alt que al tram de baix. Els anàlisis d'ordenació separen clarament els dos trams en el primer eix, mentre que els dos trams presentaven una pauta temporal similar en el segon eix. La similaritat entre els dos trams en composició taxonòmica, densitats i biomasses després de les crescudes d'abril i maig de 2002, indiquen que les pertorbacions del cabal poden actuar com a un mecanisme de reinici de la comunitat bentònica i jugar un paper important per a la restauració d'ecosistemes fluvials. Els dos trams presentaven una biomassa de perifiton, plantes vasculars, CPOM i FPOM similars, mentre que clorofil·la a, algues filamentoses, molses i SPOM eren majors al tram d'aigües avall. La densitat relativa de trituradors era menor sota la font puntual mentre que col·lectors i filtradors van ser afavorits. La biomassa relativa de trituradors també era menor sota la font puntual, però la biomassa de col·lectors i depredadors va augmentar. Les relacions entre densitat de grups tròfics i els seus recursos eren rarament significatives. La relació s'explicava millor amb la biomassa de macroinvertebrats. Els dos trams compartien la mateixa relació per raspadors, col·lectors i filtradors però no per trituradors i depredadors. La densitat i la biomassa de macroinvertebrats es trobaven positivament correlacionades amb la quantitat de recursos tròfics i la complexitat d'habitat, mentre que la riquesa taxonòmica es trobava negativament relacionada amb paràmetres hidràulics. La influència dels substrats inorgànics prenia menor rellevància per a la distribució dels macroinvertebrats. Els anàlisis d'ordenació mostren com les variables del microhabitat de major rellevància eren CPOM, clorofil·la a, algues filamentoses i velocitat. La cobertura de sorra només era significativa per al tram d'aigües amunt i les molses, al d'aigües avall. El número de correlacions significatives entre macroinvertebrats i les variables del microhabitat era més elevat per al tram de dalt que per al de baix, bàsicament per diferències en composició taxonòmica. La biomassa de macroinvertebrats va aportar una informació semblant a la obtinguda per la densitat. Perifiton i molses tenien uns continguts de nutrients similars en els dos trams. Els %C i %N d'algues filamentoses també eren similars en els dos trams però el %P sota la font puntual era el doble que al tram de dalt. Les relacions estoquiomètriques en CPOM, FPOM i SPOM eren considerablement menors sota la font puntual. Els continguts elementals i relacions van ser molt variables entre taxons de macroinvertebrats però no van resultar significativament diferents entre els dos trams. Dípters, tricòpters i efemeròpters presentaven una estoquiometria similar, mentre que el C i el N eren inferiors en moluscs i el P en coleòpters. Els depredadors presentaven un contingut en C i N més elevat que la resta de grups tròfics, mentre que el P era major en els filtradors. Els desequilibris elementals entre consumidors i recursos eren menors en el tram d'aigües avall. A la tardor i l'hivern la major font de nutrients va ser la BOM mentre que a la primavera i a l'estiu va ser el perifiton.