Distributions on the simplex

The simplex, the sample space of compositional data, can be structured as a real Euclidean space. This fact allows to work with the coefficients with respect to an orthonormal basis. Over these coefficients we apply standard real analysis, inparticular, we define two different laws of probability trought the density function and we study their main properties

Meaning of the λ parameter of skew–normal and log–skew normal distributions in fluid geochemistry

The literature related to skew–normal distributions has grown rapidly in recent yearsbut at the moment few applications concern the description of natural phenomena withthis type of probability models, as well as the interpretation of their parameters. Theskew–normal distributions family represents an extension of the normal family to whicha parameter (λ) has been added to regulate the skewness. The development of this theoreticalfield has followed the general tendency in Statistics towards more flexible methodsto represent features of the data, as adequately as possible, and to reduce unrealisticassumptions as the normality that underlies most methods of univariate and multivariateanalysis. In this paper an investigation on the shape of the frequency distribution of thelogratio ln(Cl−/Na+) whose components are related to waters composition for 26 wells,has been performed. Samples have been collected around the active center of Vulcanoisland (Aeolian archipelago, southern Italy) from 1977 up to now at time intervals ofabout six months. Data of the logratio have been tentatively modeled by evaluating theperformance of the skew–normal model for each well. Values of the λ parameter havebeen compared by considering temperature and spatial position of the sampling points.Preliminary results indicate that changes in λ values can be related to the nature ofenvironmental processes affecting the data

Methods for testing association between uncertain genotypes and quantitative traits.

Interpretability and power of genome-wide association studies can be increased by imputing unobserved genotypes, using a reference panel of individuals genotyped at higher marker density. For many markers, genotypes cannot be imputed with complete certainty, and the uncertainty needs to be taken into account when testing for association with a given phenotype. In this paper, we compare currently available methods for testing association between uncertain genotypes and quantitative traits. We show that some previously described methods offer poor control of the false-positive rate (FPR), and that satisfactory performance of these methods is obtained only by using ad hoc filtering rules or by using a harsh transformation of the trait under study. We propose new methods that are based on exact maximum likelihood estimation and use a mixture model to accommodate nonnormal trait distributions when necessary. The new methods adequately control the FPR and also have equal or better power compared to all previously described methods. We provide a fast software implementation of all the methods studied here; our new method requires computation time of less than one computer-day for a typical genome-wide scan, with 2.5 M single nucleotide polymorphisms and 5000 individuals.

Matrix-assisted laser desorption ionization time-of-flight mass spectrometry and PCR-based rapid diagnosis of Staphylococcus aureus bacteraemia.

Effective empirical treatment is of paramount importance to improve the outcome of patients with Staphylococcus aureus bacteraemia. We aimed to evaluate a PCR-based rapid diagnosis of methicillin resistance (GeneXpert MRSA) after early detection of S. aureus bacteraemia using matrix-assisted laser desorption ionization time-of-flight mass spectrometry (MALDI-TOF MS). Patients with a first episode of S. aureus bacteraemia identified using MALDI-TOF MS were randomized in a prospective interventional open study between October 2010 and August 2012. In the control group, antibiotic susceptibility testing was performed after MALDI-TOF MS identification on blood culture pellets. In the intervention group, a GeneXpert MRSA was performed after S. aureus identification. The primary outcome was the performance of GeneXpert MRSA directly on blood cultures. We then assessed the impact of early diagnosis of methicillin resistance on the empirical treatment. In all, 197 episodes of S. aureus bacteraemia were included in the study, of which 106 were included in the intervention group. Median time from MALDI-TOF MS identification to GeneXpert MRSA result was 97 min (range 25-250). Detection of methicillin resistance using GeneXpert MRSA had a sensitivity of 99% and a specificity of 100%. There was less unnecessary coverage of MRSA in the intervention group (17.1% versus 29.2%, p 0.09). GeneXpert MRSA was highly reliable in diagnosing methicillin resistance when performed directly on positive blood cultures. This could help to avoid unnecessary prescriptions of anti-MRSA agents and promote the introduction of earlier adequate coverage in unsuspected cases.

A compositional statistical analysis of capital stock

Most of economic literature has presented its analysis under the assumption of homogeneous capital stock.However, capital composition differs across countries. What has been the pattern of capital compositionassociated with World economies? We make an exploratory statistical analysis based on compositional datatransformed by Aitchinson logratio transformations and we use tools for visualizing and measuring statisticalestimators of association among the components. The goal is to detect distinctive patterns in the composition.As initial findings could be cited that:1. Sectorial components behaved in a correlated way, building industries on one side and , in a lessclear view, equipment industries on the other.2. Full sample estimation shows a negative correlation between durable goods component andother buildings component and between transportation and building industries components.3. Countries with zeros in some components are mainly low income countries at the bottom of theincome category and behaved in a extreme way distorting main results observed in the fullsample.4. After removing these extreme cases, conclusions seem not very sensitive to the presence ofanother isolated cases

A survey addressing the fundamental matrix estimation problem

Epipolar geometry is a key point in computer vision and the fundamental matrix estimation is the only way to compute it. This article surveys several methods of fundamental matrix estimation which have been classified into linear methods, iterative methods and robust methods. All of these methods have been programmed and their accuracy analysed using real images. A summary, accompanied with experimental results, is given

Is the ability of biocontrol fluorescent pseudomonads to produce the antifungal metabolite 2,4-diacetylphloroglucinol really synonymous with higher plant protection?

The antifungal compound 2,4-diacetylphloroglucinol (Phl) contributes to biocontrol in pseudomonads, but whether or not Phl(+) biocontrol pseudomonads display higher plant-protecting activity than Phl(-) biocontrol pseudomonads remains to be demonstrated. This issue was addressed by assessing 230 biocontrol fluorescent pseudomonads selected from a collection of 3132 bacterial isolates obtained from 63 soils worldwide. One-third of the biocontrol pseudomonads were Phl(+) and almost all Phl(+) isolates also produced hydrogen cyanide (HCN). The only Phl(+) HCN(-) strain did harbor hcn genes, but with the deletion of a 134 bp hcnC fragment corresponding to an ADP-binding motif. Statistical analysis of biocontrol isolate distributions indicated that Phl production ability was associated with superior disease suppression activity in the Pythium-cucumber and Fusarium-tomato pathosystems, but this was also the case with HCN production ability. However, HCN significance was not as strong, as indicated both by the comparison of Phl(-) HCN(+) and Phl(-) HCN(-) strains and by correlation analyses. This is the first population-level demonstration of the higher plant-protecting activity of Phl(+) biocontrol pseudomonads in comparison with Phl(-) biocontrol pseudomonads.

Predicting spatial patterns of plant species richness: a comparison of direct macroecological and species stacking approaches

Aim This study compares the direct, macroecological approach (MEM) for modelling species richness (SR) with the more recent approach of stacking predictions from individual species distributions (S-SDM). We implemented both approaches on the same dataset and discuss their respective theoretical assumptions, strengths and drawbacks. We also tested how both approaches performed in reproducing observed patterns of SR along an elevational gradient.Location Two study areas in the Alps of Switzerland.Methods We implemented MEM by relating the species counts to environmental predictors with statistical models, assuming a Poisson distribution. S-SDM was implemented by modelling each species distribution individually and then stacking the obtained prediction maps in three different ways - summing binary predictions, summing random draws of binomial trials and summing predicted probabilities - to obtain a final species count.Results The direct MEM approach yields nearly unbiased predictions centred around the observed mean values, but with a lower correlation between predictions and observations, than that achieved by the S-SDM approaches. This method also cannot provide any information on species identity and, thus, community composition. It does, however, accurately reproduce the hump-shaped pattern of SR observed along the elevational gradient. The S-SDM approach summing binary maps can predict individual species and thus communities, but tends to overpredict SR. The two other S-SDM approaches the summed binomial trials based on predicted probabilities and summed predicted probabilities - do not overpredict richness, but they predict many competing end points of assembly or they lose the individual species predictions, respectively. Furthermore, all S-SDM approaches fail to appropriately reproduce the observed hump-shaped patterns of SR along the elevational gradient.Main conclusions Macroecological approach and S-SDM have complementary strengths. We suggest that both could be used in combination to obtain better SR predictions by following the suggestion of constraining S-SDM by MEM predictions.

Modern statistical models for forensic fingerprint examinations : a critical review

Over the last decade, the development of statistical models in support of forensic fingerprint identification has been the subject of increasing research attention, spurned on recently by commentators who claim that the scientific basis for fingerprint identification has not been adequately demonstrated. Such models are increasingly seen as useful tools in support of the fingerprint identification process within or in addition to the ACE-V framework. This paper provides a critical review of recent statistical models from both a practical and theoretical perspective. This includes analysis of models of two different methodologies: Probability of Random Correspondence (PRC) models that focus on calculating probabilities of the occurrence of fingerprint configurations for a given population, and Likelihood Ratio (LR) models which use analysis of corresponding features of fingerprints to derive a likelihood value representing the evidential weighting for a potential source.

CODAPACK3D. A new version of Compositional Data Package

The statistical analysis of compositional data should be treated using logratios of parts,which are difficult to use correctly in standard statistical packages. For this reason afreeware package, named CoDaPack was created. This software implements most of thebasic statistical methods suitable for compositional data.In this paper we describe the new version of the package that now is calledCoDaPack3D. It is developed in Visual Basic for applications (associated with Excel©),Visual Basic and Open GL, and it is oriented towards users with a minimum knowledgeof computers with the aim at being simple and easy to use.This new version includes new graphical output in 2D and 3D. These outputs could bezoomed and, in 3D, rotated. Also a customization menu is included and outputs couldbe saved in jpeg format. Also this new version includes an interactive help and alldialog windows have been improved in order to facilitate its use.To use CoDaPack one has to access Excel© and introduce the data in a standardspreadsheet. These should be organized as a matrix where Excel© rows correspond tothe observations and columns to the parts. The user executes macros that returnnumerical or graphical results. There are two kinds of numerical results: new variablesand descriptive statistics, and both appear on the same sheet. Graphical output appearsin independent windows. In the present version there are 8 menus, with a total of 38submenus which, after some dialogue, directly call the corresponding macro. Thedialogues ask the user to input variables and further parameters needed, as well aswhere to put these results. The web site http://ima.udg.es/CoDaPack contains thisfreeware package and only Microsoft Excel© under Microsoft Windows© is required torun the software.Kew words: Compositional data Analysis, Software

The ability of haemolysins expressed by atypical enteropathogenic Escherichia coli to bind to extracellular matrix components

Typical and atypical enteropathogenic Escherichia coli (EPEC) are considered important bacterial causes of diarrhoea. Considering the repertoire of virulence genes, atypical EPEC (aEPEC) is a heterogeneous group, harbouring genes that are found in other diarrheagenic E. coli pathotypes, such as those encoding haemolysins. Haemolysins are cytolytic toxins that lyse host cells disrupting the function of the plasma membrane. In addition, these cytolysins mediate a connection to vascular tissue and/or blood components, such as plasma and cellular fibronectin. Therefore, we investigated the haemolytic activity of 72 aEPEC isolates and determined the correlation of this phenotype with the presence of genes encoding enterohaemolysins (Ehly) and cytolysin A (ClyA). In addition, the correlation between the expression of haemolysins and the ability of these secreted proteins to adhere to extracellular matrix (ECM) components was also assessed in this study. Our findings demonstrate that a subset of aEPEC presents haemolytic activity due to the expression of Ehlys and/or ClyA and that this activity is closely related to the ability of these isolates to bind to ECM components.

Applications of stochastic analysis in finance and statistical inference

First: A continuous-time version of Kyle's model (Kyle 1985), known as the Back's model (Back 1992), of asset pricing with asymmetric information, is studied. A larger class of price processes and of noise traders' processes are studied. The price process, as in Kyle's model, is allowed to depend on the path of the market order. The process of the noise traders' is an inhomogeneous Lévy process. Solutions are found by the Hamilton-Jacobi-Bellman equations. With the insider being risk-neutral, the price pressure is constant, and there is no equilibirium in the presence of jumps. If the insider is risk-averse, there is no equilibirium in the presence of either jumps or drifts. Also, it is analised when the release time is unknown. A general relation is established between the problem of finding an equilibrium and of enlargement of filtrations. Random announcement time is random is also considered. In such a case the market is not fully efficient and there exists equilibrium if the sensitivity of prices with respect to the global demand is time decreasing according with the distribution of the random time. Second: Power variations. it is considered, the asymptotic behavior of the power variation of processes of the form _integral_0^t u(s-)dS(s), where S_ is an alpha-stable process with index of stability 0&alpha&2 and the integral is an Itô integral. Stable convergence of corresponding fluctuations is established. These results provide statistical tools to infer the process u from discrete observations. Third: A bond market is studied where short rates r(t) evolve as an integral of g(t-s)sigma(s) with respect to W(ds), where g and sigma are deterministic and W is the stochastic Wiener measure. Processes of this type are particular cases of ambit processes. These processes are in general not of the semimartingale kind.

Disruption of the peritrophic matrix by exogenous chitinase feeding reduces fecundity in Lutzomyia longipalpis females

Lutzomyia longipalpis is the most important vector of visceral leishmaniasis in Brazil. When female sandflies feed on blood, a peritrophic matrix (PM) is formed around the blood bolus. The PM is secreted by midgut cells and composed of proteins, glycoproteins and chitin microfibrils. The PM functions as both a physical barrier against pathogens present in the food bolus and blood meal digestion regulator. Previous studies of mosquitoes and sandflies have shown that the absence of a PM, resulting from adding an exogenous chitinase to the blood meal, accelerates digestion. In the present study, we analysed biological factors associated with the presence of a PM in L. longipalpis females. Insects fed blood containing chitinase (BCC) accelerated egg-laying relative to a control group fed blood without chitinase. However, in the BCC-fed insects, the number of females that died without laying eggs was higher and the number of eggs laid per female was lower. The eggs in both groups were viable and generated adults. Based on these data, we suggest that the absence of a PM accelerates nutrient acquisition, which results in premature egg production and oviposition; however, the absence of a PM reduces the total number of eggs laid per female. Reduced fecundity in the absence of a PM may be due to inefficient nutrient conversion or the loss of the protective role of the PM.

Extracellular matrix molecules enhance the neurotrophic effect of schwann cell-like differentiated adipose-derived stem cells and increase cell survival under stress conditions.

Since the first reports of induction of adipose-derived stem cells (ASC) into neuronal and glial cell phenotypes, expectations have increased regarding their use in tissue engineering applications for nerve repair. Cell adhesion to extracellular matrix (ECM) is a basic feature of survival, differentiation, and migration of Schwann cells (SC) during nerve regeneration, and fibronectin and laminin are two key molecules of this process. Interaction between ECM and SC-like differentiated ASC (dASC) could potentially improve the neurotrophic potential of the stem cells. We have investigated the effect of ECM molecules on SC-like dASC in terms of proliferation, adhesion, and cell viability. Fibronectin and laminin did not affect the proliferation of dASC when compared with cell adherent tissue culture plastic, but significantly improved viability and cell attachment when dASC were exposed to apoptotic conditions. To assess the influence of the ECM molecules on dASC neurotrophic activity, dASC were seeded onto ECM-coated culture inserts suspended above dorsal root ganglia (DRG) sensory neurons. Neurite outgrowth of DRG neurons was enhanced when dASC were seeded on fibronectin and laminin when compared with controls. When DRG neurons and dASC were in direct contact on the various surfaces there was significantly enhanced neurite outgrowth and coculture with laminin-conditioned dASC produced the longest neurites. Compared with primary SCs, dASC grown on laminin produced similar levels of neurite outgrowth in the culture insert experiments but neurite length was shorter in the direct contact groups. Anti β1 integrin blocking antibody could inhibit baseline and dASC evoked neurite elongation but had no effect on outgrowth mediated by laminin-conditioned dASC. ECM molecules had no effect on the levels of nerve growth factor and brain-derived neurotrophic factor secretion from dASC. The results of the study suggest that ECM molecules can significantly improve the potential of dASC for nerve regeneration.