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 computing and graphics. The environment in which many classical and modern statistical techniques have been implemented, but many are supplied as packages. There are 8 standard packages and many more are 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 is a MixeR package for compositional data analysis that provides support for operations on compositions: perturbation and power multiplication, subcomposition with or without residuals, 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 of subsets 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 simple and 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

Multilevel bootstrap analysis with assumptions violated

Resumen tomado de la publicación

Family-based association analysis with ordered categorical phenotypes, covariates and interactions

Genetic association analyses of family-based studies with ordered categorical phenotypes are often conducted using methods either for quantitative or for binary traits, which can lead to suboptimal analyses. Here we present an alternative likelihood-based method of analysis for single nucleotide polymorphism (SNP) genotypes and ordered categorical phenotypes in nuclear families of any size. Our approach, which extends our previous work for binary phenotypes, permits straightforward inclusion of covariate, gene-gene and gene-covariate interaction terms in the likelihood, incorporates a simple model for ascertainment and allows for family-specific effects in the hypothesis test. Additionally, our method produces interpretable parameter estimates and valid confidence intervals. We assess the proposed method using simulated data, and apply it to a polymorphism in the c-reactive protein (CRP) gene typed in families collected to investigate human systemic lupus erythematosus. By including sex interactions in the analysis, we show that the polymorphism is associated with anti-nuclear autoantibody (ANA) production in females, while there appears to be no effect in males.

A practical approach to a multi-level analysis with a sparse binary outcome within a large surgical trial

Objectives: To assess the potential source of variation that surgeon may add to patient outcome in a clinical trial of surgical procedures. Methods: Two large (n = 1380) parallel multicentre randomized surgical trials were undertaken to compare laparoscopically assisted hysterectomy with conventional methods of abdominal and vaginal hysterectomy; involving 43 surgeons. The primary end point of the trial was the occurrence of at least one major complication. Patients were nested within surgeons giving the data set a hierarchical structure. A total of 10% of patients had at least one major complication, that is, a sparse binary outcome variable. A linear mixed logistic regression model (with logit link function) was used to model the probability of a major complication, with surgeon fitted as a random effect. Models were fitted using the method of maximum likelihood in SAS((R)). Results: There were many convergence problems. These were resolved using a variety of approaches including; treating all effects as fixed for the initial model building; modelling the variance of a parameter on a logarithmic scale and centring of continuous covariates. The initial model building process indicated no significant 'type of operation' across surgeon interaction effect in either trial, the 'type of operation' term was highly significant in the abdominal trial, and the 'surgeon' term was not significant in either trial. Conclusions: The analysis did not find a surgeon effect but it is difficult to conclude that there was not a difference between surgeons. The statistical test may have lacked sufficient power, the variance estimates were small with large standard errors, indicating that the precision of the variance estimates may be questionable.

NOx and O3 above a tropical rainforest: an analysis with a global and box model

A cross-platform field campaign, OP3, was conducted in the state of Sabah in Malaysian Borneo between April and July of 2008. Among the suite of observations recorded, the campaign included measurements of NOx and O3 – crucial outputs of any model chemistry mechanism. We describe the measurements of these species made from both the ground site and aircraft. We then use the output from two resolutions of the chemistry transport model p-TOMCAT to illustrate the ability of a global model chemical mechanism to capture the chemistry at the rainforest site. The basic model performance is good for NOx and poor for ozone. A box model containing the same chemical mechanism is used to explore the results of the global model in more depth and make comparisons between the two. Without some parameterization of the nighttime boundary layer – free troposphere mixing (i.e. the use of a dilution parameter), the box model does not reproduce the observations, pointing to the importance of adequately representing physical processes for comparisons with surface measurements. We conclude with a discussion of box model budget calculations of chemical reaction fluxes, deposition and mixing, and compare these results to output from p-TOMCAT. These show the same chemical mechanism behaves similarly in both models, but that emissions and advection play particularly strong roles in influencing the comparison to surface measurements.

Direct Computation of Stochastic Flow in Reservoirs with Uncertain Parameters

A direct method is presented for determining the uncertainty in reservoir pressure, flow, and net present value (NPV) using the time-dependent, one phase, two- or three-dimensional equations of flow through a porous medium. The uncertainty in the solution is modelled as a probability distribution function and is computed from given statistical data for input parameters such as permeability. The method generates an expansion for the mean of the pressure about a deterministic solution to the system equations using a perturbation to the mean of the input parameters. Hierarchical equations that define approximations to the mean solution at each point and to the field covariance of the pressure are developed and solved numerically. The procedure is then used to find the statistics of the flow and the risked value of the field, defined by the NPV, for a given development scenario. This method involves only one (albeit complicated) solution of the equations and contrasts with the more usual Monte-Carlo approach where many such solutions are required. The procedure is applied easily to other physical systems modelled by linear or nonlinear partial differential equations with uncertain data.

The forms of bioavailability of phosphorus in integrated vertical flow constructed wetland with earthworms and different substrates

A sequential extraction method was utilized to analyze seven forms of P in an integrated vertical-flow constructed wetland (IVFCW) containing earthworms and different substrates. The aluminum-bound P (Al-P) content was found to be lower, and the occluded P (Oc-P) content was higher in the IVFCW. The addition of earthworms into the influent chamber of IVFCW increased the exchange P (Ex-P), iron-bound P (Fe-P), calcium bound P (Ca-P), Oc-P, detritus-bound (De-P) and organic P (Org-P) content in the influent chamber, and also enhanced P content uptake by wetland plants. A significantly positive correlation between P content of above-ground wetland plants and the Ex-P, Fe-P, Oc-P and Org-P content in the rhizosphere was found (P < 0.05), which indicated that the Ex-P, Fe-P, Oc-P and Org-P could be bio-available P. The Ex-P, Fe-P, De-P, Oc-P and Ca-P content of the influent chamber was higher where the substrate contained a mixture of Qing sand and river sand rather than only river sand. Also the IVFCW with earthworms and both Qing sand and river sand had a higher removal efficiency of P, which was related to higher P content uptake by wetland plants and P retained in IVFCW. These findings suggest that addition of earthworms in IVFCW increases the bioavailable P content, resulting in enhanced P content uptake by wetland plants.

Immobilization of glucose oxidase enzyme (GOD) in large pore ordered mesoporous cage-like FDU-1 silica

Large pore ordered mesoporous silica FDU-1 with three-dimensional (3D) face-centered cubic, Fm3m arrangement of rnesopores, was synthesized under strong acid media using B-50-6600 poly(ethylene oxide)-poly(butylene oxide)-poly(ethylene oxide) triblock copolymer (EO(39)BO(47)EO(39)), tetraethyl orthosilicate (TEOS) and trimethyl-benzene (TMB). Large pore FDU-1 silica was obtained by using the following gel composition 1TEOS:0.00735B50-6600:0.00735TMB:6HCl:155H(2)O. The pristine material exhibited a BET specific surface area of 684 m(2) g(-1), total pore volume of 0.89 cm(3) g(-1), external surface area of 49 m(2) g(-1) and microporous volume of 0.09 cm(3) g(-1). The enzyme activity was determined by the Flow Injection Analysis-Chemiluminescence (FIA-CL) method. For GOD immobilized on the FDU-1 silica, GOD supernatant and GOD solution, the FIA-CL results were 9.0, 18.6 and 34.0 U, respectively. The value obtained for the activity of the GOD solution with FIA-CL method is in agreement with the 35 U, obtained by spectrophotometry. (C) 2011 Elsevier B.V. All rights reserved.

STABILITY ANALYSIS WITH APPLICATIONS OF A TWO-DIMENSIONAL DYNAMICAL SYSTEM ARISING FROM A STOCHASTIC MODEL FOR AN ASSET MARKET

We analyze the stability properties of equilibrium solutions and periodicity of orbits in a two-dimensional dynamical system whose orbits mimic the evolution of the price of an asset and the excess demand for that asset. The construction of the system is grounded upon a heterogeneous interacting agent model for a single risky asset market. An advantage of this construction procedure is that the resulting dynamical system becomes a macroscopic market model which mirrors the market quantities and qualities that would typically be taken into account solely at the microscopic level of modeling. The system`s parameters correspond to: (a) the proportion of speculators in a market; (b) the traders` speculative trend; (c) the degree of heterogeneity of idiosyncratic evaluations of the market agents with respect to the asset`s fundamental value; and (d) the strength of the feedback of the population excess demand on the asset price update increment. This correspondence allows us to employ our results in order to infer plausible causes for the emergence of price and demand fluctuations in a real asset market. The employment of dynamical systems for studying evolution of stochastic models of socio-economic phenomena is quite usual in the area of heterogeneous interacting agent models. However, in the vast majority of the cases present in the literature, these dynamical systems are one-dimensional. Our work is among the few in the area that construct and study analytically a two-dimensional dynamical system and apply it for explanation of socio-economic phenomena.

Trypanosoma cruzi: establishment of permeable cells for RNA processing analysis with drugs

Pre-mRNA maturation in trypanosomatids occurs through a process called trans-splicing which involves excision of introns and union of exons in two independent transcripts. For the first time, we present the standardization of Trypanosoma cruzi permeable cells (Y strain) as a model for trans-splicing study of mRNAs in trypanosomes, following by RNase protection reaction, which localizes the SL exon and intron. This trans-splicing reaction in vitro was also used to analyze the influence of NFOH-121, a nitrofurazone-derivative, on this mechanism. The results suggested that the prodrug affects the RNA processing in these parasites, but the trans-splicing reaction still occurred.

The Temperature and Humidity in a Low-Flow Anesthesia Workstation With and Without a Heat and Moisture Exchanger

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Synchrotron structural characterization of electrochemically synthesized hexacyanoferrates containing K(+): A revisited analysis of electrochemical redox

The presumably soluble KFe(+3)[Fe(2+)(CN)(6)] structure of electrochemically synthesized hexacyanoferrate materials (Prussian Blue) containing K(+) ions was determined for the first time in this study. Prior to drawing conclusions from a structural analysis, the main goal was to make a precise analysis of the inferred soluble structure, that is, KFe(+3) [Fe(2+)(CN)(6)], which is frequently referred to in the literature as the final stable electrochemically synthesized structure. Indeed, a successful X-ray powder diffraction experiment using X-ray synchrotron radiation was made of a powder placed in a 0.5 mm diameter borosilicate glass capillary, which was obtained by removing sixty 90 nm thin films from the substrates on which they were prepared. However, the conclusions were highly unexpected, because the structure showed that the [Fe(CN)61 group was absent from similar to 25% of the structure, invalidating the previously presumed soluble KFe(+3)[Fe(2+)(CN)(6)] structure. This information led to the conclusion that the real structure of Prussian Blue electrochemically synthesized after the stabilization process is Fe(4)[Fe(CN)(6)](3)center dot mH(2)O containing a certain fraction of inserted K(+) ions. In fact, based on an electrogravimetric analysis (Gimenez-Romero et al., J. Phys. Chem. B 2006, 110, 2715 and 19352) complemented by the Fourier maps. it is possible to affirm that the K(+) was part of the water crystalline substructure. Therefore, the interplay mechanism was reexamined considering more precisely the role played by the water crystalline substructure and the K+ alkali metal ion. As a final conclusion, it is proposed that the most precise way to represent the structure of electrochemically synthesized and stabilized hexacyanoferrate materials is Fe(4)(3+) Fe(2+)(CN)(6)](3)center dot[K(h)(+)center dot OH(h)(-)center dot mH(2)O]. The importance of this result is that the widespread use of the terms soluble and insoluble in the electrochemical literature could be reconsidered. Indeed, only one type of structure is insoluble, and that is Fe(4)[Fe(CN)(6)](3)center dot mH(2)O hence, the use of the terms soluble and insoluble is inappropriate from a structural point of view. The result of the presence of the [Fe(CN)61 vacancy a, roup is that the water Substructure cannot be ignored in the ionic interplay mechanism which controls the intercalation and redox process, as was previously confirmed by electrogravimetric analyses (Gimenez-Romero et al., J. Phys. Chem. B 2006, 110, 2715 Garcia-Jareno et al., Electrochim. Acta 1998, 44, 395: Kulesza, Inorg. Chem. 1990, 29, 2395).