An R Library for Compositional Data Analysis in Archaeometry

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

A compositional data analysis package for R providing multiple approaches

���compositions��� is a new R-package for the analysis of compositional and positive data.It contains four classes corresponding to the four different types of compositional andpositive geometry (including the Aitchison geometry). It provides means for computation,plotting and high-level multivariate statistical analysis in all four geometries.These geometries are treated in an fully analogous way, based on the principle of workingin coordinates, and the object-oriented programming paradigm of R. In this way,called functions automatically select the most appropriate type of analysis as a functionof the geometry. The graphical capabilities include ternary diagrams and tetrahedrons,various compositional plots (boxplots, barplots, piecharts) and extensive graphical toolsfor principal components. Afterwards, ortion and proportion lines, straight lines andellipses in all geometries can be added to plots. The package is accompanied by ahands-on-introduction, documentation for every function, demos of the graphical capabilitiesand plenty of usage examples. It allows direct and parallel computation inall four vector spaces and provides the beginner with a copy-and-paste style of dataanalysis, while letting advanced users keep the functionality and customizability theydemand of R, as well as all necessary tools to add own analysis routines. A completeexample is included in the appendix

Exploration of geological variability and possible processes through the use of compositional data analysis: an example using scottish metamorphosed

Developments in the statistical analysis of compositional data over the last twodecades have made possible a much deeper exploration of the nature of variability,and the possible processes associated with compositional data sets from manydisciplines. In this paper we concentrate on geochemical data sets. First we explainhow hypotheses of compositional variability may be formulated within the naturalsample space, the unit simplex, including useful hypotheses of subcompositionaldiscrimination and specific perturbational change. Then we develop through standardmethodology, such as generalised likelihood ratio tests, statistical tools to allow thesystematic investigation of a complete lattice of such hypotheses. Some of these tests are simple adaptations of existing multivariate tests but others require specialconstruction. We comment on the use of graphical methods in compositional dataanalysis and on the ordination of specimens. The recent development of the conceptof compositional processes is then explained together with the necessary tools for astaying- in-the-simplex approach, namely compositional singular value decompositions. All these statistical techniques are illustrated for a substantial compositional data set, consisting of 209 major-oxide and rare-element compositions of metamorphosed limestones from the Northeast and Central Highlands of Scotland.Finally we point out a number of unresolved problems in the statistical analysis ofcompositional processes

Determinants and burden of chronic kidney disease in the population-based CoLaus study: a cross-sectional analysis.

BACKGROUND: Chronic kidney disease (CKD) represents an increasing health burden. We present the population-based prevalence of CKD and compare the CKD Epidemiology collaboration (CKD-EPI) and modification of diet in renal disease (MDRD) equations to estimate the glomerular filtration rate, using the revised CKD classification with three albuminuria classes. We also explore factors associated with CKD. METHODS: The Swiss population-based, cross-sectional CoLaus study conducted in Lausanne (2003-2006) included 2810 men and 3111 women aged 35-75. CKD prevalence was assessed using CKD-EPI and MDRD equations and albuminuria estimated by the albumin-to-creatinine ratio in spot morning urine. Multivariate logistic regression was used to analyse determinants of CKD. RESULTS: Prevalence [95% confidence interval (CI)] of all stages CKD was 10.0% (9.2-10.8%) with CKD-EPI and 13.8% (12.9-14.6%) with MDRD. Using the revised CKD classification, the prevalence of low-, medium-, high- and very high-risk groups was 90.0, 8.46, 1.18 and 0.35% with CKD-EPI, respectively. With MDRD, the corresponding values were 86.24, 11.86, 1.55 and 0.35%. Using the revised classification, CKD-EPI systematically reclassified people in a lower risk category than MDRD. Age and obesity were more strongly associated with CKD in men [odds ratio (95% CI): 2.23(1.95; 2.56) per 10 years and 3.05(2.08;4.47), respectively] than in women [1.46 (1.29; 1.65) and 1.78 (1.30;2.44), respectively]. Hypertension, type 2 diabetes, serum homocysteine and uric acid were positively independently associated with CKD in men and women. CONCLUSIONS: One in 10 adults suffers from CKD in the population of Lausanne. CKD-EPI systematically reclassifies people in a lower CKD risk category than MDRD. Serum homocysteine and uric acid levels are associated with CKD independently of classical risk factors such as age, hypertension and diabetes.

Improving prediction of recanalization in acute large-vessel occlusive stroke.

BACKGROUND: Recanalization in acute ischemic stroke with large-vessel occlusion is a potent indicator of good clinical outcome. OBJECTIVE: To identify easily available clinical and radiologic variables predicting recanalization at various occlusion sites. METHODS: All consecutive, acute stroke patients from the Acute STroke Registry and Analysis of Lausanne (2003-2011) who had a large-vessel occlusion on computed tomographic angiography (CTA) (<��12��h) were included. Recanalization status was assessed at 24��h (range: 12-48��h) with CTA, magnetic resonance angiography, or ultrasonography. Complete and partial recanalization (corresponding to the modified Treatment in Cerebral Ischemia scale��2-3) were grouped together. Patients were categorized according to occlusion site and treatment modality. RESULTS: Among 439 patients, 51% (224) showed complete or partial recanalization. In multivariate analysis, recanalization of any occlusion site was most strongly associated with endovascular treatment, including bridging therapy (odds ratio [OR]��7.1, 95% confidence interval [CI]��2.2-23.2), and less so with intravenous thrombolysis (OR��1.6, 95%��CI��1.0-2.6) and recanalization treatments performed beyond guidelines (OR��2.6, 95%��CI��1.2-5.7). Clot location (large vs. intermediate) and tandem pathology (the combination of intracranial occlusion and symptomatic extracranial stenosis) were other variables discriminating between recanalizers and non-recanalizers. For patients with intracranial occlusions, the variables significantly associated with recanalization after 24��h were: baseline National Institutes of Health Stroke Scale (NIHSS) (OR��1.04, 95%��CI��1.02-1.1), Alberta Stroke Program Early CT Score (ASPECTS) on initial computed tomography (OR��1.2, 95%��CI��1.1-1.3), and an altered level of consciousness (OR��0.2, 95%��CI��0.1-0.5). CONCLUSIONS: Acute endovascular treatment is the single most important factor promoting recanalization in acute ischemic stroke. The presence of extracranial vessel stenosis or occlusion decreases recanalization rates. In patients with intracranial occlusions, higher NIHSS score and ASPECTS and normal vigilance facilitate recanalization. Clinical use of these predictors could influence recanalization strategies in individual patients.

Methods for the ex vivo characterization of human CD8+ T subsets based on gene expression and replicative history analysis.

The generation of an antigen-specific T-lymphocyte response is a complex multi-step process. Upon T-cell receptor-mediated recognition of antigen presented by activated dendritic cells, naive T-lymphocytes enter a program of proliferation and differentiation, during the course of which they acquire effector functions and may ultimately become memory T-cells. A major goal of modern immunology is to precisely identify and characterize effector and memory T-cell subpopulations that may be most efficient in disease protection. Sensitive methods are required to address these questions in exceedingly low numbers of antigen-specific lymphocytes recovered from clinical samples, and not manipulated in vitro. We have developed new techniques to dissect immune responses against viral or tumor antigens. These allow the isolation of various subsets of antigen-specific T-cells (with major histocompatibility complex [MHC]-peptide multimers and five-color FACS sorting) and the monitoring of gene expression in individual cells (by five-cell reverse transcription-polymerase chain reaction [RT-PCR]). We can also follow their proliferative life history by flow-fluorescence in situ hybridization (FISH) analysis of average telomere length. Recently, using these tools, we have identified subpopulations of CD8+ T-lymphocytes with distinct proliferative history and partial effector-like properties. Our data suggest that these subsets descend from recently activated T-cells and are committed to become differentiated effector T-lymphocytes.

Meat consumption and mortality - results from the European Prospective Investigation into Cancer and Nutrition.

BACKGROUND Recently, some US cohorts have shown a moderate association between red and processed meat consumption and mortality supporting the results of previous studies among vegetarians. The aim of this study was to examine the association of red meat, processed meat, and poultry consumption with the risk of early death in the European Prospective Investigation into Cancer and Nutrition (EPIC). METHODS Included in the analysis were 448,568 men and women without prevalent cancer, stroke, or myocardial infarction, and with complete information on diet, smoking, physical activity and body mass index, who were between 35 and 69 years old at baseline. Cox proportional hazards regression was used to examine the association of meat consumption with all-cause and cause-specific mortality. RESULTS As of June 2009, 26,344 deaths were observed. After multivariate adjustment, a high consumption of red meat was related to higher all-cause mortality (hazard ratio (HR) = 1.14, 95% confidence interval (CI) 1.01 to 1.28, 160+ versus 10 to 19.9 g/day), and the association was stronger for processed meat (HR = 1.44, 95% CI 1.24 to 1.66, 160+ versus 10 to 19.9 g/day). After correction for measurement error, higher all-cause mortality remained significant only for processed meat (HR = 1.18, 95% CI 1.11 to 1.25, per 50 g/d). We estimated that 3.3% (95% CI 1.5% to 5.0%) of deaths could be prevented if all participants had a processed meat consumption of less than 20 g/day. Significant associations with processed meat intake were observed for cardiovascular diseases, cancer, and 'other causes of death'. The consumption of poultry was not related to all-cause mortality. CONCLUSIONS The results of our analysis support a moderate positive association between processed meat consumption and mortality, in particular due to cardiovascular diseases, but also to cancer.

A survey strategy for human respiratory syncytial virus detection among haematopoietic stem cell transplant patients: epidemiological and methodological analysis

Human respiratory syncytial virus (HRSV) causes severe infections among children and immunocompromised patients. We compared HRSV infections among Haematopoietic Stem Cell Transplant program (HSCT) patients and children using direct immunofluorescence (DFA), point-of-care RSV Bio Easy�� and a polymerase chain reaction (PCR) assay. Overall, 102 samples from HSCT patients and 128 from children obtained positivity rate of 18.6% and 14.1% respectively. PCR sensitivity was highest mainly on samples collected after five days of symptoms onset. A combination of both DFA and reverse transcriptase-PCR methods for HSCT high-risk patients is the best diagnostic flow for HRSV diagnosis among these patients.

Transcriptome analysis of the mobile genome ICEclc in Pseudomonas knackmussii B13.

Background: Integrative and conjugative elements (ICE) form a diverse group of DNA elements that are integrated in the chromosome of the bacterial host, but can occasionally excise and horizontally transfer to a new host cell. ICE come in different families, typically with a conserved core for functions controlling the element's behavior and a variable region providing auxiliary functions to the host. The ICEclc element of Pseudomonas knackmussii strain B13 is representative for a large family of chromosomal islands detected by genome sequencing approaches. It provides the host with the capacity to degrade chloroaromatics and 2-aminophenol. Results: Here we study the transcriptional organization of the ICEclc core region. By northern hybridizations, reverse-transcriptase polymerase chain reaction (RT-PCR) and Rapid Amplification of cDNA Ends (5'-RACE) fifteen transcripts were mapped in the core region. The occurrence and location of those transcripts were further confirmed by hybridizing labeled cDNA to a semi-tiling micro-array probing both strands of the ICEclc core region. Dot blot and semi-tiling array hybridizations demonstrated most of the core transcripts to be upregulated during stationary phase on 3-chlorobenzoate, but not on succinate or glucose. Conclusions: The transcription analysis of the ICEclc core region provides detailed insights in the mode of regulatory organization and will help to further understand the complex mode of behavior of this class of mobile elements. We conclude that ICEclc core transcription is concerted at a global level, more reminiscent of a phage program than of plasmid conjugation.

Effectiveness of holistic neuropsychological rehabilitation for Spanish population with acquired brain injury measured using Rasch analysis.

INTRODUCTION The Rasch model is increasingly used in the field of rehabilitation because it improves the accuracy of measurements of patient status and their changes after therapy. OBJECTIVE To determine the long-term effectiveness of a holistic neuropsychological rehabilitation program for Spanish outpatients with acquired brain injury (ABI) using Rasch analysis. METHODS Eighteen patients (ten with long evolution - patients who started the program > 6 months after ABI- and eight with short evolution) and their relatives attended the program for 6 months. Patients' and relatives' answers to the European Brain Injury Questionnaire and the Frontal Systems Behavior Scale at 3 time points (pre-intervention. post-intervention and 12 month follow-up) were transformed into linear measures called logits. RESULTS The linear measures revealed significant improvements with large effects at the follow-up assessment on cognitive and executive functioning, social and emotional self-regulation, apathy and mood. At follow-up, the short evolution group achieved greater improvements in mood and cognitive functioning than the long evolution patients. CONCLUSIONS The program showed long-term effectiveness for most of the variables, and it was more effective for mood and cognitive functioning when patients were treated early. Relatives played a key role in the effectiveness of the rehabilitation program.

The TT genotype of the STAT4 rs7574865 polymorphism is associated with high disease activity and disability in patients with early arthritis

BACKGROUND The number of copies of the HLA-DRB1 shared epitope, and the minor alleles of the STAT4 rs7574865 and the PTPN22 rs2476601 polymorphisms have all been linked with an increased risk of developing rheumatoid arthritis. In the present study, we investigated the effects of these genetic variants on disease activity and disability in patients with early arthritis. METHODOLOGY AND RESULTS We studied 640 patients with early arthritis (76% women; median age, 52 years), recording disease-related variables every 6 months during a 2-year follow-up. HLA-DRB1 alleles were determined by PCR-SSO, while rs7574865 and rs2476601 were genotyped with the Taqman 5' allelic discrimination assay. Multivariate analysis was performed using generalized estimating equations for repeated measures. After adjusting for confounding variables such as gender, age and ACPA, the TT genotype of rs7574865 in STAT4 was associated with increased disease activity (DAS28) as compared with the GG genotype (�� coefficient [95% confidence interval] = 0.42 [0.01-0.83], p = 0.044). Conversely, the presence of the T allele of rs2476601 in PTPN22 was associated with diminished disease activity during follow-up in a dose-dependent manner (CT genotype = -0.27 [-0.56- -0.01], p = 0.042; TT genotype = -0.68 [-1.64- -0.27], p = 0.162). After adjustment for gender, age and disease activity, homozygosity for the T allele of rs7574865 in STAT4 was associated with greater disability as compared with the GG genotype. CONCLUSIONS Our data suggest that patients with early arthritis who are homozygous for the T allele of rs7574865 in STAT4 may develop a more severe form of the disease with increased disease activity and disability.

Implementation of a nutrition education program in a handball team; consequences on nutritional status.

OBJECTIVE: To evaluate nutritional status and dietary habits after implementation of a nutritional education program in professional handball players. RESEARCH METHODS AND PROCEDURES: Longitudinal study of 14 handball players evaluated with 72-h recall, a questionnaire on food consumption and anthropometric measures during 4 months. The intervention consisted of a nutrition education program. RESULTS: Energy intake was consistently below the recommended allowances. Macronutrient intakes as a percentage of total energy intake were below the recommended allowances for carbohydrates, and above recommended allowances for fats. Nutritional education was followed by a significant increase (p < 0.01) in total energy and macronutrient intakes, with no significant changes in macronutrient or micronutrient intakes after adjustment for energy intake. DISCUSSION: The imbalance in nutrient intake in handball players suggests that detailed re-analysis is needed to determine specific recommendations for this population. Nutritional education with continuous follow-up to monitor athletes' dietary habits may lead them to adopt appropriate nutritional habits to optimize dietary intakes. The lack of specific recommendations for micronutrient intakes in athletes leads to confusion regarding appropriate intakes; biochemical tests that yield normal values (albeit approaching cut-off values for deficiency) may disguise deficient status for some nutrients when strenuous exercise is involved. CONCLUSION: In-depth studies with nutrition education programs that include long-term follow-up are advisable to avoid deficiencies that can lead to irreversible damage in competitive athletes.

Obesity and survival in operable breast cancer patients treated with adjuvant anthracyclines and taxanes according to pathological subtypes: a pooled analysis.

INTRODUCTION Obesity is an unfavorable prognostic factor in breast cancer (BC) patients regardless of menopausal status and treatment received. However, the association between obesity and survival outcome by pathological subtype requires further clarification. METHODS We performed a retrospective analysis including 5,683 operable BC patients enrolled in four randomized clinical trials (GEICAM/9906, GEICAM/9805, GEICAM/2003-02, and BCIRG 001) evaluating anthracyclines and taxanes as adjuvant treatments. Our primary aim was to assess the prognostic effect of body mass index (BMI) on disease recurrence, breast cancer mortality (BCM), and overall mortality (OM). A secondary aim was to detect differences of such prognostic effects by subtype. RESULTS Multivariate survival analyses adjusting for age, tumor size, nodal status, menopausal status, surgery type, histological grade, hormone receptor status, human epidermal growth factor receptor 2 (HER2) status, chemotherapy regimen, and under-treatment showed that obese patients (BMI 30.0 to 34.9) had similar prognoses to that of patients with a BMI���<���25 (reference group) in terms of recurrence (Hazard Ratio [HR]���=���1.08, 95% Confidence Interval [CI]���=���0.90 to 1.30), BCM (HR���=���1.02, 0.81 to 1.29), and OM (HR���=���0.97, 0.78 to 1.19). Patients with severe obesity (BMI���������35) had a significantly increased risk of recurrence (HR���=���1.26, 1.00 to 1.59, P���=���0.048), BCM (HR���=���1.32, 1.00 to 1.74, P���=���0.050), and OM (HR���=���1.35, 1.06 to 1.71, P���=���0.016) compared to our reference group. The prognostic effect of severe obesity did not vary by subtype. CONCLUSIONS Severely obese patients treated with anthracyclines and taxanes present a worse prognosis regarding recurrence, BCM, and OM than patients with BMI���<���25. The magnitude of the harmful effect of BMI on survival-related outcomes was similar across subtypes.

Economic and environmental effects of the CO2 taxation: an input-output analysis for Spain

CO2 emissions induced by human activities are the major cause of climate change; hence, strong environmental policy that limits the growing dependence on fossil fuel is indispensable. Tradable permits and environmental taxes are the usual tools used in CO2 reduction strategies. Such economic tools provide incentives to polluting industries to reduce their emissions through market signals. The aim of this work is to investigate the direct and indirect effects of an environmental tax on Spanish products and services. We apply an environmentally extended input-output (EIO) model to identify CO2 emission intensities of products and services and, accordingly, we estimate the tax proportional to these intensities. The short-term price effects are analyzed using an input-output price model. The effect of tax introduction on consumption prices and its influence on consumers��� welfare are determined. We also quantify the environmental impacts of such taxation in terms of the reduction in CO2 emissions. The results, based on the Spanish economy for the year 2007, show that sectors with relatively poor environmental profile are subjected to high environmental tax rates. And consequently, applying a CO2 tax on these sectors, increases production prices and induces a slight increase in consumer price index and a decrease in private welfare. The revenue from the tax could be used to counter balance the negative effects on social welfare and also to stimulate the increase of renewable energy shares in the most impacting sectors. Finally, our analysis highlights that the environmental and economic goals cannot be met at the same time with the environmental taxation and this shows the necessity of finding other (complementary or alternative) measures to ensure both the economic and ecological efficiencies. Keywords: CO2 emissions; environmental tax; input-output model, effects of environmental taxation.

Positive deviance study to inform a Chagas disease control program in southern Ecuador

Chagas disease is caused by Trypanosoma cruzi, which is mainly transmitted by the faeces of triatomine insects that find favourable environments in poorly constructed houses. Previous studies have documented persistent triatomine infestation in houses in the province of Loja in southern Ecuador despite repeated insecticide and educational interventions. We aim to develop a sustainable strategy for the interruption of Chagas disease transmission by promoting living environments that are designed to prevent colonisation of rural houses by triatomines. This study used positive deviance to inform the design of an anti-triatomine prototype house by identifying knowledge, attitudes and practices used by families that have remained triatomine-free (2010-2012). Positive deviants reported practices that included maintenance of structural elements of the house, fumigation of dwellings and animal shelters, sweeping with "insect repellent" plants and relocation of domestic animals away from the house, among others. Participants favoured construction materials that do not drastically differ from those currently used (adobe walls and tile roofs). They also expressed their belief in a clear connection between a clean house and health. The family's economic dynamics affect space use and must be considered in the prototype's design. Overall, the results indicate a positive climate for the introduction of housing improvements as a protective measure against Chagas disease in this region.