Possible solution of some essential zero problems in compositional data analysis

One of the tantalising remaining problems in compositional data analysis lies in how to deal with data sets in which there are components which are essential zeros. By anessential zero we mean a component which is truly zero, not something recorded as zero simply because the experimental design or the measuring instrument has not been sufficiently sensitive to detect a trace of the part. Such essential zeros occur inmany compositional situations, such as household budget patterns, time budgets,palaeontological zonation studies, ecological abundance studies. Devices such as nonzero replacement and amalgamation are almost invariably ad hoc and unsuccessful insuch situations. From consideration of such examples it seems sensible to build up amodel in two stages, the first determining where the zeros will occur and the secondhow the unit available is distributed among the non-zero parts. In this paper we suggest two such models, an independent binomial conditional logistic normal model and a hierarchical dependent binomial conditional logistic normal model. The compositional data in such modelling consist of an incidence matrix and a conditional compositional matrix. Interesting statistical problems arise, such as the question of estimability of parameters, the nature of the computational process for the estimation of both the incidence and compositional parameters caused by the complexity of the subcompositional structure, the formation of meaningful hypotheses, and the devising of suitable testing methodology within a lattice of such essential zero-compositional hypotheses. The methodology is illustrated by application to both simulated and real compositional data

Efectos de la aplicación de cenizas de caldera de biomasa en el modelo jerárquico de agregación de un suelo forestal bajo condiciones oceánicas

Se ha estudiado los efectos que tiene la aplicación de cenizas de caldera de biomasa, en el modelo jerárquico de agregación (Tisdall y Oades, 1982) y en la estabilización del C orgánico en un suelo forestal situado en la zona templada del País Vasco. Para ello, se aplicaron 3 tratamientos con diferentes dosis de ceniza en muestras de suelo tamizadas a 250 μm procedentes de un huerto semillero de pinus radiata. Estas muestras fueron incubadas durante 3 meses y fraccionadas los días 29, 44, 64, 78 y 92 con la intención de separar los macroagregados grandes (LMagg), de los macroagregados pequeños (Magg), microagregados (magg), limos (silt) y arcillas (clay). Todas las fracciones fueron analizadas para determinar su contenido en C orgánico. Los resultados mostraron que con la aplicación de cenizas, el modelo jerárquico de agregación de suelo se cumple para las fracciones LMagg, Magg y magg, pero no para la fracción silt. Además se ha observado que las cenizas promueven la formación de microagregados aumentado así la capacidad de secuestro de C del suelo, pero en cambio disminuye la proporción de macroagregados, hecho que podría acarrear una disminución en la calidad estructural del suelo.

Efectos de la aplicación de biochar en el modelo jerárquico de agregación de un suelo forestal bajo condiciones oceánicas

El objetivo del presente proyecto ha sido estudiar los efectos del biochar en el modelo jerárquico de agregación del suelo. Para estudiar estas propiedades se ha realizado una incubación de tres tratamientos de biochar Miscanthus (biochar 10t/ha (B10), 2t/ha (B2) y 10t/ha+Nitrógeno (BN), durante 91 días, realizando fraccionamientos en húmedo de las muestras los días 28, 43, 63, 77 y 91. En estos fraccionamientos se han separado las fracciones macroagregrado (LMagg+Magg), microagregados (magg), limos (Silt) y arcillas (Clay). Por último se determinó el C orgánico de todas las fracciones de los tratamientos obtenidos con el objeto de obtener información sobre el efecto que puede tener la aplicación del biochar sobre la estabilización de C en las fracciones del suelo. Los resultados obtenidos indican que con la aplicación de biochar, la teoría jerárquica de agregación propuesta por Tisdall y Oades (1982), no se cumple en la dinámica de agregación del suelo, si bien se puede observar una jerarquía en la estabilización del C orgánico en los macroagregados y microagregados. Además las dosis altas de biochar producen mayor cantidad de microagregados, favoreciendo el secuestro de C en el suelo. Por otro lado, con dosis altas de biochar se aprecia una disminución de la fracción de macroagregados, hecho que puede influir negativamente en la estructura del suelo.

Ab initio benchmark study for the oxidative addition of CH4 to Pd: importance of basis-set flexibility and polarization

To obtain a state-of-the-art benchmark potential energy surface (PES) for the archetypal oxidative addition of the methane C-H bond to the palladium atom, we have explored this PES using a hierarchical series of ab initio methods (Hartree-Fock, second-order Møller-Plesset perturbation theory, fourth-order Møller-Plesset perturbation theory with single, double and quadruple excitations, coupled cluster theory with single and double excitations (CCSD), and with triple excitations treated perturbatively [CCSD(T)]) and hybrid density functional theory using the B3LYP functional, in combination with a hierarchical series of ten Gaussian-type basis sets, up to g polarization. Relativistic effects are taken into account either through a relativistic effective core potential for palladium or through a full four-component all-electron approach. Counterpoise corrected relative energies of stationary points are converged to within 0.1-0.2 kcal/mol as a function of the basis-set size. Our best estimate of kinetic and thermodynamic parameters is -8.1 (-8.3) kcal/mol for the formation of the reactant complex, 5.8 (3.1) kcal/mol for the activation energy relative to the separate reactants, and 0.8 (-1.2) kcal/mol for the reaction energy (zero-point vibrational energy-corrected values in parentheses). This agrees well with available experimental data. Our work highlights the importance of sufficient higher angular momentum polarization functions, f and g, for correctly describing metal-d-electron correlation and, thus, for obtaining reliable relative energies. We show that standard basis sets, such as LANL2DZ+ 1f for palladium, are not sufficiently polarized for this purpose and lead to erroneous CCSD(T) results. B3LYP is associated with smaller basis set superposition errors and shows faster convergence with basis-set size but yields relative energies (in particular, a reaction barrier) that are ca. 3.5 kcal/mol higher than the corresponding CCSD(T) values

Gender inequality in health among elderly people in a combined framework of socioeconomic position, family characteristics and social support

This study analyses gender inequalities in health among elderly people in Catalonia (Spain) by adopting a conceptual framework that globally considers three dimensions of health determinants : socio-economic position, family characteristics and social support. Data came from the 2006 Catalonian Health Survey. For the purposes of this study a sub-sample of people aged 65–85 years with no paid job was selected (1,113 men and 1,484 women). The health outcomes analysed were self-perceived health status, poor mental health status and long-standing limiting illness. Multiple logistic regression models separated by sex were fitted and a hierarchical model was fitted in three steps. Health status among elderly women was poorer than among the men for the three outcomes analysed. Whereas living with disabled people was positively related to the three health outcomes and confidant social support was negatively associated with all of them in both sexes, there were gender differences in other social determinants of health. Our results emphasise the importance of using an integrated approach for the analysis of health inequalities among elderly people, simultaneously considering socio-economic position, family characteristics and social support, as well as different health indicators, in order fully to understand the social determinants of the health status of older men and women.

Assessment and reuse of contents in the competence-based educational platform InterMediActor

This paper describes a failure alert system and a methodology for content reuse in a new instructional design system called InterMediActor (IMA). IMA provides an environment for instructional content design, production and reuse, and for students’ evaluation based in content specification through a hierarchical structure of competences. The student assessment process and information extraction process for content reuse are explained.

GeneID in "Drosophila"

GeneID is a program to predict genes in anonymous genomic sequences designed with a hierarchical structure. In the first step, splice sites, and start and stop codons are predicted and scored along the sequence using position weight matrices (PWMs). In the second step, exons are built from the sites. Exons are scored as the sum of the scores of the defining sites, plus the log-likelihood ratio of a Markov model for coding DNA. In the last step, from the set of predicted exons, the gene structure is assembled, maximizing the sum of the scores of the assembled exons. In this paper we describe the obtention of PWMs for sites, and the Markov model of coding DNA in Drosophila melanogaster. We also compare other models of coding DNA with the Markov model. Finally, we present and discuss the results obtained when GeneID is used to predict genes in the Adh region. These results show that the accuracy of GeneID predictions compares currently with that of other existing tools but that GeneID is likely to be more efficient in terms of speed and memory usage.

Incorporating features of distribution and progression for automatic Makam classification

Automatic classification of makams from symbolic data is a rarely studied topic. In this paper, first a review of an n-gram based approach is presented using various representations of the symbolic data. While a high degree of precision can be obtained, confusion happens mainly for makams using (almost) the same scale and pitch hierarchy but differ in overall melodic progression, seyir. To further improve the system, first n-gram based classification is tested for various sections of the piece to take into account a feature of the seyir that melodic progression starts in a certain region of the scale. In a second test, a hierarchical classification structure is designed which uses n-grams and seyir features in different levels to further improve the system.

Drug Utilization Studies and Data Registries in Primary Care.

This article reviews the methodology of the studies on drug utilization with particular emphasis on primary care. Population based studies of drug inappropriateness can be done with microdata from Health Electronic Records and e-prescriptions. Multilevel models estimate the influence of factors affecting the appropriateness of drug prescription at different hierarchical levels: patient, doctor, health care organization and regulatory environment. Work by the GIUMAP suggest that patient characteristics are the most important factor in the appropriateness of prescriptions with significant effects at the general practicioner level.

Paral·lelització especulativa de programes sobre un clúster

Estudi sobre la millora de rendiment (en temps d’execució) al'algorisme de gràfics Fast Multipath Radiosity Using Hierarchical Subscenes gràcies a l’execució paral•lela especulada que ens permet obtenir el motor d'especulació per a clústers desenvolupat en el grup de recerca BCDS de la Universitat de Girona

Emotion and reason in everyday risk perception

Although research has documented the importance of emotion in risk perception, little is knownabout its prevalence in everyday life. Using the Experience Sampling Method, 94 part-timestudents were prompted at random via cellular telephones to report on mood state and threeemotions and to assess risk on thirty occasions during their working hours. The emotions valence, arousal, and dominance were measured using self-assessment manikins (Bradley &Lang, 1994). Hierarchical linear models (HLM) revealed that mood state and emotions explainedsignificant variance in risk perception. In addition, valence and arousal accounted for varianceover and above reason (measured by severity and possibility of risks). Six risks were reassessedin a post-experimental session and found to be lower than their real-time counterparts.The study demonstrates the feasibility and value of collecting representative samples of data withsimple technology. Evidence for the statistical consistency of the HLM estimates is provided inan Appendix.

A simple permutation test for clusteredness

Hierarchical clustering is a popular method for finding structure in multivariate data,resulting in a binary tree constructed on the particular objects of the study, usually samplingunits. The user faces the decision where to cut the binary tree in order to determine the numberof clusters to interpret and there are various ad hoc rules for arriving at a decision. A simplepermutation test is presented that diagnoses whether non-random levels of clustering are presentin the set of objects and, if so, indicates the specific level at which the tree can be cut. The test isvalidated against random matrices to verify the type I error probability and a power study isperformed on data sets with known clusteredness to study the type II error.

A multilevel analysis on the determinants of regional health care expenditure. A note

We apply a multilevel hierarchical model to explore whether anaggregation fallacy exists in estimating the income elasticity of healthexpenditure by ignoring the regional composition of national healthexpenditure figures. We use data for 110 regions in eight OECD countriesin 1997: Australia, Canada, France, Germany, Italy, Spain, Sweden andUnited Kingdom. In doing this we have tried to identify two sources ofrandom variation: within countries and between-countries. Our resultsshow that: 1- Variability between countries amounts to (SD) 0.5433, andjust 13% of that can be attributed to income elasticity and the remaining87% to autonomous health expenditure; 2- Within countries, variabilityamounts to (SD) 1.0249; and 3- The intra-class correlation is 0.5300. Weconclude that we have to take into account the degree of fiscaldecentralisation within countries in estimating income elasticity ofhealth expenditure. Two reasons lie behind this: a) where there isdecentralisation to the regions, policies aimed at emulating diversitytend to increase national health care expenditure; and b) without fiscaldecentralisation, central monitoring of finance tends to reduce regionaldiversity and therefore decrease national health expenditure.

Why England? Demand, growth and inequality during the Industrial Revolution

Why was England first? And why Europe? We present a probabilistic model that builds on big-push models by Murphy, Shleifer and Vishny (1989), combined with hierarchical preferences. The interaction of exogenous demographic factors (in particular the English low-pressure variant of the European marriage pattern)and redistributive institutions such as the old Poor Law combined to make an Industrial Revolution more likely. Essentially, industrialization is the result of having a critical mass of consumers that is rich enough to afford (potentially) mass-produced goods. Our model is then calibrated to match the main characteristics of the English economy in 1750 and the observed transition until 1850.This allows us to address explicitly one of the key features of the British IndustrialRevolution unearthed by economic historians over the last three decades the slowness of productivity and output change. In our calibration, we find that the probability of Britain industrializing is 5 times larger than France s. Contrary to the recent argument by Pomeranz, China in the 18th century had essentially no chance to industrialize at all. This difference is decomposed into a demographic and a policy component, with the former being far more important than the latter.

Forecasting and turning point predictions in a Bayesian panel VAR model

We provide methods for forecasting variables and predicting turning points in panel Bayesian VARs. We specify a flexible model which accounts for both interdependencies in the cross section and time variations in the parameters. Posterior distributions for the parameters are obtained for a particular type of diffuse, for Minnesota-type and for hierarchical priors. Formulas for multistep, multiunit point and average forecasts are provided. An application to the problem of forecasting the growth rate of output and of predicting turning points in the G-7 illustrates the approach. A comparison with alternative forecasting methods is also provided.