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

Y-90 Ibritumomab Tiuxetan (Zevalin (R)) Consolidation of First Remission In Advanced Stage Follicular Non Hodgkin's Lymphoma Updated Results After a Median Follow up of 662 Months From the International, Randomized, Phase III First Line Indolent Trial (FIT) In 414 Patients

The FIT trial was conducted to evaluate the safety and efficacy of 90Y-ibritumomab tiuxetan (0.4 mCi/kg; maximum dose 32 mCi) when used as consolidation of first complete or partial remission in patients with previously untreated, advanced-stage follicular lymphoma (FL). Patients were randomly assigned to either 90Y-ibritumomab treatment (n = 207) or observation (n = 202) within 3 months (mo) of completing initial induction therapy (chemotherapy only: 86%; rituximab in combination with chemotherapy: 14%). Response status prior to randomization did not differ between the groups: 52% complete response (CR)/CR unconfirmed (CRu) to induction therapy and 48% partial response (PR) in the 90Y-ibritumomab arm vs 53% CR/CRu and 44% PR in the control arm. The primary endpoint was progression-free survival (PFS) of the intent-to-treat (ITT) population. Results from the first extended follow-up after a median of 3.5 years revealed a significant improvement in PFS from the time of randomization with 90Y-ibritumomab consolidation compared with control (36.5 vs 13.3 mo, respectively; P < 0.0001; Morschhauser et al. JCO. 2008; 26:5156-5164). Here we report a median follow-up of 66.2 mo (5.5 years). Five-year PFS was 47% in the 90Y-ibritumomab group and 29% in the control group (hazard ratio (HR) = 0.51, 95% CI 0.39-0.65; P < 0.0001). Median PFS in the 90Y-ibritumomab group was 49 mo vs 14 mo in the control group. In patients achieving a CR/CRu after induction, 5-year PFS was 57% in the 90Y-ibritumomab group, and the median had not yet been reached at 92 months, compared with a 43% 5-year PFS in the control group and a median of 31 mo (HR = 0.61, 95% CI 0.42-0.89). For patients in PR after induction, the 5-year PFS was 38% in the 90Y-ibritumomab group with a median PFS of 30 mo vs 14% in the control group with a median PFS of 6 mo (HR = 0.38, 95% CI 0.27-0.53). Patients who had received rituximab as part of induction treatment had a 5-year PFS of 64% in the 90Y-ibritumomab group and 48% in the control group (HR = 0.66, 95% CI 0.30-1.47). For all patients, time to next treatment (as calculated from the date of randomization) differed significantly between both groups; median not reached at 99 mo in the 90Y-ibritumomab group vs 35 mo in the control group (P < 0.0001). The majority of patients received rituximab-containing regimens when treated after progression (63/82 [77%] in the 90Y-ibritumomab group and 102/122 [84%] in the control group). Overall response rate to second-line treatment was 79% in the 90Y-ibritumomab group (57% CR/CRu and 22% PR) vs 78% in the control arm (59% CR/CRu, 19% PR). Five-year overall survival was not significantly different between the groups; 93% and 89% in the 90Y-ibritumomab and control groups, respectively (P = 0.561). To date, 40 patients have died; 18 in the 90Y-ibritumomab group and 22 in the control group. Secondary malignancies were diagnosed in 16 patients in the 90Y-ibritumomab arm vs 9 patients in the control arm (P = 0.19). There were 6 (3%) cases of myelodysplastic syndrome (MDS)/acute myelogenous leukemia (AML) in the 90Y-ibritumomab arm vs 1 MDS in the control arm (P = 0.063). In conclusion, this extended follow-up of the FIT trial confirms the benefit of 90Y-ibritumomab consolidation with a nearly 3 year advantage in median PFS. A significant 5-year PFS improvement was confirmed for patients with a CR/CRu or a PR after induction. Effective rescue treatment with rituximab-containing regimens may explain the observed no difference in overall survival between both patient groups who were - for the greater part - rituximab-naïve.

Compositional Data Analysis with R

R from http://www.r-project.org/ is ‘GNU S’ – a language and environment for statistical computingand graphics. The environment in which many classical and modern statistical techniques havebeen implemented, but many are supplied as packages. There are 8 standard packages and many moreare available through the cran family of Internet sites http://cran.r-project.org .We started to develop a library of functions in R to support the analysis of mixtures and our goal isa MixeR package for compositional data analysis that provides support foroperations on compositions: perturbation and power multiplication, subcomposition with or withoutresiduals, centering of the data, computing Aitchison’s, Euclidean, Bhattacharyya distances,compositional Kullback-Leibler divergence etc.graphical presentation of compositions in ternary diagrams and tetrahedrons with additional features:barycenter, geometric mean of the data set, the percentiles lines, marking and coloring ofsubsets of the data set, theirs geometric means, notation of individual data in the set . . .dealing with zeros and missing values in compositional data sets with R procedures for simpleand multiplicative replacement strategy,the time series analysis of compositional data.We’ll present the current status of MixeR development and illustrate its use on selected data sets

Systematic survey of clonal complexity in tuberculosis at a populational level and detailed characterization of the isolates involved.

Clonally complex infections by Mycobacterium tuberculosis are progressively more accepted. Studies of their dimension in epidemiological scenarios where the infective pressure is not high are scarce. Our study systematically searched for clonally complex infections (mixed infections by more than one strain and simultaneous presence of clonal variants) by applying mycobacterial interspersed repetitive-unit (MIRU)-variable-number tandem-repeat (VNTR) analysis to M. tuberculosis isolates from two population-based samples of respiratory (703 cases) and respiratory-extrapulmonary (R+E) tuberculosis (TB) cases (71 cases) in a context of moderate TB incidence. Clonally complex infections were found in 11 (1.6%) of the respiratory TB cases and in 10 (14.1%) of those with R+E TB. Among the 21 cases with clonally complex TB, 9 were infected by 2 independent strains and the remaining 12 showed the simultaneous presence of 2 to 3 clonal variants. For the 10 R+E TB cases with clonally complex infections, compartmentalization (different compositions of strains/clonal variants in independent infected sites) was found in 9 of them. All the strains/clonal variants were also genotyped by IS6110-based restriction fragment length polymorphism analysis, which split two MIRU-defined clonal variants, although in general, it showed a lower discriminatory power to identify the clonal heterogeneity revealed by MIRU-VNTR analysis. The comparative analysis of IS6110 insertion sites between coinfecting clonal variants showed differences in the genes coding for a cutinase, a PPE family protein, and two conserved hypothetical proteins. Diagnostic delay, existence of previous TB, risk for overexposure, and clustered/orphan status of the involved strains were analyzed to propose possible explanations for the cases with clonally complex infections. Our study characterizes in detail all the clonally complex infections by M. tuberculosis found in a systematic survey and contributes to the characterization that these phenomena can be found to an extent higher than expected, even in an unselected population-based sample lacking high infective pressure.

Effects of highly active antiretroviral therapy with nelfinavir in vertically HIV-1 infected children: 3 years of follow-up. Long-term response to nelfinavir in children.

BACKGROUND Antiretroviral treatment (ART) in children has special features and consequently, results obtained from clinical trials with antiretroviral drugs in adults may not be representative of children. Nelfinavir (NFV) is an HIV-1 Protease Inhibitor (PI) which has become as one of the first choices of PI for ART in children. We studied during a 3-year follow-up period the effects of highly active antiretroviral therapy with nelfinavir in vertically HIV-1 infected children. METHODS Forty-two vertically HIV-infected children on HAART with NFV were involved in a multicentre prospective study. The children were monitored at least every 3 months with physical examinations, and blood sample collection to measure viral load (VL) and CD4+ cell count. We performed a logistic regression analysis to determinate the odds ratio of baseline characteristics on therapeutic failure. RESULTS Very important increase in CD4+ was observed and VL decreased quickly and it remained low during the follow-up study. Children with CD4+ <25% at baseline achieved CD4+ >25% at 9 months of follow-up. HIV-infected children who achieved undetectable viral load (uVL) were less than 40% in each visit during follow-up. Nevertheless, HIV-infected children with VL >5000 copies/ml were less than 50% during the follow-up study. Only baseline VL was an important factor to predict VL control during follow-up. Virological failure at defined end-point was confirmed in 30/42 patients. Along the whole of follow-up, 16/42 children stopped HAART with NFV. Baseline characteristics were not associated with therapeutic change. CONCLUSION NFV is a safe drug with a good profile and able to achieve an adequate response in children.

Angiogenesis-related gene expression profile with independent prognostic value in advanced ovarian carcinoma.

BACKGROUND Ovarian carcinoma is the most important cause of gynecological cancer-related mortality in Western societies. Despite the improved median overall survival in patients receiving chemotherapy regimens such as paclitaxel and carboplatin combination, relapse still occurs in most advanced diseased patients. Increased angiogenesis is associated with rapid recurrence and decreased survival in ovarian cancer. This study was planned to identify an angiogenesis-related gene expression profile with prognostic value in advanced ovarian carcinoma patients. METHODOLOGY/PRINCIPAL FINDINGS RNAs were collected from formalin-fixed paraffin-embedded samples of 61 patients with III/IV FIGO stage ovarian cancer who underwent surgical cytoreduction and received a carboplatin plus paclitaxel regimen. Expression levels of 82 angiogenesis related genes were measured by quantitative real-time polymerase chain reaction using TaqMan low-density arrays. A 34-gene-profile which was able to predict the overall survival of ovarian carcinoma patients was identified. After a leave-one-out cross validation, the profile distinguished two groups of patients with different outcomes. Median overall survival and progression-free survival for the high risk group was 28.3 and 15.0 months, respectively, and was not reached by patients in the low risk group at the end of follow-up. Moreover, the profile maintained an independent prognostic value in the multivariate analysis. The hazard ratio for death was 2.3 (95% CI, 1.5 to 3.2; p<0.001). CONCLUSIONS/SIGNIFICANCE It is possible to generate a prognostic model for advanced ovarian carcinoma based on angiogenesis-related genes using formalin-fixed paraffin-embedded samples. The present results are consistent with the increasing weight of angiogenesis genes in the prognosis of ovarian carcinoma.

News from compositions, the R package

The R-package “compositions”is a tool for advanced compositional analysis. Its basicfunctionality has seen some conceptual improvement, containing now some facilitiesto work with and represent ilr bases built from balances, and an elaborated subsys-tem for dealing with several kinds of irregular data: (rounded or structural) zeroes,incomplete observations and outliers. The general approach to these irregularities isbased on subcompositions: for an irregular datum, one can distinguish a “regular” sub-composition (where all parts are actually observed and the datum behaves typically)and a “problematic” subcomposition (with those unobserved, zero or rounded parts, orelse where the datum shows an erratic or atypical behaviour). Systematic classificationschemes are proposed for both outliers and missing values (including zeros) focusing onthe nature of irregularities in the datum subcomposition(s).To compute statistics with values missing at random and structural zeros, a projectionapproach is implemented: a given datum contributes to the estimation of the desiredparameters only on the subcompositon where it was observed. For data sets withvalues below the detection limit, two different approaches are provided: the well-knownimputation technique, and also the projection approach.To compute statistics in the presence of outliers, robust statistics are adapted to thecharacteristics of compositional data, based on the minimum covariance determinantapproach. The outlier classification is based on four different models of outlier occur-rence and Monte-Carlo-based tests for their characterization. Furthermore the packageprovides special plots helping to understand the nature of outliers in the dataset.Keywords: coda-dendrogram, lost values, MAR, missing data, MCD estimator,robustness, rounded zeros

Evaluation of the GlideScope (R) for tracheal intubation in patients with cervical spine immobilization by a semi-rigid collar

Background: In patients with cervical spine injury, a cervical collar may prevent cervical spine movements but renders tracheal intubation with a standard laryngoscope difficult if not impossible. We hypothesized that despite the presence of a semi-rigid cervical collar and with the patient's head taped to the trolley, we would be able to intubate all patients with the GlideScopeR and its dedicated stylet. Methods: 50 adult patients (ASA 1 or 2, BMI ≤35 kg/m2) scheduled for elective surgical procedures requiring tracheal intubation were included. After standardized induction of general anesthesia and neuromuscular blockade, the neck was immobilized with an appropriately sized semi-rigid Philadelphia Patriot® cervical collar, the head was taped to the trolley. Laryngoscopy was attempted using a Macintosh laryngoscope blade 4 and the modified Cormack Lehane grade was noted. Subsequently, laryngoscopy with the GlideScopeR was graded and followed by oro-tracheal intubation. Results: All patients were successfully intubated with the GlideScopeR and its dedicated stylet. The median intubation time was 50 sec [43; 61]. The modified Cormack Lehane grade was 3 or 4 at direct laryngoscopy. It was significantly reduced with the GlideScopeR (p <0.0001), reaching 2a in most of patients. Maximal mouth opening was significantly reduced with the cervical collar applied, 4.5 cm [4.5; 5.0] vs. 2.0 cm [1.8; 2.0] (p <0.0001). Conclusions: The GlideScope® allows oro-tracheal intubation in patients having their cervical spine immobilized by a semi-rigid collar and their head taped to the trolley. It furthermore decreases significantly the modified Cormack Lehane grade.

Executive functions profile in extreme eating/weight conditions: from anorexia nervosa to obesity.

BACKGROUND Extreme weight conditions (EWC) groups along a continuum may share some biological risk factors and intermediate neurocognitive phenotypes. A core cognitive trait in EWC appears to be executive dysfunction, with a focus on decision making, response inhibition and cognitive flexibility. Differences between individuals in these areas are likely to contribute to the differences in vulnerability to EWC. The aim of the study was to investigate whether there is a common pattern of executive dysfunction in EWC while comparing anorexia nervosa patients (AN), obese subjects (OB) and healthy eating/weight controls (HC). METHODS Thirty five AN patients, fifty two OB and one hundred thirty seven HC were compared using the Wisconsin Card Sorting Test (WCST); Stroop Color and Word Test (SCWT); and Iowa Gambling Task (IGT). All participants were female, aged between 18 and 60 years. RESULTS There was a significant difference in IGT score (F(1.79); p<.001), with AN and OB groups showing the poorest performance compared to HC. On the WCST, AN and OB made significantly more errors than controls (F(25.73); p<.001), and had significantly fewer correct responses (F(2.71); p<.001). Post hoc analysis revealed that the two clinical groups were not significantly different from each other. Finally, OB showed a significant reduced performance in the inhibition response measured with the Stroop test (F(5.11); p<.001) compared with both AN and HC. CONCLUSIONS These findings suggest that EWC subjects (namely AN and OB) have similar dysfunctional executive profile that may play a role in the development and maintenance of such disorders.

Kallikrein genes are associated with lupus and glomerular basement membrane-specific antibody-induced nephritis in mice and humans.

Immune-mediated nephritis contributes to disease in systemic lupus erythematosus, Goodpasture syndrome (caused by antibodies specific for glomerular basement membrane [anti-GBM antibodies]), and spontaneous lupus nephritis. Inbred mouse strains differ in susceptibility to anti-GBM antibody-induced and spontaneous lupus nephritis. This study sought to clarify the genetic and molecular factors that maybe responsible for enhanced immune-mediated renal disease in these models. When the kidneys of 3 mouse strains sensitive to anti-GBM antibody-induced nephritis were compared with those of 2 control strains using microarray analysis, one-fifth of the underexpressed genes belonged to the kallikrein gene family,which encodes serine esterases. Mouse strains that upregulated renal and urinary kallikreins exhibited less evidence of disease. Antagonizing the kallikrein pathway augmented disease, while agonists dampened the severity of anti-GBM antibody-induced nephritis. In addition, nephritis-sensitive mouse strains had kallikrein haplotypes that were distinct from those of control strains, including several regulatory polymorphisms,some of which were associated with functional consequences. Indeed, increased susceptibility to anti-GBM antibody-induced nephritis and spontaneous lupus nephritis was achieved by breeding mice with a genetic interval harboring the kallikrein genes onto a disease-resistant background. Finally, both human SLE and spontaneous lupus nephritis were found to be associated with kallikrein genes, particularly KLK1 and the KLK3 promoter, when DNA SNPs from independent cohorts of SLE patients and controls were compared. Collectively, these studies suggest that kallikreins are protective disease-associated genes in anti-GBM antibody-induced nephritis and lupus.

Measles outbreak in Andalusia, Spain, January to August 2011.

On 7 January 2011, a six year-old child living in a Roma community near Seville, southern Spain, was hospitalised with measles. Contact tracing identified a probable index case with onset of symptoms on 20 December 2011 and several unreported cases among children under the age of 15 years in the same town. The outbreak initially spread in districts in the city of Seville with a high proportion of Roma residents, and later to other cities and towns in Andalusia. While some towns experienced wide spread of the disease with significant clusters of cases, most of the affected locations saw non-clustered cases or very few secondary cases. The outbreak resulted in 1,759 confirmed or probable cases of which 393 (19%) required hospitalisation. Measles virus of genotype D4 was diagnosed in more than half of the cases. Significant differences (p<0.0001) by age group were found between clustered and non-clustered cases. The highest proportion of clustered cases occurred in the age group of 5-14 yearolds, while the highest proportion of non-clustered cases was seen in those older than 29 years. The last confirmed case related to this outbreak was reported on 20 August 2011.

Dietary fibre intake and risks of cancers of the colon and rectum in the European prospective investigation into cancer and nutrition (EPIC).

BACKGROUND Earlier analyses within the EPIC study showed that dietary fibre intake was inversely associated with colorectal cancer risk, but results from some large cohort studies do not support this finding. We explored whether the association remained after longer follow-up with a near threefold increase in colorectal cancer cases, and if the association varied by gender and tumour location. METHODOLOGY/PRINCIPAL FINDINGS After a mean follow-up of 11.0 years, 4,517 incident cases of colorectal cancer were documented. Total, cereal, fruit, and vegetable fibre intakes were estimated from dietary questionnaires at baseline. Hazard ratios (HRs) and 95% confidence intervals (CIs) were estimated using Cox proportional hazards models stratified by age, sex, and centre, and adjusted for total energy intake, body mass index, physical activity, smoking, education, menopausal status, hormone replacement therapy, oral contraceptive use, and intakes of alcohol, folate, red and processed meats, and calcium. After multivariable adjustments, total dietary fibre was inversely associated with colorectal cancer (HR per 10 g/day increase in fibre 0.87, 95% CI: 0.79-0.96). Similar linear associations were observed for colon and rectal cancers. The association between total dietary fibre and risk of colorectal cancer risk did not differ by age, sex, or anthropometric, lifestyle, and dietary variables. Fibre from cereals and fibre from fruit and vegetables were similarly associated with colon cancer; but for rectal cancer, the inverse association was only evident for fibre from cereals. CONCLUSIONS/SIGNIFICANCE Our results strengthen the evidence for the role of high dietary fibre intake in colorectal cancer prevention.

Implementació dels mòduls SD i MM del sistema SAP R/3 en una empresa del sector metal·lúrgic

SAP R/3 és un software estàndard de gestió empresarial que engloba tots els processos iactivitats d’una companyia i permet el disseny de l’aplicació a les necessitats de l’empresa.Diferents avantatges, com ara disposar d’una base de dades única, han fet escollir aquestsistema a una empresa del sector metal·lúrgic com a eina de gestió:Aquest projecte té com a objecte la implantació de dos dels mòduls que inclou el sistemaSAP R/3, concretament el mòdul SD (Sales and Distribution), amb el qual es controlen lesvendes de l’empresa, i el mòdul MM (Material Management), mòdul amb el qual es controlatot el que fa referència als materials, la gestió d’estocs i tota la gestió de compres del’empresa