Regression models for grouped survival data: Estimation and sensitivity analysis

In this study, regression models are evaluated for grouped survival data when the effect of censoring time is considered in the model and the regression structure is modeled through four link functions. The methodology for grouped survival data is based on life tables, and the times are grouped in k intervals so that ties are eliminated. Thus, the data modeling is performed by considering the discrete models of lifetime regression. The model parameters are estimated by using the maximum likelihood and jackknife methods. To detect influential observations in the proposed models, diagnostic measures based on case deletion, which are denominated global influence, and influence measures based on small perturbations in the data or in the model, referred to as local influence, are used. In addition to those measures, the local influence and the total influential estimate are also employed. Various simulation studies are performed and compared to the performance of the four link functions of the regression models for grouped survival data for different parameter settings, sample sizes and numbers of intervals. Finally, a data set is analyzed by using the proposed regression models. (C) 2010 Elsevier B.V. All rights reserved.

The exponentiated generalized gamma distribution with application to lifetime data

A four-parameter extension of the generalized gamma distribution capable of modelling a bathtub-shaped hazard rate function is defined and studied. The beauty and importance of this distribution lies in its ability to model monotone and non-monotone failure rate functions, which are quite common in lifetime data analysis and reliability. The new distribution has a number of well-known lifetime special sub-models, such as the exponentiated Weibull, exponentiated generalized half-normal, exponentiated gamma and generalized Rayleigh, among others. We derive two infinite sum representations for its moments. We calculate the density of the order statistics and two expansions for their moments. The method of maximum likelihood is used for estimating the model parameters and the observed information matrix is obtained. Finally, a real data set from the medical area is analysed.

Double generalized linear model for tissue culture proportion data: a Bayesian perspective

Joint generalized linear models and double generalized linear models (DGLMs) were designed to model outcomes for which the variability can be explained using factors and/or covariates. When such factors operate, the usual normal regression models, which inherently exhibit constant variance, will under-represent variation in the data and hence may lead to erroneous inferences. For count and proportion data, such noise factors can generate a so-called overdispersion effect, and the use of binomial and Poisson models underestimates the variability and, consequently, incorrectly indicate significant effects. In this manuscript, we propose a DGLM from a Bayesian perspective, focusing on the case of proportion data, where the overdispersion can be modeled using a random effect that depends on some noise factors. The posterior joint density function was sampled using Monte Carlo Markov Chain algorithms, allowing inferences over the model parameters. An application to a data set on apple tissue culture is presented, for which it is shown that the Bayesian approach is quite feasible, even when limited prior information is available, thereby generating valuable insight for the researcher about its experimental results.

Evaluation of FAO Penman-Monteith and alternative methods for estimating reference evapotranspiration with missing data in Southern Ontario, Canada

Grass reference evapotranspiration (ETo) is an important agrometeorological parameter for climatological and hydrological studies, as well as for irrigation planning and management. There are several methods to estimate ETo, but their performance in different environments is diverse, since all of them have some empirical background. The FAO Penman-Monteith (FAD PM) method has been considered as a universal standard to estimate ETo for more than a decade. This method considers many parameters related to the evapotranspiration process: net radiation (Rn), air temperature (7), vapor pressure deficit (Delta e), and wind speed (U); and has presented very good results when compared to data from lysimeters Populated with short grass or alfalfa. In some conditions, the use of the FAO PM method is restricted by the lack of input variables. In these cases, when data are missing, the option is to calculate ETo by the FAD PM method using estimated input variables, as recommended by FAD Irrigation and Drainage Paper 56. Based on that, the objective of this study was to evaluate the performance of the FAO PM method to estimate ETo when Rn, Delta e, and U data are missing, in Southern Ontario, Canada. Other alternative methods were also tested for the region: Priestley-Taylor, Hargreaves, and Thornthwaite. Data from 12 locations across Southern Ontario, Canada, were used to compare ETo estimated by the FAD PM method with a complete data set and with missing data. The alternative ETo equations were also tested and calibrated for each location. When relative humidity (RH) and U data were missing, the FAD PM method was still a very good option for estimating ETo for Southern Ontario, with RMSE smaller than 0.53 mm day(-1). For these cases, U data were replaced by the normal values for the region and Delta e was estimated from temperature data. The Priestley-Taylor method was also a good option for estimating ETo when U and Delta e data were missing, mainly when calibrated locally (RMSE = 0.40 mm day(-1)). When Rn was missing, the FAD PM method was not good enough for estimating ETo, with RMSE increasing to 0.79 mm day(-1). When only T data were available, adjusted Hargreaves and modified Thornthwaite methods were better options to estimate ETo than the FAO) PM method, since RMSEs from these methods, respectively 0.79 and 0.83 mm day(-1), were significantly smaller than that obtained by FAO PM (RMSE = 1.12 mm day(-1). (C) 2009 Elsevier B.V. All rights reserved.

Spatio-temporal modeling of agricultural yield data with an application to pricing crop insurance contracts

This article presents a statistical model of agricultural yield data based on a set of hierarchical Bayesian models that allows joint modeling of temporal and spatial autocorrelation. This method captures a comprehensive range of the various uncertainties involved in predicting crop insurance premium rates as opposed to the more traditional ad hoc, two-stage methods that are typically based on independent estimation and prediction. A panel data set of county-average yield data was analyzed for 290 counties in the State of Parana (Brazil) for the period of 1990 through 2002. Posterior predictive criteria are used to evaluate different model specifications. This article provides substantial improvements in the statistical and actuarial methods often applied to the calculation of insurance premium rates. These improvements are especially relevant to situations where data are limited.

Progeny evaluation for resistance to Phaeosphaeria leaf spot in tropical maize

Maize breeding programmes in Brazil and elsewhere seek reliable methods to identify genotypes resistant to Phaeosphaeria leaf spot. The area under the disease progress curve (AUDPC) is an accurate method to evaluate the severity of foliar diseases. However, at least three data points are required to calculate the AUDPC, which is unfeasible when there are thousands of genotypes to be assessed. The aim of this work was to estimate the heritability of disease resistance, evaluate disease severity at different times using a nine-point scale in comparison to the AUDPC, and establish the most suitable phenological period for disease assessment. A repeated experiment was conducted in a 11 x 11 lattice experimental design with three replications. Disease assessments were carried out at flowering, 15 and 30 days post-anthesis for the parental lines DS95, DAS21, the F1 generation and 118 F2:3 progenies. Then, the AUDPC was obtained and results compared with the single-point evaluations used to calculate it. Individual and joint analyses of variance were conducted to obtain heritabiliy estimates. The assessments performed after the flowering stage gave higher estimates of heritability and correlation with AUDPC. We concluded that one assessment between the 15th and 30th day after flowering could provide enough information to distinguish maize genotypes for their resistance to Phaeosphaeria leaf spot under tropical conditions.

Quantitative Trait Loci Mapping and The Genetic Basis of Heterosis in Maize and Rice

Despite its importance to agriculture, the genetic basis of heterosis is still not well understood. The main competing hypotheses include dominance, overdominance, and epistasis. NC design III is an experimental design that. has been used for estimating the average degree of dominance of quantitative trait 106 (QTL) and also for studying heterosis. In this study, we first develop a multiple-interval mapping (MIM) model for design III that provides a platform to estimate the number, genomic positions, augmented additive and dominance effects, and epistatic interactions of QTL. The model can be used for parents with any generation of selling. We apply the method to two data sets, one for maize and one for rice. Our results show that heterosis in maize is mainly due to dominant gene action, although overdominance of individual QTL could not completely be ruled out due to the mapping resolution and limitations of NC design III. For rice, the estimated QTL dominant effects could not explain the observed heterosis. There is evidence that additive X additive epistatic effects of QTL could be the main cause for the heterosis in rice. The difference in the genetic basis of heterosis seems to be related to open or self pollination of the two species. The MIM model for NC design III is implemented in Windows QTL Cartographer, a freely distributed software.

Nutritional Aspects of Second Generation Soy Foods

Samples of 15 second generation soy-based products (n = 3), commercially available, were analyzed for their protein and isoflavone contents and in vitro antioxidant activity, by means of the Folin-Ciocalteu reducing ability, DPPH radical scavenging capacity, and oxygen radical absorbance capacity. Isoflavone identification and quantification were performed by high-performance liquid chromatography. Products containing soy and/or soy-based ingredients represent important sources of protein in addition to the low fat amounts. However, a large variation in isoflavone content and in vitro antioxidant capacity was observed. The isoflavone content varied from 2.4 to 18.1 mg/100 g (FW), and soy kibe and soy sausage presented the highest amounts. Chocolate had the highest antioxidant capacity, but this fact was probably associated with the addition of cocoa liquor, a well-known source of polyphenolics. This study showed that the soy-based foods do not present a significant content of isoflavones when compared with the grain, and their in vitro antioxidant capacity is not related with these compounds but rather to the presence of other phenolics and synthetic antioxidants, such as sodium erythorbate. However, they may represent alternative sources and provide soy protein, isoflavones, and vegetable fat for those who are not ready to eat traditional soy foods.

Effect of Brazil nut supplementation on the blood levels of selenium and glutathione peroxidase in hemodialysis patients

Objective: In patients who have undergone hemodialysis, large amounts of reactive oxygen species (ROS) are produced and, at higher concentrations, ROS are thought to be involved in the pathogenesis of cardiovascular disease. It has been proposed that selenium (Se) may exert an anti-atherogenic influence by reducing oxidative stress. The richest known food source of selenium is the Brazil nut (Bertholletia excelsa, family Lecythidaceae), found in the Amazon region. We evaluated the effect of Brazil nut supplementation on blood levels of Se and glutathione peroxidase (GSH-Px) activity in patients on hemodialysis. Methods: A total of 81 patients on hemodialysis (52.0 +/- 15.2 y old, average time on dialysis 82.3 +/- 91.4 mo, body mass index 24.9 +/- 4.4 kg/m(2)) from the RenalCor and RenalVida Clinics in Rio de Janeiro, Brazil, were studied. All patients received one nut (around 5 g, averaging 58.1 mu g Se/g) a day for 3 mo. The Se concentrations in the nuts and in plasma and erythrocytes were determined by atomic absorption spectrophotometry with hydride generation (Hitachi, Z-500). GSH-Px levels were measured using Randox commercial kits. Results: Plasma Se (18.8 +/- 17.4 mu g/L) and erythrocyte (72.4 +/- 37.9 mg/L) levels were below the normal, range before nut supplementation. After supplementation, the plasma level increased to 104.0 +/- 65.0 mu g/L and erythrocytes to 244.1 +/- 119.5 mg/L (P<0.0001). The activity of GSH-Px also increased after supplementation, from 46.6 +/- 14.9 to 55.9 +/- 23.6 U/g of hemoglobin (P<0.0001). Before supplementation, 11% of patients had GSH-Px activity below the normal range (27.5-73.6 U/g of hemoglobin). After supplementation, all patients showed GSH-Px activity within the normal range. Conclusion: The data revealed that the investigated patients presented Se deficiency and that the consumption of only one Brazil nut a day (5 g) during 3 mo was effective to increase the Se concentration and GSH-Px activity in these patients, thus improving their antioxidant status. (C) 2010 Elsevier Inc. All rights reserved.

Method development and optimization for the determination of selenium in bean and soil samples using hydride generation electrothermal atomic absorption spectrometry

The present investigation is the first part of an initiative to prepare a regional map of the natural abundance of selenium in various areas of Brazil, based on the analysis of bean and soil samples. Continuous-flow hydride generation electrothermal atomic absorption spectrometry (HG-ET AAS) with in situ trapping on an iridium-coated graphite tube has been chosen because of the high sensitivity and relative simplicity. The microwave-assisted acid digestion for bean and soil samples was tested for complete recovery of inorganic and organic selenium compounds (selenomethionine). The reduction of Se(VI) to Se(IV) was optimized in order to guarantee that there is no back-oxidation, which is of importance when digested samples are not analyzed immediately after the reduction step. The limits of detection and quantification of the method were 30 ng L(-1) Se and 101 ng L(-1) Se, respectively, corresponding to about 3 ng g(-1) and 10 ng g(-1), respectively, in the solid samples, considering a typical dilution factor of 100 for the digestion process. The results obtained for two certified food reference materials (CRM), soybean and rice, and for a soil and sediment CRM confirmed the validity of the investigated method. The selenium content found in a number of selected bean samples varied between 5.5 +/- 0.4 ng g(-1) and 1726 +/- 55 ng g(-1), and that in soil samples varied between 113 +/- 6.5 ng g(-1) and 1692 +/- 21 ng g(-1). (C) 2011 Elsevier B.V. All rights reserved.

Genetic population data of 12 STR loci of the PowerPlex (R) Y system in the state of Sao Paulo population (Southeast of Brazil)

Allele frequency distributions and population data for 12 Y-chromosomal short tandem repeats (STRs) included in the PowerPlex (R) Y Systems (Promega) were obtained for a sample of 200 healthy unrelated males living in S (a) over tildeo Paulo State (Southeast of Brazil). A total of 192 haplotypes were identified, of which 184 were unique and 8 were found in 2 individuals. The average gene diversity of the 12 Y-STR was 0.6746 and the haplotype diversity was 0.9996. Pairwise analysis confirmed that our population is more similar with the Italy, North Portugal and Spain, being more distant of the Japan. (c) 2007 Elsevier Ireland Ltd. All rights reserved.

Effect of Simvastatin on the Generation of Autoantibodies Against Oxidized LDL and Progression of Atherosclerosis in Rabbits

Oxidized Low-Density Lipoproteins (oxLDL) and autoantibodies against oxLDL are important in the development of atherosclerotic lesions. Statins are efficacious in the control of dyslipidemia and prevention of atherosclerosis; however, many questions concerning the mechanism of action of such drugs remain unknown. This work investigated the effect of simvastatin on generation of autoantibodies against oxLDL and development of atherosclerosis in rabbits. The animals were divided into three groups: control, hypercholesterolemic, and hypercholesterolemic simvastatin (3.0 mg simvastatin/ kg body weight). Concentrations of autoantibodies against oxLDL were determined on days 0,30 and 60 of the experiment and the atherosclerotic lesions were evaluated at the end of the study. Simvastatin reduced intimal proliferation in the thoracic region, prevented arterial calcification and inhibited the generation of autoantibodies against oxLDL. In conclusion, daily administration of simvastatin slows down atherosclerotic lesion development in rabbits with induced hypercholesterolemia and inhibition on generation of autoantibodies against oxLDL contributes to the cardioprotective effect observed.

Generation and Analysis of Interaction Energy Maps of p-Substituted Benzoic Acid [(5-Nitro-thiophen-2-yl)-methylene]-Hydrazides Active against MRSA

In this preliminary study eighteen p-substituted benzoic acid [(5-nitro-thiophen-2-yl)-methylene]-hydrazides with antimicrobial activity were evaluated against multidrug-resistant Staphylococcus aureus, correlating the three-dimensional characteristics of the ligands with their respective bioactivities. The computer programs Sybyl and CORINA were used, respectively, for the design and three-dimensional conversion of the ligands. Molecular interaction fields were calculated using GRID program. Calculations using Volsurf resulted in a statistically consistent model with 48 structural descriptors showing that hydrophobicity is a fundamental property in the analyzed biological response.

Brazilian Network of Food Data Systems and LATINFOODS Regional Technical Compilation Committee: Food composition activities (2006-2009)

The Brazilian Network of Food Data Systems (BRASILFOODS) has been keeping the Brazilian Food Composition Database-USP (TBCA-USP) (http://www.fcf.usp.br/tabela) since 1998. Besides the constant compilation, analysis and update work in the database, the network tries to innovate through the introduction of food information that may contribute to decrease the risk for non-transmissible chronic diseases, such as the profile of carbohydrates and flavonoids in foods. In 2008, data on carbohydrates, individually analyzed, of 112 foods, and 41 data related to the glycemic response produced by foods widely consumed in the country were included in the TBCA-USP. Data (773) about the different flavonoid subclasses of 197 Brazilian foods were compiled and the quality of each data was evaluated according to the USDAs data quality evaluation system. In 2007, BRASILFOODS/USP and INFOODS/FAO organized the 7th International Food Data Conference ""Food Composition and Biodiversity"". This conference was a unique opportunity for interaction between renowned researchers and participants from several countries and it allowed the discussion of aspects that may improve the food composition area. During the period, the LATINFOODS Regional Technical Compilation Committee and BRASILFOODS disseminated to Latin America the Form and Manual for Data Compilation, version 2009, ministered a Food Composition Data Compilation course and developed many activities related to data production and compilation. (C) 2010 Elsevier Inc. All rights reserved.

Expression patterns of cell wall-modifying genes from banana during fruit ripening and in relationship with finger drop

Few molecular studies have been devoted to the finger drop process that occurs during banana fruit ripening. Recent studies revealed the involvement of changes in the properties of cell wall polysaccharides in the pedicel rupture area. In this study, the expression of cell-wall modifying genes was monitored in peel tissue during post-harvest ripening of Cavendish banana fruit, at median area (control zone) and compared with that in the pedicel rupture area (drop zone). To this end, three pectin methylesterase (PME) and seven xyloglucan endotransglycosylase/hydrolase (XTH) genes were isolated. The accumulation of their mRNAs and those of polygalaturonase, expansin, and pectate lyase genes already isolated from banana were examined. During post-harvest ripening, transcripts of all genes were detected in both zones, but accumulated differentially. MaPME1, MaPG1, and MaXTH4 mRNA levels did not change in either zone. Levels of MaPME3 and MaPG3 mRNAs increased greatly only in the control zone and at the late ripening stages. For other genes, the main molecular changes occurred 1-4 d after ripening induction. MaPME2, MaPEL1, MaPEL2, MaPG4, MaXTH6, MaXTH8, MaXTH9, MaEXP1, MaEXP4, and MaEXP5 accumulated highly in the drop zone, contrary to MaXTH3 and MaXTH5, and MaEXP2 throughout ripening. For MaPG2, MaXET1, and MaXET2 genes, high accumulation in the drop zone was transient. The transcriptional data obtained from all genes examined suggested that finger drop and peel softening involved similar mechanisms. These findings also led to the proposal of a sequence of molecular events leading to finger drop and to suggest some candidates.