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

A comparison of the alr and ilr transformations for kernel density estimation of compositional data

In a seminal paper, Aitchison and Lauder (1985) introduced classical kernel densityestimation techniques in the context of compositional data analysis. Indeed, they gavetwo options for the choice of the kernel to be used in the kernel estimator. One ofthese kernels is based on the use the alr transformation on the simplex SD jointly withthe normal distribution on RD-1. However, these authors themselves recognized thatthis method has some deficiencies. A method for overcoming these dificulties based onrecent developments for compositional data analysis and multivariate kernel estimationtheory, combining the ilr transformation with the use of the normal density with a fullbandwidth matrix, was recently proposed in Martín-Fernández, Chacón and Mateu-Figueras (2006). Here we present an extensive simulation study that compares bothmethods in practice, thus exploring the finite-sample behaviour of both estimators

CODAMAT: a modern analogue techinque for compositional data

The quantitative estimation of Sea Surface Temperatures from fossils assemblages is afundamental issue in palaeoclimatic and paleooceanographic investigations. TheModern Analogue Technique, a widely adopted method based on direct comparison offossil assemblages with modern coretop samples, was revised with the aim ofconforming it to compositional data analysis. The new CODAMAT method wasdeveloped by adopting the Aitchison metric as distance measure. Modern coretopdatasets are characterised by a large amount of zeros. The zero replacement was carriedout by adopting a Bayesian approach to the zero replacement, based on a posteriorestimation of the parameter of the multinomial distribution. The number of modernanalogues from which reconstructing the SST was determined by means of a multipleapproach by considering the Proxies correlation matrix, Standardized Residual Sum ofSquares and Mean Squared Distance. This new CODAMAT method was applied to theplanktonic foraminiferal assemblages of a core recovered in the Tyrrhenian Sea.Kew words: Modern analogues, Aitchison distance, Proxies correlation matrix,Standardized Residual Sum of Squares

Pandemic 2009 A(H1N1) infection requiring hospitalization of elderly Spanish adults.

To describe the clinical presentation and prognosis of elderly adults hospitalized with pandemic 2009 A(H1N1) influenza infection and to compare these data with those of younger patients. DESIGN: Prospective, observational, multicenter study. SETTING: Thirteen hospitals in Spain. PARTICIPANTS: Adults admitted to the hospital with confirmed pandemic 2009 A(H1N1) influenza infection. MEASUREMENTS: Demographic, clinical, laboratory, radiological, and outcome variables. RESULTS: Between June 12 and November 10, 2009, 585 adults with confirmed 2009 A(H1N1) influenza were hospitalized, of whom 50 (8.5%) were aged 65 and older (median age 72, range 65-87). Older adults (≥ 65) were more likely to have associated comorbidities (88.0% vs 51.2%; P < .001), primarily chronic pulmonary diseases (46.0% vs 27.3%; P < .001). Lower respiratory tract symptoms and signs such as dyspnea (60.0% vs 45.6%) and wheezing (46.0% vs 27.8%; P = .007) were also more common in these elderly adults, although pulmonary infiltrates were present in just 14 (28.0%) of the older adults, compared with 221 (41.3%) of the younger adults (P = .06). Multilobar involvement was less frequent in elderly adults with pulmonary infiltrates than younger adults with pulmonary infiltrates (21.4% vs 60.0%; P = .05). Rhinorrhea (4.0% vs 21.9%; P = .003), myalgias (42.0% vs 59.1%; P = .01), and sore throat (14.0% vs 29.2%; P = .02) were more frequent in younger adults. Early antiviral therapy (<48 hours) was similar in the two groups (34.0% vs 37.9%; P = .58). Two older adults (4.0%) died during hospitalization, compared with 11 (2.1%) younger adults (P = .30). CONCLUSION: Elderly adults with 2009 A(H1N1) influenza had fewer viral-like upper respiratory symptoms than did younger adults. Pneumonia was more frequent in younger adults. No significant differences were observed in hospital mortality.

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 thetransfer 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 isessential in cases where origin is disputed. A robust method for discrimination of thismaterial through the use of elemental analysis and compositional data analysis is required.Initial studies have characterised the variability within a given nephrite source. This hasincluded investigation of both in situ outcrops and alluvial material. Methods for thediscrimination of two geographically close nephrite sources are being developed.Key Words: forensic, jade, nephrite, laser ablation, inductively coupled plasma massspectrometry, multivariate analysis, elemental analysis, compositional data analysis

A Factor analysis of hidrochemical composition of Llobregat river basin

Hydrogeological research usually includes some statistical studies devised to elucidate mean background state, characterise relationships among different hydrochemical parameters, and show the influence of human activities. These goals are achieved either by means of a statistical approach or by mixing modelsbetween end-members. Compositional data analysis has proved to be effective with the first approach, but there is no commonly accepted solution to the end-member problem in a compositional framework.We present here a possible solution based on factor analysis of compositions illustrated with a case study.We find two factors on the compositional bi-plot fitting two non-centered orthogonal axes to the most representative variables. Each one of these axes defines a subcomposition, grouping those variables thatlay nearest to it. With each subcomposition a log-contrast is computed and rewritten as an equilibrium equation. These two factors can be interpreted as the isometric log-ratio coordinates (ilr) of three hiddencomponents, that can be plotted in a ternary diagram. These hidden components might be interpreted as end-members.We have analysed 14 molarities in 31 sampling stations all along the Llobregat River and its tributaries, with a monthly measure during two years. We have obtained a bi-plot with a 57% of explained totalvariance, from which we have extracted two factors: factor G, reflecting geological background enhanced by potash mining; and factor A, essentially controlled by urban and/or farming wastewater. Graphicalrepresentation of these two factors allows us to identify three extreme samples, corresponding to pristine waters, potash mining influence and urban sewage influence. To confirm this, we have available analysisof diffused and widespread point sources identified in the area: springs, potash mining lixiviates, sewage, and fertilisers. Each one of these sources shows a clear link with one of the extreme samples, exceptfertilisers due to the heterogeneity of their composition.This approach is a useful tool to distinguish end-members, and characterise them, an issue generally difficult to solve. It is worth note that the end-member composition cannot be fully estimated but only characterised through log-ratio relationships among components. Moreover, the influence of each endmember in a given sample must be evaluated in relative terms of the other samples. These limitations areintrinsic to the relative nature of compositional data

The Life Course Determinants of Vulnerability in Late Careers

Late career is often seen as a more vulnerable life-stage in the labour market, in which workers may experience a deterioration in job quality. Using a life course perspective and longitudinal data, this article analyses the vulnerability associated with late career by focusing on four occupational dimensions: working-time, career continuity, retirement timing and income change. The research is carried out using data from Switzerland, a country where the age profile of the labour force is an increasing issue. The paper also adopts a cumulative disadvantage perspective to examine the impact of previous work and family life experiences on work life vulnerability at older age. Our data come from the Survey of Health, Ageing and Retirement in Europe (SHARELIFE). The paper uses cluster analysis, sequence analysis and ordered logistic regression. Results show that women with previous family responsibilities resulting in long-term unemployment or caring, often with health complications, are more likely to be vulnerable to deterioration in job quality in late career. This suggests that experiences in the last period of the working life may be just as gendered as earlier periods.

Computational anatomy for studying use-dependant brain plasticity.

In this article we provide a comprehensive literature review on the in vivo assessment of use-dependant brain structure changes in humans using magnetic resonance imaging (MRI) and computational anatomy. We highlight the recent findings in this field that allow the uncovering of the basic principles behind brain plasticity in light of the existing theoretical models at various scales of observation. Given the current lack of in-depth understanding of the neurobiological basis of brain structure changes we emphasize the necessity of a paradigm shift in the investigation and interpretation of use-dependent brain plasticity. Novel quantitative MRI acquisition techniques provide access to brain tissue microstructural properties (e.g., myelin, iron, and water content) in-vivo, thereby allowing unprecedented specific insights into the mechanisms underlying brain plasticity. These quantitative MRI techniques require novel methods for image processing and analysis of longitudinal data allowing for straightforward interpretation and causality inferences.

La implicació de les famílies immigrades en els processos educatius dels seus fills/es: estratègies davant l'èxit i el fracàs escolar

En el present projecte hem analitzat els determinants de les trajectòries educatives dels i les adolescents d'origen immigrant, centrant I'atenció en el paper de les seves famílies davant de I ‘èxit o fracàs escolar del seu fillla. Amb aquest objectiu, I'estudi combina tècniques quantitatives i qualitatives. Per una banda hem analitzat les dades longitudinals del Panel de Famílies i lnfancia, que ens permeten fer un seguiment de les trajectòries educatives i personals de 248 alumnes d'origen immigrant que al 2006 estudiaven I'ESO al llarg de la seva adolescència, i identificar els factors socials responsables de la seva diversificació. Els resultats indiquen que malgrat presentar actituds bastant favorables als estudis i I'assoliment educatiu, concentren diverses situacions de vulnerabilitat a la llar (dificultats socioeconòmiques, estructures familiars atípiques, i erosió de capital social), que incideixen negativament sobre els seus rendiments acadèmics. Per altra, hem realitzat 59 entrevistes semi-estructurades per a complementar i facilitar la interpretació dels resultats obtinguts a la recerca quantitativa i copsar les narratives dels propis protagonistes. Aquestes entrevistes s'han realitzat a: una submostra de les famílies d'aquests alumnes, seleccionades en funció de perfils d’èxit o fracàs educatiu de la trajectòria del menor (46), una submostra d'estudiants resilients (a), i una sèrie d'agents educatius i socials, que inclou membres d'equips directius de centres escolars, AMPA i entitats dedicades a I'atenció a la infància i les famílies (5). El projecte que presentem té una clara vocació de servei públic. L'objectiu és incrementar el coneixement de factors "extraescolars" que poden condicionar I ‘èxit escolar dels estudiants d'origen immigrant. Aquest coneixement constitueix la base per al disseny i orientació de programes d'acompanyament a les famílies dels infants en situació de risc. La nostra voluntat (que reflecteix el principal objectiu de I'lnstitut d’infància i Món Urbà, instituci6 que impulsa el projecte) és contribuir a la transferència de coneixement que pugui ser d'utilitat pels agents que treballen directament sobre les qüestions que estudiem.

Italy's Path to Very Low Fertility : The Adequacy of Economic and Second Demographic Transition Theories

The deep drop of the fertility rate in Italy to among the lowest in the world challenges contemporary theories of childbearing and family building. Among high-income countries, Italy was presumed to have characteristics of family values and female labor force participation that would favor higher fertility than its European neighbors to the north. We test competing economic and cultural explanations, drawing on new nationally representative, longitudinal data to examine first union, first birth, and second birth. Our event history analysis finds some support for economic determinants of family formation and fertility, but the clear importance of regional differences and of secularization suggests that such an explanation is at best incomplete and that cultural and ideational factors must be considered.

Late bloomer or language impaired? Screening the speech of a Spanish-English bilingual toddler for early signs of language impairment

This case study presents corpus data gathered from a Spanish-English bilingual child with expressive language delay. Longitudinal data on the child’s linguistic development was collected from the onset of productive speech at age 1;1 until age 4 over the course of 28 video-taped sessions with the child’s principal caregivers. A literature review focused on the relationship between language delay and persisting disorders—including a discussion of the frequent difficulty in distinguishing between the two at early stages of bilingual development—is followed by an analysis of the child’s productive development in 2 distinct phases. An attempt is made to assess the child’s speech at age 4 for preliminary signs of SLI and to consider techniques for identifying ‘at risk’ bilingual children (that is, those with productive language delay, poor oral fluency, and family history of language problems) based on samples of recorded and transcribed speech.

Marijuana use, diet, body mass index, and cardiovascular risk factors (from the CARDIA study).

Marijuana use has been associated with increased appetite, high caloric diet, acute increase in blood pressure, and decreases in high-density lipoprotein cholesterol and triglycerides. Marijuana is the most commonly used illicit drug in the United States, but its long-term effects on body mass index (BMI) and cardiovascular risk factors are unknown. Using 15 years of longitudinal data from 3,617 black and white young adults participating in the Coronary Artery Risk Development in Young Adults (CARDIA) study, we assessed whether marijuana use was associated with caloric intake, BMI, and cardiovascular risk factors. Of the 3,617 participants, 1,365 (38%) reported ever using marijuana. Marijuana use was associated with male gender, tobacco smoking, and other illicit drug use. More extensive marijuana use was associated with a higher caloric intake (2,746 kcal/day in never users to 3,365 kcal/day in those who used marijuana for > or = 1,800 days over 15 years) and alcohol intake (3.6 to 10.8 drinks/week), systolic blood pressure (112.7 to 116.5 mm Hg), and triglyceride levels (84 to 100 mg/dl or 0.95 to 1.13 mmol/L, all p values for trend < 0.001), but not with higher BMI and lipid and glucose levels. In multivariate analysis, the associations between marijuana use and systolic blood pressure and triglycerides disappeared, having been mainly confounded by greater alcohol use in marijuana users. In conclusion, although marijuana use was not independently associated with cardiovascular risk factors, it was associated with other unhealthy behaviors, such as high caloric diet, tobacco smoking, and other illicit drug use, which all have long-term detrimental effects on health.

Rules of origin : from analysis to reform

General Introduction This thesis can be divided into two main parts :the first one, corresponding to the first three chapters, studies Rules of Origin (RoOs) in Preferential Trade Agreements (PTAs); the second part -the fourth chapter- is concerned with Anti-Dumping (AD) measures. Despite wide-ranging preferential access granted to developing countries by industrial ones under North-South Trade Agreements -whether reciprocal, like the Europe Agreements (EAs) or NAFTA, or not, such as the GSP, AGOA, or EBA-, it has been claimed that the benefits from improved market access keep falling short of the full potential benefits. RoOs are largely regarded as a primary cause of the under-utilization of improved market access of PTAs. RoOs are the rules that determine the eligibility of goods to preferential treatment. Their economic justification is to prevent trade deflection, i.e. to prevent non-preferred exporters from using the tariff preferences. However, they are complex, cost raising and cumbersome, and can be manipulated by organised special interest groups. As a result, RoOs can restrain trade beyond what it is needed to prevent trade deflection and hence restrict market access in a statistically significant and quantitatively large proportion. Part l In order to further our understanding of the effects of RoOs in PTAs, the first chapter, written with Pr. Olivier Cadot, Celine Carrère and Pr. Jaime de Melo, describes and evaluates the RoOs governing EU and US PTAs. It draws on utilization-rate data for Mexican exports to the US in 2001 and on similar data for ACP exports to the EU in 2002. The paper makes two contributions. First, we construct an R-index of restrictiveness of RoOs along the lines first proposed by Estevadeordal (2000) for NAFTA, modifying it and extending it for the EU's single-list (SL). This synthetic R-index is then used to compare Roos under NAFTA and PANEURO. The two main findings of the chapter are as follows. First, it shows, in the case of PANEURO, that the R-index is useful to summarize how countries are differently affected by the same set of RoOs because of their different export baskets to the EU. Second, it is shown that the Rindex is a relatively reliable statistic in the sense that, subject to caveats, after controlling for the extent of tariff preference at the tariff-line level, it accounts for differences in utilization rates at the tariff line level. Finally, together with utilization rates, the index can be used to estimate total compliance costs of RoOs. The second chapter proposes a reform of preferential Roos with the aim of making them more transparent and less discriminatory. Such a reform would make preferential blocs more "cross-compatible" and would therefore facilitate cumulation. It would also contribute to move regionalism toward more openness and hence to make it more compatible with the multilateral trading system. It focuses on NAFTA, one of the most restrictive FTAs (see Estevadeordal and Suominen 2006), and proposes a way forward that is close in spirit to what the EU Commission is considering for the PANEURO system. In a nutshell, the idea is to replace the current array of RoOs by a single instrument- Maximum Foreign Content (MFC). An MFC is a conceptually clear and transparent instrument, like a tariff. Therefore changing all instruments into an MFC would bring improved transparency pretty much like the "tariffication" of NTBs. The methodology for this exercise is as follows: In step 1, I estimate the relationship between utilization rates, tariff preferences and RoOs. In step 2, I retrieve the estimates and invert the relationship to get a simulated MFC that gives, line by line, the same utilization rate as the old array of Roos. In step 3, I calculate the trade-weighted average of the simulated MFC across all lines to get an overall equivalent of the current system and explore the possibility of setting this unique instrument at a uniform rate across lines. This would have two advantages. First, like a uniform tariff, a uniform MFC would make it difficult for lobbies to manipulate the instrument at the margin. This argument is standard in the political-economy literature and has been used time and again in support of reductions in the variance of tariffs (together with standard welfare considerations). Second, uniformity across lines is the only way to eliminate the indirect source of discrimination alluded to earlier. Only if two countries face uniform RoOs and tariff preference will they face uniform incentives irrespective of their initial export structure. The result of this exercise is striking: the average simulated MFC is 25% of good value, a very low (i.e. restrictive) level, confirming Estevadeordal and Suominen's critical assessment of NAFTA's RoOs. Adopting a uniform MFC would imply a relaxation from the benchmark level for sectors like chemicals or textiles & apparel, and a stiffening for wood products, papers and base metals. Overall, however, the changes are not drastic, suggesting perhaps only moderate resistance to change from special interests. The third chapter of the thesis considers whether Europe Agreements of the EU, with the current sets of RoOs, could be the potential model for future EU-centered PTAs. First, I have studied and coded at the six-digit level of the Harmonised System (HS) .both the old RoOs -used before 1997- and the "Single list" Roos -used since 1997. Second, using a Constant Elasticity Transformation function where CEEC exporters smoothly mix sales between the EU and the rest of the world by comparing producer prices on each market, I have estimated the trade effects of the EU RoOs. The estimates suggest that much of the market access conferred by the EAs -outside sensitive sectors- was undone by the cost-raising effects of RoOs. The chapter also contains an analysis of the evolution of the CEECs' trade with the EU from post-communism to accession. Part II The last chapter of the thesis is concerned with anti-dumping, another trade-policy instrument having the effect of reducing market access. In 1995, the Uruguay Round introduced in the Anti-Dumping Agreement (ADA) a mandatory "sunset-review" clause (Article 11.3 ADA) under which anti-dumping measures should be reviewed no later than five years from their imposition and terminated unless there was a serious risk of resumption of injurious dumping. The last chapter, written with Pr. Olivier Cadot and Pr. Jaime de Melo, uses a new database on Anti-Dumping (AD) measures worldwide to assess whether the sunset-review agreement had any effect. The question we address is whether the WTO Agreement succeeded in imposing the discipline of a five-year cycle on AD measures and, ultimately, in curbing their length. Two methods are used; count data analysis and survival analysis. First, using Poisson and Negative Binomial regressions, the count of AD measures' revocations is regressed on (inter alia) the count of "initiations" lagged five years. The analysis yields a coefficient on measures' initiations lagged five years that is larger and more precisely estimated after the agreement than before, suggesting some effect. However the coefficient estimate is nowhere near the value that would give a one-for-one relationship between initiations and revocations after five years. We also find that (i) if the agreement affected EU AD practices, the effect went the wrong way, the five-year cycle being quantitatively weaker after the agreement than before; (ii) the agreement had no visible effect on the United States except for aone-time peak in 2000, suggesting a mopping-up of old cases. Second, the survival analysis of AD measures around the world suggests a shortening of their expected lifetime after the agreement, and this shortening effect (a downward shift in the survival function postagreement) was larger and more significant for measures targeted at WTO members than for those targeted at non-members (for which WTO disciplines do not bind), suggesting that compliance was de jure. A difference-in-differences Cox regression confirms this diagnosis: controlling for the countries imposing the measures, for the investigated countries and for the products' sector, we find a larger increase in the hazard rate of AD measures covered by the Agreement than for other measures.

Overview of hospitalizations by ambulatory care sensitive conditions in the municipality of Cotia, Brazil

Objective To describe the profile of Hospitalizations by Amulatory Care Sensitive Conditions (HACSC), in the Municipality of Cotia, from 2008 to 2012. Method ecological, exploratory, longitudinal study with a quantitative approach. Data on HACSC, by age group and sex, were obtained from the Department of the Unified Health System. For data analysis descriptive statistics were used. Results During the period, there were 46,676 admissions, excluding deliveries, 7,753 (16.61%) by HACSC. The main causes were cerebrovascular diseases, 16.96%, heart failure, 15.50%, hypertension, 10.80% and infection of the kidney and urinary tract, 10.51%. Regarding gender, HACSC occurred predominantly in males. There was a greater number of HACSC at extreme age ranges, especially in the elderly. Conclusion Chronic diseases predominate among the leading causes of HACSC and there was no significant difference between sex.





Compositional Data in Biomedical Research

Modern methods of compositional data analysis are not well known in biomedical research.Moreover, there appear to be few mathematical and statistical researchersworking on compositional biomedical problems. Like the earth and environmental sciences,biomedicine has many problems in which the relevant scienti c information isencoded in the relative abundance of key species or categories. I introduce three problemsin cancer research in which analysis of compositions plays an important role. Theproblems involve 1) the classi cation of serum proteomic pro les for early detection oflung cancer, 2) inference of the relative amounts of di erent tissue types in a diagnostictumor biopsy, and 3) the subcellular localization of the BRCA1 protein, and it'srole in breast cancer patient prognosis. For each of these problems I outline a partialsolution. However, none of these problems is \solved". I attempt to identify areas inwhich additional statistical development is needed with the hope of encouraging morecompositional data analysts to become involved in biomedical research