Compositional data analysis with R

Examples of compositional data. The simplex, a suitable sample space for compositional data and Aitchison's geometry. R, a free language and environment for statistical computing and graphics

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

Suture repair of umbilical hernia during caesarean section: a case-control study.

PURPOSE: The objective of this study was to investigate the additional burdens in terms of pain, prolongation of surgery and morbidity which is added to elective caesarean section if umbilical hernia suture repair is performed simultaneously. Secondly, patient's satisfaction and hernia recurrence rate were assessed. METHODS: Consecutive women with symptomatic umbilical hernia undergoing internal or external suture repair during elective caesarean were included in this retrospective cohort-control study. Data on post-operative pain, duration of surgery and morbidity of a combined procedure were collected. These patients were matched 1:10 to women undergoing caesarean section only. Additionally, two subgroups were assessed separately: external and internal suture hernia repair. These subgroups were compared for patient's satisfaction, cosmesis, body image and recurrence rate. RESULTS: Fourteen patients with a mean age of 37 years were analysed. Internal suture repair (n = 7) prolonged caesarean section by 20 min (p = 0.001) and external suture repair (n = 7) by 34 min (p < 0.0001). Suture repair did not increase morphine use (0.38 ± 0.2 vs. 0.4 ± 02 mg/kg body weight), had no procedure-related morbidity and prolonged hospitalization by 0.5 days (p = 0.01). At a median follow-up of 37 (5-125) months, two recurrences in each surgical technique, internal and external suture repair, occurred (28 %). Body image and cosmesis score showed a higher level of functioning in internal suture repair (p = 0.02; p = 0.04). DISCUSSION: Despite a high recurrence rate, internal suture repair of a symptomatic umbilical hernia during elective caesarean section should be offered to women if requested. No additional morbidity or scar is added to caesarean section. Internal repair is faster, and cosmetic results are better, additional skin or fascia dissection is avoided, and it seems to be as effective as an external approach. Yet, women must be informed on the high recurrence rate.

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

Biased thinking assessed by external observers in borderline personality disorder

Objectives. Biased thinking (to some extent overlapping with the concepts of cognitive distortions and cognitive errors) is a key concept in cognitive therapy of Borderline Personality Disorder (BPD). Specific contents and cognitive processes related to BPD functioning are known. However, most studies are based on self-report measures which present a number of important limitations, in particular the difficulty in assessing non-conscious processes infused by affect. So far, no studies were conducted using valid observer-rated methodology addressing the question of biased thinking in BPD as it unfolds spontaneously in session. Design. This is a controlled interview study comparing two matched groups, BPD patients and healthy controls. Methods. A total of N= 25 clinical dynamic interviews with patients presenting with BPD were transcribed and rated using the Cognitive Errors Rating Scale (Drapeau, Perry, & Dunkley, 2008); their cognitive profiles were compared to those of N= 25 healthy controls who underwent the same procedure. Results. Overall, results indicated that no between-group difference in the frequency of specific biases was found. However, heightened levels of negative cognitive biases, in particular over-generalizing and fortune-telling, were associated with BPD. Furthermore, negative over-generalizing was associated with the number of BPD symptoms. Conclusions. These results have high levels of ecological validity and are promising for the refinement of cognitive theory of BPD. Clinical implications for assessment and intervention are discussed.

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.

A new perspective on the zoogeography of the sibling mouse-eared bat species Myotis myotis and Myotis blythii: Morphological, genetical and ecological evidence

The actual geographic distribution of the two sibling mouse-eared bat species Myotis myotis and Myotis blythii, which occur widely sympatrically in the western Palaearctic region, remains largely controversial. This concerns particularly the specific attribution of marginal populations from the Mediterranean islands and from adjacent areas of North Africa and Asia, which are morphologically intermediate between continental M. myotis and M. blythii from Europe. This study attempts to clarify this question by using four different approaches: cranial morphology, external morphology, genetics and trophic ecology. The three latter methods show unambiguously that North Africa, Malta, Sardinia and Corsica are presently inhabited by monospecific populations of M. myotis. In contrast, cranial morphometrics do not yield conclusive results. These results contradict all recent studies, which attribute North African and Maltese mouse-eared bats to M. blythii and consider that Sardinia and Corsica harbour sympatric populations of the two species. As concerns south-eastern populations, doubts are also expressed about the attribution of the subspecific taxon omari which may actually refer to M. myotis instead of M. blythii. Protein electrophoresis is presently the only absolute method available for determining M. myotis and M. blythii throughout their distribution ranges. However, species identification may be approached by relying on less sophisticated morphometrical methods as presented in this study. Species-specific habitat specializations are probably responsible for the differences observed between the geographic distributions of M. myotis and M. blythii, as they provide a logical groundwork for a coherent model of speciation for these two bat species.

Altered adaptive but not veridical decision-making in substance dependent individuals.

Drug addiction is associated with impaired judgment in unstructured situations in which success depends on self-regulation of behavior according to internal goals (adaptive decision-making). However most executive measures are aimed at assessing decision-making in structured scenarios, in which success is determined by external criteria inherent to the situation (veridical decision-making). The aim of this study was to examine the performance of Substance Abusers (SA, n = 97) and Healthy Comparison participants (HC, n = 81) in two behavioral tasks that mimic the uncertainty inherent in real-life decision-making: the Cognitive Bias Task (CB) and the Iowa Gambling Task (IGT) (administered only to SA). A related goal was to study the interdependence between performances on both tasks. We conducted univariate analyses of variance (ANOVAs) to contrast the decision-making performance of both groups; and used correlation analyses to study the relationship between both tasks. SA showed a marked context-independent decision-making strategy on the CB's adaptive condition, but no differences were found on the veridical conditions in a subsample of SA (n = 34) and HC (n = 22). A high percentage of SA (75%) also showed impaired performance on the IGT. Both tasks were only correlated when no impaired participants were selected. Results indicate that SA show abnormal decision-making performance in unstructured situations, but not in veridical situations.

A First approximation to the realism/idealism debate about the external world

El text intenta fer una primera aproximació al debat contemporani entre realistes i anti-realistes sobre el món empíric, centrant-se en les posicions de Putnam i Nagel. El seu objectiu principal és el d'entendre les motivacions de les posicions i l'estructura actual del debat, i el d'establir les característiques que hauria de tenir qualsevol posició satisfactòria

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

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