Assessing influence in survival data with a cure fraction and covariates

Diagnostic methods have been an important tool in regression analysis to detect anomalies, such as departures from error assumptions and the presence of outliers and influential observations with the fitted models. Assuming censored data, we considered a classical analysis and Bayesian analysis assuming no informative priors for the parameters of the model with a cure fraction. A Bayesian approach was considered by using Markov Chain Monte Carlo Methods with Metropolis-Hasting algorithms steps to obtain the posterior summaries of interest. Some influence methods, such as the local influence, total local influence of an individual, local influence on predictions and generalized leverage were derived, analyzed and discussed in survival data with a cure fraction and covariates. The relevance of the approach was illustrated with a real data set, where it is shown that, by removing the most influential observations, the decision about which model best fits the data is changed.

High-throughput SNP genotyping in the highly heterozygous genome of Eucalyptus: assay success, polymorphism and transferability across species

Background: High-throughput SNP genotyping has become an essential requirement for molecular breeding and population genomics studies in plant species. Large scale SNP developments have been reported for several mainstream crops. A growing interest now exists to expand the speed and resolution of genetic analysis to outbred species with highly heterozygous genomes. When nucleotide diversity is high, a refined diagnosis of the target SNP sequence context is needed to convert queried SNPs into high-quality genotypes using the Golden Gate Genotyping Technology (GGGT). This issue becomes exacerbated when attempting to transfer SNPs across species, a scarcely explored topic in plants, and likely to become significant for population genomics and inter specific breeding applications in less domesticated and less funded plant genera. Results: We have successfully developed the first set of 768 SNPs assayed by the GGGT for the highly heterozygous genome of Eucalyptus from a mixed Sanger/454 database with 1,164,695 ESTs and the preliminary 4.5X draft genome sequence for E. grandis. A systematic assessment of in silico SNP filtering requirements showed that stringent constraints on the SNP surrounding sequences have a significant impact on SNP genotyping performance and polymorphism. SNP assay success was high for the 288 SNPs selected with more rigorous in silico constraints; 93% of them provided high quality genotype calls and 71% of them were polymorphic in a diverse panel of 96 individuals of five different species. SNP reliability was high across nine Eucalyptus species belonging to three sections within subgenus Symphomyrtus and still satisfactory across species of two additional subgenera, although polymorphism declined as phylogenetic distance increased. Conclusions: This study indicates that the GGGT performs well both within and across species of Eucalyptus notwithstanding its nucleotide diversity >= 2%. The development of a much larger array of informative SNPs across multiple Eucalyptus species is feasible, although strongly dependent on having a representative and sufficiently deep collection of sequences from many individuals of each target species. A higher density SNP platform will be instrumental to undertake genome-wide phylogenetic and population genomics studies and to implement molecular breeding by Genomic Selection in Eucalyptus.

A role for mammalian target of rapamycin (mTOR) pathway in non alcoholic steatohepatitis related-cirrhosis

Non-alcoholic fatty liver disease (NAFLD) encompasses the whole spectrum of steatosis, nonalcoholic steatohepatitis (NASH), and NASH-related cirrhosis (NASH/Cir). Although molecular advances have been made in this field, the pathogenesis of NAFLD is not completely understood. The gene expression profiling associated to NASH/Cir was assessed, in an attempt to better characterize the pathways involved in its etiopathogenesis. Methods: In the first step, we used cDNA microarray to evaluate the gene expression profiles in normal liver (n=3) and NASH/Cir samples (n=3) by GeneSifter (TM) analysis to identify differentially expressed genes and biological pathways. Second, tissue microarray was used to determine immunohistochemical expression of phosphorylated mTOR and 4E-BP1 in 11 normal liver samples, 10 NASH/Cir samples and in 37 samples of cirrhosis of other etiologies to further explore the involvement of the mTOR pathway evidenced by the gene expression analysis. Results: 138 and 106 genes were, respectively, up and down regulated in NASH/Cir in comparison to normal liver. Among the 9 pathways identified as significantly modulated in NASH/Cir, the participation of the mTOR pathway was confirmed, since expression of cytoplasmic and membrane phospho-mTOR were higher in NASH/Cir in comparison to cirrhosis of other etiologies and to normal liver. Conclusions: Recent findings have suggested a role for the cellular ""nutrient sensor"" mTOR in NAFLD and the present study corroborates the participation of this pathway in NASH/Cir. Phospho-mTOR evaluation might be of clinical utility as a potential marker for identification of NASH/Cir in cases mistakenly considered as cryptogenic cirrhosis owing to paucity of clinical data.

Disequilibrium Coefficient: A Bayesian Perspective

Hardy-Weinberg Equilibrium (HWE) is an important genetic property that populations should have whenever they are not observing adverse situations as complete lack of panmixia, excess of mutations, excess of selection pressure, etc. HWE for decades has been evaluated; both frequentist and Bayesian methods are in use today. While historically the HWE formula was developed to examine the transmission of alleles in a population from one generation to the next, use of HWE concepts has expanded in human diseases studies to detect genotyping error and disease susceptibility (association); Ryckman and Williams (2008). Most analyses focus on trying to answer the question of whether a population is in HWE. They do not try to quantify how far from the equilibrium the population is. In this paper, we propose the use of a simple disequilibrium coefficient to a locus with two alleles. Based on the posterior density of this disequilibrium coefficient, we show how one can conduct a Bayesian analysis to verify how far from HWE a population is. There are other coefficients introduced in the literature and the advantage of the one introduced in this paper is the fact that, just like the standard correlation coefficients, its range is bounded and it is symmetric around zero (equilibrium) when comparing the positive and the negative values. To test the hypothesis of equilibrium, we use a simple Bayesian significance test, the Full Bayesian Significance Test (FBST); see Pereira, Stern andWechsler (2008) for a complete review. The disequilibrium coefficient proposed provides an easy and efficient way to make the analyses, especially if one uses Bayesian statistics. A routine in R programs (R Development Core Team, 2009) that implements the calculations is provided for the readers.

Prognostic index expression of cyclin-D1, cerbB-2 and VEGF: metastases vs corresponding primary cancers and metastatic vs non-metastatic adenocarcinomas

The prognostic relevance of different molecular markers in lung cancer is a crucial issue still worth investigating, and the specimens collected and analyzed represent a valuable source of material. Cyclin-D1, c-erbB-2 and vascular endothelial growth factor (VEGF) have shown to be promising as prognosticators in human cancer. In this study, we sought to examine the importance of Cyclin-D1, c-erbB-2 and VEGF, and to study the quantitative relationship among these factors and disease progression in metastases vs corresponding primary cancer, and metastatic vs non metastatic cancers. Material and Methods: We used immunohistochemistry and morphometric analysis to evaluate the amount of tumour staining for Cyclin-D1, c-erbB-2 and VEGF in 52 patients with surgically excised ademocarcinoma of the lung, and the outcome for our study was survival time until death from hematogenic metastases. Results: Metastasis presented lower c-erbB-2 expression than corresponding primary cancers (p=0.02). Cyclin-D1 and VEGF expression were also lower in metastases than in corresponding primary cancers, but this difference did not achieve statistical significance. Non-metastatic cancers also presented significantly lower Cyclin-D1 and c-erbB-2 expression than metastatic cancers (p<0.01 and p<0.01, respectively). Equally significant was the difference between higher c-erbB-2 expression by metastatic cancers compared to non-metastatic cancers (p=0.02). Considering survival in Kaplan-Maier analysis, Cyclin-D1 (p=0.04), c-erbB-2 (p=0.04) and VEGF (p<0.01) were important predictors of survival in metastatic cancers.

Chaotic model and memory in single calcium-activated potassium channel kinetics

Ion channels are pores formed by proteins and responsible for carrying ion fluxes through cellular membranes. The ion channels can assume conformational states thereby controlling ion flow. Physically, the conformational transitions from one state to another are associated with energy barriers between them and are dependent on stimulus, such as, electrical field, ligands, second messengers, etc. Several models have been proposed to describe the kinetics of ion channels. The classical Markovian model assumes that a future transition is independent of the time that the ion channel stayed in a previous state. Others models as the fractal and the chaotic assume that the rate of transitions between the states depend on the time that the ionic channel stayed in a previous state. For the calcium activated potassium channels of Leydig cells the R/S Hurst analysis has indicated that the channels are long-term correlated with a Hurst coefficient H around 0.7, showing a persistent memory in this kinetic. Here, we applied the R/S analysis to the opening and closing dwell time series obtained from simulated data from a chaotic model proposed by L. Liebovitch and T. Toth [J. Theor. Biol. 148, 243 (1991)] and we show that this chaotic model or any model that treats the set of channel openings and closings as independent events is inadequate to describe the long-term correlation (memory) already described for the experimental data. (C) 2008 American Institute of Physics.

Energy gain induced by boundary crisis

We consider a nonlinear system and show the unexpected and surprising result that, even for high dissipation, the mean energy of a particle can attain higher values than when there is no dissipation in the system. We reconsider the time-dependent annular billiard in the presence of inelastic collisions with the boundaries. For some magnitudes of dissipation, we observe the phenomenon of boundary crisis, which drives the particles to an asymptotic attractive fixed point located at a value of energy that is higher than the mean energy of the nondissipative case and so much higher than the mean energy just before the crisis. We should emphasize that the unexpected results presented here reveal the importance of a nonlinear dynamics analysis to explain the paradoxical strategy of introducing dissipation in the system in order to gain energy.

Constraints on the infrared behavior of the ghost propagator in Yang-Mills theories

We present rigorous upper and lower bounds for the momentum-space ghost propagator G(p) of Yang-Mills theories in terms of the smallest nonzero eigenvalue (and of the corresponding eigenvector) of the Faddeev-Popov matrix. We apply our analysis to data from simulations of SU(2) lattice gauge theory in Landau gauge, using the largest lattice sizes to date. Our results suggest that, in three and in four space-time dimensions, the Landau gauge ghost propagator is not enhanced as compared to its tree-level behavior. This is also seen in plots and fits of the ghost dressing function. In the two-dimensional case, on the other hand, we find that G(p) diverges as p(-2-2 kappa) with kappa approximate to 0.15, in agreement with A. Maas, Phys. Rev. D 75, 116004 (2007). We note that our discussion is general, although we make an application only to pure gauge theory in Landau gauge. Our simulations have been performed on the IBM supercomputer at the University of Sao Paulo.

Greenhouse gas emissions and energy balance of palm oil biofuel

The search for alternatives to fossil fuels is boosting interest in biodiesel production. Among the crops used to produce biodiesel, palm trees stand out due to their high productivity and positive energy balance. This work assesses life cycle emissions and the energy balance of biodiesel production from palm oil in Brazil. The results are compared through a meta-analysis to previous published studies: Wood and Corley (1991) [Wood BJ, Corley RH. The energy balance of oil palm cultivation. In: PORIM intl. palm oil conference agriculture; 1991.], Malaysia; Yusoff and Hansen (2005) [Yusoff S. Hansen SB. Feasibility study of performing an life cycle assessment on crude palm oil production in Malaysia. International Journal of Life Cycle Assessment 2007;12:50-8], Malaysia; Angarita et al. (2009) [Angarita EE, Lora EE, Costa RE, Torres EA. The energy balance in the palm oil-derived methyl ester (PME) life cycle for the cases in Brazil and Colombia. Renewable Energy 2009;34:2905-13], Colombia; Pleanjai and Gheewala (2009) [Pleanjai S. Gheewala SH. Full chain energy analysis of biodiesel production from palm oil in Thailand. Applied Energy 2009;86:S209-14], Thailand; and Yee et al. (2009) [Yee KF, Tan KT, Abdullah AZ, Lee la. Life cycle assessment of palm biodiesel: revealing facts and benefits for sustainability. Applied Energy 2009;86:S189-96], Malaysia. In our study, data for the agricultural phase, transport, and energy content of the products and co-products were obtained from previous assessments done in Brazil. The energy intensities and greenhouse gas emission factors were obtained from the Simapro 7.1.8. software and other authors. These factors were applied to the inputs and outputs listed in the selected studies to render them comparable. The energy balance for our study was 1:5.37. In comparison the range for the other studies is between 1:3.40 and 1:7.78. Life cycle emissions determined in our assessment resulted in 1437 kg CO(2)e/ha, while our analysis based on the information provided by other authors resulted in 2406 kg CO(2)e/ha, on average. The Angarita et al. (2009) [Angarita EE, Lora EE, Costa RE, Torres EA. The energy balance in the palm oil-derived methyl ester (PME) life cycle for the cases in Brazil and Colombia. Renewable Energy 2009:34:2905-13] study does not report emissions. When compared to diesel on a energy basis, avoided emissions due to the use of biodiesel account for 80 g CO(2)e/MJ. Thus, avoided life Cycle emissions associated with the use of biodiesel yield a net reduction of greenhouse gas emissions. We also assessed the carbon balance between a palm tree plantation, including displaced emissions from diesel, and a natural ecosystem. Considering the carbon balance outcome plus life cycle emissions the payback time for a tropical forest is 39 years. The result published by Gibbs et al. (2008) [Gibbs HK, Johnston M, Foley JA, Holloway T, Monfreda C, Ramankutty N, et al., Carbon payback times for crop-based biofuel expansion in the tropics: the effects of changing yield and technology. Environmental Research Letters 2008;3:10], which ignores life cycle emissions, determined a payback range for biodiesel production between 30 and 120 years. Crown Copyright (C) 2010 Published by Elsevier Ltd. All rights reserved.

Failure assessments of corroded pipelines with axial defects using stress-based criteria: Numerical studies and verification analyses

Conventional procedures used to assess the integrity of corroded piping systems with axial defects generally employ simplified failure criteria based upon a plastic collapse failure mechanism incorporating the tensile properties of the pipe material. These methods establish acceptance criteria for defects based on limited experimental data for low strength structural steels which do not necessarily address specific requirements for the high grade steels currently used. For these cases, failure assessments may be overly conservative or provide significant scatter in their predictions, which lead to unnecessary repair or replacement of in-service pipelines. Motivated by these observations, this study examines the applicability of a stress-based criterion based upon plastic instability analysis to predict the failure pressure of corroded pipelines with axial defects. A central focus is to gain additional insight into effects of defect geometry and material properties on the attainment of a local limit load to support the development of stress-based burst strength criteria. The work provides an extensive body of results which lend further support to adopt failure criteria for corroded pipelines based upon ligament instability analyses. A verification study conducted on burst testing of large-diameter pipe specimens with different defect length shows the effectiveness of a stress-based criterion using local ligament instability in burst pressure predictions, even though the adopted burst criterion exhibits a potential dependence on defect geometry and possibly on material`s strain hardening capacity. Overall, the results presented here suggests that use of stress-based criteria based upon plastic instability analysis of the defect ligament is a valid engineering tool for integrity assessments of pipelines with axial corroded defects. (C) 2008 Elsevier Ltd. All rights reserved.

Outsourcing decision-making aspects considered by IT departments in Brazilian companies

The study of Information Technology (IT) outsourcing is relevant because companies are outsourcing their activities more than ever. An important IT outsourcing research area is the decision-making process. In other words, the comprehension of how companies decide about outsourcing their IT operations is relevant from research point of view. Therefore, the objective of this study is to understand the decision-making process used by Brazilian companies when outsourcing their IT operations. An analysis of the literature that refers to this subject showed that six aspects are usually considered by companies on the evaluation of IT outsourcing service alternatives. This research verified how these six aspects are considered by Brazilian companies on IT outsourcing decisions. The survey showed that Brazilian companies consider all the six aspects, but each of them has a different level of importance. The research also grouped the aspects according to their level of importance and interdependency, using factorial analysis to understand the logic behind IT outsourcing decision process. (C) 2009 Elsevier B.V. All rights reserved.

Electro-optically sensitive diamond-like carbon thin films deposited by reactive magnetron sputtering for electronic device applications

The goal of this work is to study and relate electrical and optical properties of diamond-like carbon (DLC) thin films for applications in electronic devices. DLC films were deposited in a reactive RF magnetron sputtering system on p-type silicon and glass substrates. The target was a 99.9999% pure, 6 in. diameter graphite plate and methane was used as processing gas. Eight DLC films were produced for each substrate, varying deposition time, the reactor pressure between 5 mTorr and 10 mTorr while the RF power was applied at 13.56 MHz and varied between 100, 150, 200 and 250W. After deposition, the films were analyzed by I-V and C-V measurements (Cheng et al. (2004) [1]) in order to determine the electric resistivity, photo-current response and dielectric constant, optical transmittance, used to find the optical gap by the Tauc method; and by photoluminescence analysis to determine the photoemission and confirm the optical band gap. These characteristics are compared and the influence of the deposition parameters is discussed. (C) 2011 Published by Elsevier B.V.

Characterization of fungal soil communities by F-RISA and arbuscular mycorrhizal fungi from Araucaria angustifolia forest soils after replanting and wildfire disturbances

The Fungal Ribosomal Intergenic Spacer Analysis (F-RISA) was used to characterize soil fungal communities from three ecosystems of Araucaria angustifolia from Brazil: a native forest and two replanted forest ecosystems, one of them with a past history of wildfire. The arbuscular mycorrhizal fungi (AMF) infection was evaluated in Araucaria roots of 18-month-old axenic plants previously inoculated with soils collected from those areas in a greenhouse experiment. The principal component analysis of F-RISA profiles showed different soil fungal community between the three studied areas. Sixty three percent of F-RISA fragments amplified in the soil and the substrate samples presented lengths between 500 and 700 bp. The number of Operational Taxonomic Units (OTUs) was 34 for soil and 38 for substrate, however, more fragments were detected in soil (214) than in substrate (163). An in silico F-RISA analysis to compare our data with ITS1-5.8S-ITS2 sequences from NCBI database showed the presence of Ascomycota, Basidiomycota and Glomeromycota among the soil and substrate fungal communities. AMF infection was higher in plants inoculated with soil from the native forest and the replanted forest with wildfire, both presenting similar chemical characteristics but with different disturbance levels. These results indicate that soil chemical composition may influence the soil fungal community structures rather than the anthropogenic or fire disturbances.

Estimating Amino Acid Requirement of Brazilian Freshwater Fish from Muscle Amino Acid Profile

Information on nutritional requirement of some Brazilian farmed fish species, especially essential amino acids (EAA) requirements, is scarce. The estimation of amino acids requirements based on amino acid composition of fish is a fast and reliable alternative. Matrinxa, Brycon amazonicus, and curimbata, Prochilodus lineatus, are two important Brazilian fish with potential for aquaculture. The objective of the present study was to estimate amino acid requirements of these species and analyze similarities among amino acid composition of different fish species by cluster analysis. To estimate amino acid requirement, the following formula was used: amino acid requirement = [(amount of an individual amino acid in fish muscle tissue) x (average totalEAA requirement among channel catfish, Ictalurus punctatus, Nile tilapia, Oreochromis niloticus, and common carp, Cyprinus carpio)]/(average fish muscle totalEAA). Most values found lie within the range of requirements determined for other omnivorous fish species, in exception of leucine requirement estimated for both species, and arginine requirement estimated for matrinxa alone. Rather than writing off the need for regular dose-response assays under the ideal protein concept to determine EAA requirements of curimbata and matrinxa, results set solid base for the study of tropical species dietary amino acids requirements.

A note on the multiple partners assignment game

In the assignment game of Shapley and Shubik [Shapley, L.S., Shubik, M., 1972. The assignment game. I. The core, International journal of Game Theory 1, 11-130] agents are allowed to form one partnership at most. That paper proves that, in the context of firms and workers, given two stable payoffs for the firms there is a stable payoff which gives each firm the larger of the two amounts and also one which gives each of them the smaller amount. Analogous result applies to the workers. Sotomayor [Sotomayor, M., 1992. The multiple partners game. In: Majumdar, M. (Ed.), Dynamics and Equilibrium: Essays in Honor to D. Gale. Mcmillian, pp. 322-336] extends this analysis to the case where both types of agents may form more than one partnership and an agent`s payoff is multi-dimensional. Instead, this note concentrates in the total payoff of the agents. It is then proved the rather unexpected result that again the maximum of any pair of stable payoffs for the firms is stable but the minimum need not be, even if we restrict the multiplicity of partnerships to one of the sides. (C) 2009 Elsevier B.V. All rights reserved.