Intrinsic test of independence in contingency tables

A condition needed for testing nested hypotheses from a Bayesianviewpoint is that the prior for the alternative model concentratesmass around the small, or null, model. For testing independencein contingency tables, the intrinsic priors satisfy this requirement.Further, the degree of concentration of the priors is controlled bya discrete parameter m, the training sample size, which plays animportant role in the resulting answer regardless of the samplesize.In this paper we study robustness of the tests of independencein contingency tables with respect to the intrinsic priors withdifferent degree of concentration around the null, and comparewith other “robust” results by Good and Crook. Consistency ofthe intrinsic Bayesian tests is established.We also discuss conditioning issues and sampling schemes,and argue that conditioning should be on either one margin orthe table total, but not on both margins.Examples using real are simulated data are given

Compositional Time Series: An Application

The composition of the labour force is an important economic factor for a country.Often the changes in proportions of different groups are of interest.I this paper we study a monthly compositional time series from the Swedish LabourForce Survey from 1994 to 2005. Three models are studied: the ILR-transformed series,the ILR-transformation of the compositional differenced series of order 1, and the ILRtransformationof the compositional differenced series of order 12. For each of thethree models a VAR-model is fitted based on the data 1994-2003. We predict the timeseries 15 steps ahead and calculate 95 % prediction regions. The predictions of thethree models are compared with actual values using MAD and MSE and the predictionregions are compared graphically in a ternary time series plot.We conclude that the first, and simplest, model possesses the best predictive power ofthe three models

Compositional evolution with mass transfer in closed systems

Evolution of compositions in time, space, temperature or other covariates is frequentin practice. For instance, the radioactive decomposition of a sample changes its composition with time. Some of the involved isotopes decompose into other isotopes of thesample, thus producing a transfer of mass from some components to other ones, butpreserving the total mass present in the system. This evolution is traditionally modelledas a system of ordinary di erential equations of the mass of each component. However,this kind of evolution can be decomposed into a compositional change, expressed interms of simplicial derivatives, and a mass evolution (constant in this example). A rst result is that the simplicial system of di erential equations is non-linear, despiteof some subcompositions behaving linearly.The goal is to study the characteristics of such simplicial systems of di erential equa-tions such as linearity and stability. This is performed extracting the compositional differential equations from the mass equations. Then, simplicial derivatives are expressedin coordinates of the simplex, thus reducing the problem to the standard theory ofsystems of di erential equations, including stability. The characterisation of stabilityof these non-linear systems relays on the linearisation of the system of di erential equations at the stationary point, if any. The eigenvelues of the linearised matrix and theassociated behaviour of the orbits are the main tools. For a three component system,these orbits can be plotted both in coordinates of the simplex or in a ternary diagram.A characterisation of processes with transfer of mass in closed systems in terms of stability is thus concluded. Two examples are presented for illustration, one of them is aradioactive decay

A Log-linear model of grain size influence on the geochemistry of sediments

Sediment composition is mainly controlled by the nature of the source rock(s), and chemical (weathering) and physical processes (mechanical crushing, abrasion, hydrodynamic sorting) during alteration and transport. Although the factors controlling these processes are conceptually well understood, detailed quantification of compositional changes induced by a single process are rare, as are examples where the effects of several processes can be distinguished. The present study was designed to characterize the role of mechanical crushing and sorting in the absence of chemical weathering. Twenty sediment samples were taken from Alpine glaciers that erode almost pure granitoid lithologies. For each sample, 11 grain-size fractions from granules to clay (ø grades &-1 to &9) were separated, and each fraction was analysed for its chemical composition.The presence of clear steps in the box-plots of all parts (in adequate ilr and clr scales) against ø is assumed to be explained by typical crystal size ranges for the relevant mineral phases. These scatter plots and the biplot suggest a splitting of the full grain size range into three groups: coarser than ø=4 (comparatively rich in SiO2, Na2O, K2O, Al2O3, and dominated by “felsic” minerals like quartz and feldspar), finer than ø=8 (comparatively rich in TiO2, MnO, MgO, Fe2O3, mostly related to “mafic” sheet silicates like biotite and chlorite), and intermediate grains sizes (4≤ø &8; comparatively rich in P2O5 and CaO, related to apatite, some feldspar).To further test the absence of chemical weathering, the observed compositions were regressed against three explanatory variables: a trend on grain size in ø scale, a step function for ø≥4, and another for ø≥8. The original hypothesis was that the trend could be identified with weathering effects, whereas each step function would highlight those minerals with biggest characteristic size at its lower end. Results suggest that this assumption is reasonable for the step function, but that besides weathering some other factors (different mechanical behavior of minerals) have also an important contribution to the trend.Key words: sediment, geochemistry, grain size, regression, step function

The cover basement contact beneath the Rawil axial depression (Western Alps) - True amplitude seismic processing, petrophysical properties, and modeling

Since 1986, several near-vertical seismic reflection profiles have been recorded in Switzerland in order to map the deep geologic structure of the Alps. One objective of this endeavour has been to determine the geometries of the autochthonous basement and of the external crystalline massifs, important elements for understanding the geodynamics of the Alpine orogeny. The PNR-20 seismic line W1, located in the Rawil depression of the western Swiss Alps, provides important information on this subject. It extends northward from the `'Penninic front'' across the Helvetic nappes to the Prealps. The crystalline massifs do not outcrop along this profile. Thus, the interpretation of `'near-basement'' reflections has to be constrained by down-dip projections of surface geology, `'true amplitude'' processing, rock physical property studies and modelling. 3-D seismic modelling has been used to evaluate the seismic response of two alternative down-dip projection models. To constrain the interpretation in the southern part of the profile, `'true amplitude'' processing has provided information on the strength of the reflections. Density and velocity measurements on core samples collected up-dip from the region of the seismic line have been used to evaluate reflection coefficients of typical lithologic boundaries in the region. The cover-basement contact itself is not a source of strong reflections, but strong reflections arise from within the overlaying metasedimentary cover sequence, allowing the geometry of the top of the basement to be determined on the basis of `'near-basement'' reflections. The front of the external crystalline massifs is shown to extend beneath the Prealps, about 6 km north of the expected position. A 2-D model whose seismic response shows reflection patterns very similar to the observed is proposed.

The TRANSMED Atlas: geological-geophysical fabric of the Mediterranean region - Final report of the project

Geological research on the Mediterranean region is presently characterized by the transition from disciplinary to multidisciplinary research, as well as from national to international investigations. In order to synthesize and integrate the vast disciplinary and national datasets which are available, it is necessary to implement maximum interaction among geoscientists of different backgrounds. The creation of project-oriented task forces in universities and other research institutions, as well as the development of large international cooperation programs, is instrumental in pursuing such a multidisciplinary and supranational approach. The TRANSMED Atlas, an official publication of the 32nd International Geological Congress (Florence 2004), is the result of an international scientific cooperation program which brought together for over two years sixty-three structural geologists, geophysicists, marine geologists, petrologists, sedimentologists, stratigraphers, paleogeographers, and petroleum geologists coming from eighteen countries, and working for the petroleum industry, academia, and other institutions, both public and private. The TRANSMED Atlas provides an updated, synthetic, and coherent portrayal of the overall geological-geophysical structure of the Mediterranean domain and the surrounding areas. The initial stimulus for the Atlas came from the realization of the extremely heterogeneous nature of the existing geological-geophysical data about such domain. These data have been gathered by universities, oil companies, geological surveys and other institutions in several countries, often using different procedures and standards. In addition, much of these data are written in languages and published in outlets that are not readily accessible to the general international reader. By synthesizing and integrating a wealth of preexisting and new data derived from surficial geology, seismic sections at various scales, and mantle tomographies, the TRANSMED Atlas provides for the first time a coherent geological overview of the Mediterranean region and represents an ideal springboard for future studies.

The hypogene iron oxide copper-gold mineralization in the Mantoverde District, Northern Chile

The Mantoverde iron oxide copper-gold (IOCC) district, northern Chile, is known for its Cu production from supergene ores. Recently, exploration outlined an additional hypogene ore resource of 440 Mt with 0.56 percent Cu, and 0.12 g/t An. The hypogene sulfide mineralization occurs mainly as chalcopyrite and pyrite, typically in specularite or magnetite-cemented breccias and associated stockworks. The host rocks underwent variably intense K feldspar alteration, chloritization, sericitization, silicification, and/or carbonatization. A district scale Na(-Ca) alteration is absent. The IOCC mineralization in the district shows a strong tectonic control by northwest- to north-northwest-trending brittle structures. Large Cu sulfide-rich veins or Cu sulfide-cemented breccias are absent. Therefore, head grades of 4 percent Cu are an exception. There is a positive correlation between Cu and An grades. Gold is probably contained mostly in chalcopyrite and pyrite. Elevated concentrations of light rare-earth elements (LREE) occur locally but are attributed to redistribution of LREE within the deposits rather than to derivation from external sources. The Cu-Au ores in the Mantoverde district are low in and have relatively low contents in heavy metals that are potentially hazardous to the environment, such as As (avg 14 ppm), Hg (<5 ppm), or Cd (<0.2 ppm). The sulfur isotope ratios of chalcopyrite from the IOCC deposits lie between -5.6 and 8.9 per mil delta(34)S(VCDT). They show systematic variations within the district, which are interpreted to reflect relative distance to inferred fluid conduits and the level of deposition within the hydrothermal system. Most initial (87)Sr/(86)Sr values of altered volcanic rocks and hydrothermal calcite from the Mantoverde district are between 0.7031 and 0.7060 and are similar to those of the igneous rocks of the region. Lead isotope ratios of chalcopyrite are consistent with Pb (and by inference Cu) derived from Early Cretaceous magmatism. The sulfur, strontium, and lead isotope data of chalcopyrite, calcite gangue, or altered host rocks, respectively, are compatible with a genetic model that involves cooling of metal and sulfur-bearing magmatic-hydrothermal fluids that mix with meteoric waters or seawater at relatively shallow crustal levels. An additional exotic sulfur input is likely, though not required, for the copper mineralization. Apart from the IOCC. deposits, there are a number of smaller magnetite(-apatite) bodies in the district. These are geologically similar to the Cu-Au-bearing magnetite bodies, but are related to splays of the north-south-trending Atacama fault zone and differ in alteration and texture.

Uncertainty Quantification Of A Semi-Supervised Support Vector Regression Reservoir Model

Uncertainty quantification of petroleum reservoir models is one of the present challenges, which is usually approached with a wide range of geostatistical tools linked with statistical optimisation or/and inference algorithms. Recent advances in machine learning offer a novel approach to model spatial distribution of petrophysical properties in complex reservoirs alternative to geostatistics. The approach is based of semisupervised learning, which handles both ?labelled? observed data and ?unlabelled? data, which have no measured value but describe prior knowledge and other relevant data in forms of manifolds in the input space where the modelled property is continuous. Proposed semi-supervised Support Vector Regression (SVR) model has demonstrated its capability to represent realistic geological features and describe stochastic variability and non-uniqueness of spatial properties. On the other hand, it is able to capture and preserve key spatial dependencies such as connectivity of high permeability geo-bodies, which is often difficult in contemporary petroleum reservoir studies. Semi-supervised SVR as a data driven algorithm is designed to integrate various kind of conditioning information and learn dependences from it. The semi-supervised SVR model is able to balance signal/noise levels and control the prior belief in available data. In this work, stochastic semi-supervised SVR geomodel is integrated into Bayesian framework to quantify uncertainty of reservoir production with multiple models fitted to past dynamic observations (production history). Multiple history matched models are obtained using stochastic sampling and/or MCMC-based inference algorithms, which evaluate posterior probability distribution. Uncertainty of the model is described by posterior probability of the model parameters that represent key geological properties: spatial correlation size, continuity strength, smoothness/variability of spatial property distribution. The developed approach is illustrated with a fluvial reservoir case. The resulting probabilistic production forecasts are described by uncertainty envelopes. The paper compares the performance of the models with different combinations of unknown parameters and discusses sensitivity issues.

Sulfur and strontium isotope composition of the Llobregat River (NE Spain): Tracers of natural and anthropogenic chemicals in stream waters

The use of sulfur and strontium isotopes as tracers for the source/s of water contaminants have been applied to the water of the Llobregat River system (NE Spain). Surface water samples from June 1997 were collected from the Llobregat River and its main tributaries and creeks. The chemistry of most stream waters are controlled mainly by the weathering of Tertiary chemical sediments within the drainage basin. The largest variation in delta(34)S values were found in the small creeks with values ranging from -9.9 to 15parts per thousand, whilst in the main river channels values ranged from 6.3 to 12.4parts per thousand. The Sr-87/Sr-86 ratio for dissolved strontium ranged from 0.70795 for a non-polluted site to 0.70882 for a polluted one. Most of the waters with high NO3 and low Ca/Na ratio converge to the same Sr-87/Sr-86 value, pointing to dominant pollutant end member contribution or a mixing of pollutants with an isotopic composition around 0.7083-0.7085. Although the concentration of the natural inputs in the river for sulfate and strontium are high, as a result of the sulfate outcrops within the geology of the basin, their isotopic characteristics suggest that they can be used as a discriminating device in water pollution problems. However to establish the detailed characteristics of the isotopes as geochemical tools, specific high-resolution case studies are necessary in small areas, where the inputs are well known.

Markov chain montecarlo method applied to rounding zeros of compositional data: first approach

As stated in Aitchison (1986), a proper study of relative variation in a compositional data set should be based on logratios, and dealing with logratios excludes dealing with zeros. Nevertheless, it is clear that zero observations might be present in real data sets, either because the corresponding part is completelyabsent –essential zeros– or because it is below detection limit –rounded zeros. Because the second kind of zeros is usually understood as “a trace too small to measure”, it seems reasonable to replace them by a suitable small value, and this has been the traditional approach. As stated, e.g. by Tauber (1999) and byMartín-Fernández, Barceló-Vidal, and Pawlowsky-Glahn (2000), the principal problem in compositional data analysis is related to rounded zeros. One should be careful to use a replacement strategy that does not seriously distort the general structure of the data. In particular, the covariance structure of the involvedparts –and thus the metric properties– should be preserved, as otherwise further analysis on subpopulations could be misleading. Following this point of view, a non-parametric imputation method isintroduced in Martín-Fernández, Barceló-Vidal, and Pawlowsky-Glahn (2000). This method is analyzed in depth by Martín-Fernández, Barceló-Vidal, and Pawlowsky-Glahn (2003) where it is shown that thetheoretical drawbacks of the additive zero replacement method proposed in Aitchison (1986) can be overcome using a new multiplicative approach on the non-zero parts of a composition. The new approachhas reasonable properties from a compositional point of view. In particular, it is “natural” in the sense thatit recovers the “true” composition if replacement values are identical to the missing values, and it is coherent with the basic operations on the simplex. This coherence implies that the covariance structure of subcompositions with no zeros is preserved. As a generalization of the multiplicative replacement, in thesame paper a substitution method for missing values on compositional data sets is introduced

Modelling of Mercury’s surface composition and remote detection from the orbit with the BepiColombo Mercury Planetary Orbiter

It can be assumed that the composition of Mercury’s thin gas envelope (exosphere) is related to thecomposition of the planets crustal materials. If this relationship is true, then inferences regarding the bulkchemistry of the planet might be made from a thorough exospheric study. The most vexing of allunsolved problems is the uncertainty in the source of each component. Historically, it has been believedthat H and He come primarily from the solar wind, while Na and K originate from volatilized materialspartitioned between Mercury’s crust and meteoritic impactors. The processes that eject atoms andmolecules into the exosphere of Mercury are generally considered to be thermal vaporization, photonstimulateddesorption (PSD), impact vaporization, and ion sputtering. Each of these processes has its owntemporal and spatial dependence. The exosphere is strongly influenced by Mercury’s highly ellipticalorbit and rapid orbital speed. As a consequence the surface undergoes large fluctuations in temperatureand experiences differences of insolation with longitude. We will discuss these processes but focus moreon the expected surface composition and solar wind particle sputtering which releases material like Caand other elements from the surface minerals and discuss the relevance of composition modelling

alr approach for replacing values below the detection limit

All of the imputation techniques usually applied for replacing values below thedetection limit in compositional data sets have adverse effects on the variability. In thiswork we propose a modification of the EM algorithm that is applied using the additivelog-ratio transformation. This new strategy is applied to a compositional data set and theresults are compared with the usual imputation techniques

CoDaPack. An Excel and Visual Basic based software of compositional data analysis. Current version and discussion form upcoming versions

In the eighties, John Aitchison (1986) developed a new methodological approach for the statistical analysis of compositional data. This new methodology was implemented in Basic routines grouped under the name CODA and later NEWCODA inMatlab (Aitchison, 1997). After that, several other authors have published extensions to this methodology: Marín-Fernández and others (2000), Barceló-Vidal and others (2001), Pawlowsky-Glahn and Egozcue (2001, 2002) and Egozcue and others (2003). (...)

Major-elements trends in cenozoic volcanites of Hungary

Hungary lies entirely within the Carpatho-Pannonian Region (CPR), a dominant tectonic unit of eastern Central Europe. The CPR consists of the Pannonian Basin system, and the arc of the Carpathian Mountains surrounding the lowlands in the north, east, and southeast. In the west, the CPR is bounded by the Eastern Alps, whereas in the south, by the Dinaridic belt. (...)

Bayesian tools for zero counts in compositional data

The log-ratio methodology makes available powerful tools for analyzing compositionaldata. Nevertheless, the use of this methodology is only possible for those data setswithout null values. Consequently, in those data sets where the zeros are present, aprevious treatment becomes necessary. Last advances in the treatment of compositionalzeros have been centered especially in the zeros of structural nature and in the roundedzeros. These tools do not contemplate the particular case of count compositional datasets with null values. In this work we deal with \count zeros" and we introduce atreatment based on a mixed Bayesian-multiplicative estimation. We use the Dirichletprobability distribution as a prior and we estimate the posterior probabilities. Then weapply a multiplicative modi¯cation for the non-zero values. We present a case studywhere this new methodology is applied.Key words: count data, multiplicative replacement, composition, log-ratio analysis