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

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

Influence of the ankle arthropathies on the 3D kinematics of the lower limbs during gait: analysis with an ambulatory system.

Introduction: Ankle arthropathy is associated with a decreased motion of the ankle-hindfoot during ambulation. Ankle arthrodesis was shown to result in degeneration of the neighbour joints of the foot. Inversely, total ankle arthroplasty conceptually preserves the adjacent joints because of the residual mobility of the ankle but this has not been demonstrated yet in vivo. It has also been reported that degenerative ankle diseases, and even arthrodesis, do not result in alteration of the knee and hip joints. We present the preliminary results of a new approach of this problem based on ambulatory gait analysis. Patients and Methods: Motion analysis of the lower limbs was performed using a Physilog® (BioAGM, CH) system consisting of three-dimensional (3D) accelerometer and gyroscope, coupled to a magnetic system (Liberty©, Polhemus, USA). Both systems have been validated. Three groups of two patients were included into this pilot study and compared to healthy subjects (controls) during level walking: patients with ankle osteoarthritis (group 1), patients treated by ankle arthrodesis (group 2), patients treated by total ankle prosthesis (group 3). Results: Motion patterns of all analyzed joints over more than 20 gait cycles in each subject were highly repeatable. Motion amplitude of the ankle-hindfoot in control patients was similar to recently reported results. Ankle arthrodesis limited the motion of the ankle-hindfoot in the sagittal and horizontal planes. The prosthetic ankle allowed a more physiologic movement in the sagittal plane only. Ankle arthritis and its treatments did not influence the range of motion of the knee and hip joint during stance phase, excepted for a slight decrease of the hip flexion in groups 1 and 2. Conclusion: The reliability of the system was shown by the repeatability of the consecutive measurements. The results of this preliminary study were similar to those obtained through laboratory gait analysis. However, our system has the advantage to allow ambulatory analysis of 3D kinematics of the lower limbs outside of a gait laboratory and in real life conditions. To our knowledge this is a new concept in the analysis of ankle arthropathy and its treatments. Therefore, there is a potential to address specific questions like the difficult comparison of the benefits of ankle arthroplasty versus arthrodesis. The encouraging results of this pilot study offer the perspective to analyze the consequences of ankle arthropathy and its treatments on the biomechanics of the lower limbs ambulatory, in vivo and in daily life conditions.

Risk Evaluation for Spoofing against a Sensor Supplied with Liveness Detection

The aim of this paper is to evaluate the risks associated with the use of fake fingerprints on a livescan supplied with a method of liveness detection. The method is based on optical properties of the skin. The sensor uses several polarizations and illuminations to capture the information of the different layers of the human skin. These experiments also allow for the determination under which conditions the system is deceived and if there is an influence respectively of the nature of the fake, the mould used for the production or the individuals involved in the attack. These experiments showed that current multispectral sensors can be deceived by the use of fake fingerprints created with or without the cooperation of the subject. Fakes created from direct casts perform better than those produced by fakes created from indirect casts. The results showed that the success of the attack is influenced by two main factors. The first is the quality of the fakes, and by extension the quality of the original fingerprint. The second is the combination of the general patterns involved in the attacks since an appropriate combination can strongly increase the rates of successful attacks.

Computation of multiple correspondence analysis, with code in R

The generalization of simple correspondence analysis, for two categorical variables, to multiple correspondence analysis where they may be three or more variables, is not straighforward, both from a mathematical and computational point of view. In this paper we detail the exact computational steps involved in performing a multiple correspondence analysis, including the special aspects of adjusting the principal inertias to correct the percentages of inertia, supplementary points and subset analysis. Furthermore, we give the algorithm for joint correspondence analysis where the cross-tabulations of all unique pairs of variables are analysed jointly. The code in the R language for every step of the computations is given, as well as the results of each computation.

Perceptual mapping of practical ethics along the value chain: A multiple correspondence analysis with industry and cultural indices as supplementary variables

This paper presents findings from a study investigating a firm s ethical practices along the value chain. In so doing we attempt to better understand potential relationships between a firm s ethical stance with its customers and those of its suppliers within a supply chain and identify particular sectoral and cultural influences that might impinge on this. Drawing upon a database comprising of 667 industrial firms from 27 different countries, we found that ethical practices begin with the firm s relationship with its customers, the characteristics of which then influence the ethical stance with the firm s suppliers within the supply chain. Importantly, market structure along with some key cultural characteristics were also found to exert significant influence on the implementation of ethical policies in these firms.

Enantiomers' potential in psychopharmacology--a critical analysis with special emphasis on the antidepressant escitalopram

Stereochemistry is now influencing most areas of pharmacotherapy, with a growing awareness in the field of psychiatry and, more specifically, depression. This is due to the fact that the enantiomers of many chiral drugs may have distinct pharmacological, pharmacokinetic and/or pharmacogenetic profiles. Consequently, in some instances there may be an advantage in using a single enantiomer over the racemic form-thus providing a basis for the development of new therapeutic agents, as well as the potential to improve current treatments. This review highlights some of the potential advantages and disadvantages that using single enantiomers might offer. The principles are exemplified through reference to the stereoselective properties of several established chiral psychotropic drugs, including thioridazine, methadone, trimipramine, mianserin, mirtazapine, fluoxetine and citalopram. Emphasis is given to the treatment of depression and how the potential of one pure enantiomer-escitalopram, the S-enantiomer of the selective serotonin reuptake inhibitor citalopram-appears to be fulfilling its preclinical promise in the clinic.

Risk evaluation for spoofing against a sensor supplied with liveness detection

The aim of this paper is to evaluate the risks associated with the use of fake fingerprints on a livescan supplied with a method of liveness detection. The method is based on optical properties of the skin. The sensor uses several polarizations and illuminations to capture the information of the different layers of the human skin. These experiments also allow for the determination under which conditions the system is deceived and if there is an influence respectively of the nature of the fake, the mould used for the production or the individuals involved in the attack. These experiments showed that current multispectral sensors can be deceived by the use of fake fingerprints created with or without the cooperation of the subject. Fakes created from direct casts perform better than those produced by fakes created from indirect casts. The results showed that the success of the attack is influenced by two main factors. The first is the quality of the fakes, and by extension the quality of the original fingerprint. The second is the combination of the general patterns involved in the attacks since an appropriate combination can strongly increase the rates of successful attacks.

Synovial sarcomas usually metastasize after >5 years: a multicenter retrospective analysis with minimum follow-up of 10 years for survivors.

Background: Synovial sarcoma (SS) is a malignant soft tissue sarcoma with a poor prognosis because of late local recurrence and distant metastases. To our knowledge, no studies have minimum follow-up of 10 years that evaluate long-term outcomes for survivors. Patients and methods: Data on 62 patients who had been treated for SS from 1968 to 1999 were studied retrospectively in a multicenter study. Mean follow-up of living patients was 17.2 years and of dead patients 7.7 years. Results: Mean age at diagnosis was 35.4 years (range 6-82 years). Overall survival was 38.7%. The 5-year survival was 74.2%; 10-year survival was 61.2%; and 15-year survival was 46.5%. Fifteen patients (24%) died of disease after 10 years of follow-up. Local recurrence occurred after a mean of 3.6 years (range 0.5-14.9 years) and metastases at a mean of 5.7 years (range 0.5-16.3 years). Only four patients were treated technically correctly with a planned biopsy followed by a wide resection or amputation. Factors associated with significantly worse prognosis included larger tumor size, metastases at the time of diagnosis, high-grade histology, trunk-related disease, and lack of wide resection as primary surgical treatment. Conclusions: In SS, metastases develop late with high mortality. Patients with SS should be followed for >10 years.

High-throughput and sensitive quantitation of plasma catecholamines by ultraperformance liquid chromatography-tandem mass spectrometry using a solid phase microwell extraction plate.

Plasma catecholamines provide a reliable biomarker of sympathetic activity. The low circulating concentrations of catecholamines and analytical interferences require tedious sample preparation and long chromatographic runs to ensure their accurate quantification by HPLC with electrochemical detection. Published or commercially available methods relying on solid phase extraction technology lack sensitivity or require derivatization of catecholamine by hazardous reagents prior to tandem mass spectrometry (MS) analysis. Here, we manufactured a novel 96-well microplate device specifically designed to extract plasma catecholamines prior to their quantification by a new and highly sensitive ultraperformance liquid chromatography-tandem mass spectrometry (UPLC-MS/MS) method. Processing time, which included sample purification on activated aluminum oxide and elution, is less than 1 h per 96-well microplate. The UPLC-MS/MS analysis run time is 2.0 min per sample. This UPLC-MS/MS method does not require a derivatization step, reduces the turnaround time by 10-fold compared to conventional methods used for routine application, and allows catecholamine quantification in reduced plasma sample volumes (50-250 μL, e.g., from children and mice).

Effects of bisoprolol fumarate on left ventricular size, function, and exercise capacity in patients with heart failure: analysis with magnetic resonance myocardial tagging.

BACKGROUND: Recent data suggest that beta-blockers can be beneficial in subgroups of patients with chronic heart failure (CHF). For metoprolol and carvedilol, an increase in ejection fraction has been shown and favorable effects on the myocardial remodeling process have been reported in some studies. We examined the effects of bisoprolol fumarate on exercise capacity and left ventricular volume with magnetic resonance imaging (MRI) and applied a novel high-resolution MRI tagging technique to determine myocardial rotation and relaxation velocity. METHODS: Twenty-eight patients (mean age, 57 +/- 11 years; mean ejection fraction, 26 +/- 6%) were randomized to bisoprolol fumarate (n = 13) or to placebo therapy (n = 15). The dosage of the drugs was titrated to match that of the the Cardiac Insufficiency Bisoprolol Study protocol. Hemodynamic and gas exchange responses to exercise, MRI measurements of left ventricular end-systolic and end-diastolic volumes and ejection fraction, and left ventricular rotation and relaxation velocities were measured before the administration of the drug and 6 and 12 months later. RESULTS: After 1 year, heart rate was reduced in the bisoprolol fumarate group both at rest (81 +/- 12 before therapy versus 61 +/- 11 after therapy; P <.01) and peak exercise (144 +/- 20 before therapy versus 127 +/- 17 after therapy; P <.01), which indicated a reduction in sympathetic drive. No differences were observed in heart rate responses in the placebo group. No differences were observed within or between groups in peak oxygen uptake, although work rate achieved was higher (117.9 +/- 36 watts versus 146.1 +/- 33 watts; P <.05) and exercise time tended to be higher (9.1 +/- 1.7 minutes versus 11.4 +/- 2.8 minutes; P =.06) in the bisoprolol fumarate group. A trend for a reduction in left ventricular end-diastolic volume (-54 mL) and left ventricular end-systolic volume (-62 mL) in the bisoprolol fumarate group occurred after 1 year. Ejection fraction was higher in the bisoprolol fumarate group (25.0 +/- 7 versus 36.2 +/- 9%; P <.05), and the placebo group remained unchanged. Most changes in volume and ejection fraction occurred during the latter 6 months of treatment. With myocardial tagging, insignificant reductions in left ventricular rotation velocity were observed in both groups, whereas relaxation velocity was reduced only after bisoprolol fumarate therapy (by 39%; P <.05). CONCLUSION: One year of bisoprolol fumarate therapy resulted in an improvement in exercise capacity, showed trends for reductions in end-diastolic and end-systolic volumes, increased ejection fraction, and significantly reduced relaxation velocity. Although these results generally confirm the beneficial effects of beta-blockade in patients with chronic heart failure, they show differential effects on systolic and diastolic function.

Smooth height/age curves from stem analysis with linear programming

[Abstract]

Family Size Evolution in Drosophila Chemosensory Gene Families: A Comparative Analysis with a Critical Appraisal of Methods

Gene turnover rates and the evolution of gene family sizes are important aspects of genome evolution. Here, we use curated sequence data of the major chemosensory gene families from Drosophila-the gustatory receptor, odorant receptor, ionotropic receptor, and odorant-binding protein families-to conduct a comparative analysis among families, exploring different methods to estimate gene birth and death rates, including an ad hoc simulation study. Remarkably, we found that the state-of-the-art methods may produce very different rate estimates, which may lead to disparate conclusions regarding the evolution of chemosensory gene family sizes in Drosophila. Among biological factors, we found that a peculiarity of D. sechellia's gene turnover rates was a major source of bias in global estimates, whereas gene conversion had negligible effects for the families analyzed herein. Turnover rates vary considerably among families, subfamilies, and ortholog groups although all analyzed families were quite dynamic in terms of gene turnover. Computer simulations showed that the methods that use ortholog group information appear to be the most accurate for the Drosophila chemosensory families. Most importantly, these results reveal the potential of rate heterogeneity among lineages to severely bias some turnover rate estimation methods and the need of further evaluating the performance of these methods in a more diverse sampling of gene families and phylogenetic contexts. Using branch-specific codon substitution models, we find further evidence of positive selection in recently duplicated genes, which attests to a nonneutral aspect of the gene birth-and-death process.

Mecanização no preparo de amostras por microondas: o estado da arte

This review presents the latest advances of the mechanisation in sample preparation using microwave oven. The fundamental principles of microwaves field interaction with matter are presented, followed by the basic principles of microwave equipment construction. Hyphenation of microwave oven with flow injection analysis or robotic are discussed as well as some metodological difficulties and solutions by its implementation.

Sistema automático para análise direta de ligas metálicas por eletrodissolução anódica e espectrometria de absorção atômica

An automatic system for the direct determination of lead and tin by atomic absorption spectrometry is described. The on-line treatment of the metallic samples was obtained by anodic electrodissolution in a flow injection system. Lead was determined by flame atomic absorption spectrometry (FAAS) and tin by graphite furnace atomic absorption spectrometry (GFAAS). A computer program managed the current source and the solenoid valves that direct the fluids. Good linear correlations between absorbance and current intensity for lead and tin were observed. Results were in agreement with the certified values. Precision was always better than 5%. The recommended procedure allows the direct determination of 60 or 30 elements/h using FAAS or GFAAS, respectively.