Modeling the Potential Geographic Distribution of Black Pepper (Piper nigrum) in Asia Using GIS Tools

Known as the "king of spices", black pepper (Piper nigrum), a perennial crop of the tropics, is economically the most important and the most widely used spice crop in the world. To understand its suitable bioclimatic distribution, maximum entropy based on ecological niche modeling was used to model the bioclimatic niches of the species in its Asian range. Based on known occurrences, bioclimatic areas with higher probabilities are mainly located in the eastern and western coasts of the Indian Peninsula, the east of Sumatra Island, some areas in the Malay Archipelago, and the southeast coastal areas of China. Some undocumented places were also predicted as suitable areas. According to the jackknife procedure, the minimum temperature of the coldest month, the mean monthly temperature range, and the precipitation of the wettest month were identified as highly effective factors in the distribution of black pepper and could possibly account for the crop's distribution pattern. Such climatic requirements inhibited this species from dispersing and gaining a larger geographical range.

Characterization of aflatoxigenic and non-aflatoxigenic strains of Aspergillus section Flavi isolated from corn grains of different geographic origins in Brazil

Aflatoxins can cause great economic losses and serious risks to humans and animals health. The largest aflatoxin producers belong to Aspergillus section Flavi and can occur naturally in food commodities. Studies showed that molecular tools as well as the type of sclerotia produced by the strains could be helpful for identification of Aspergillus species and could be correlated with levels of toxin production. The purpose of this work was to characterize the genetic diversity using AFLP technique, the type of sclerotia and the ability of aflatoxin production by isolated strains from corn of different origins in Brazil, and to verify whether qPCR based on aflR and aflP genes is appropriate for estimating the level of aflatoxin production. All the 75 strains were classified as A. flavus and the AFLP technique showed a wide intraspecific variability within them. Regarding sclerotia production, 34% were classified as S and 66% as L type. Among the aflatoxin-producers, 52.8% produced aflatoxin B-1, while 47.2% aflatoxins B-1 and B-2. Statistical analysis showed no correlation between sclerotia production and aflatoxigenicty, and no correlation between the phylogenetic clusters and aflatoxin production. Concerning the relative expression of aflR and aflP, Pearson's correlation test demonstrated low positive correlation between the expression of the aflR and aflP genes and the production of AFB(1) and AFB(2), but showed high positive correlation between aflR and aflP expression. In contrast to the other reference strains, A. oryzae ATCC 7282 showed no amplification of aflR and aflP. The results highlight the need for detection of reliable and reproducible markers with a high positive correlation with aflatoxin production.

Bartlett corrections in Birnbaum-Saunders nonlinear regression models

Lemonte and Cordeiro [Birnbaum-Saunders nonlinear regression models, Comput. Stat. Data Anal. 53 (2009), pp. 4441-4452] introduced a class of Birnbaum-Saunders (BS) nonlinear regression models potentially useful in lifetime data analysis. We give a general matrix Bartlett correction formula to improve the likelihood ratio (LR) tests in these models. The formula is simple enough to be used analytically to obtain several closed-form expressions in special cases. Our results generalize those in Lemonte et al. [Improved likelihood inference in Birnbaum-Saunders regressions, Comput. Stat. DataAnal. 54 (2010), pp. 1307-1316], which hold only for the BS linear regression models. We consider Monte Carlo simulations to show that the corrected tests work better than the usual LR tests.

Classification and regression tree (CART) model to predict pulmonary tuberculosis in hospitalized patients

Background: Tuberculosis (TB) remains a public health issue worldwide. The lack of specific clinical symptoms to diagnose TB makes the correct decision to admit patients to respiratory isolation a difficult task for the clinician. Isolation of patients without the disease is common and increases health costs. Decision models for the diagnosis of TB in patients attending hospitals can increase the quality of care and decrease costs, without the risk of hospital transmission. We present a predictive model for predicting pulmonary TB in hospitalized patients in a high prevalence area in order to contribute to a more rational use of isolation rooms without increasing the risk of transmission. Methods: Cross sectional study of patients admitted to CFFH from March 2003 to December 2004. A classification and regression tree (CART) model was generated and validated. The area under the ROC curve (AUC), sensitivity, specificity, positive and negative predictive values were used to evaluate the performance of model. Validation of the model was performed with a different sample of patients admitted to the same hospital from January to December 2005. Results: We studied 290 patients admitted with clinical suspicion of TB. Diagnosis was confirmed in 26.5% of them. Pulmonary TB was present in 83.7% of the patients with TB (62.3% with positive sputum smear) and HIV/AIDS was present in 56.9% of patients. The validated CART model showed sensitivity, specificity, positive predictive value and negative predictive value of 60.00%, 76.16%, 33.33%, and 90.55%, respectively. The AUC was 79.70%. Conclusions: The CART model developed for these hospitalized patients with clinical suspicion of TB had fair to good predictive performance for pulmonary TB. The most important variable for prediction of TB diagnosis was chest radiograph results. Prospective validation is still necessary, but our model offer an alternative for decision making in whether to isolate patients with clinical suspicion of TB in tertiary health facilities in countries with limited resources.

Response plateau models fitting via isotonic regression

Within the nutritional context, the supplementation of microminerals in bird food is often made in quantities exceeding those required in the attempt to ensure the proper performance of the animals. The experiments of type dosage x response are very common in the determination of levels of nutrients in optimal food balance and include the use of regression models to achieve this objective. Nevertheless, the regression analysis routine, generally, uses a priori information about a possible relationship between the response variable. The isotonic regression is a method of estimation by least squares that generates estimates which preserves data ordering. In the theory of isotonic regression this information is essential and it is expected to increase fitting efficiency. The objective of this work was to use an isotonic regression methodology, as an alternative way of analyzing data of Zn deposition in tibia of male birds of Hubbard lineage. We considered the models of plateau response of polynomial quadratic and linear exponential forms. In addition to these models, we also proposed the fitting of a logarithmic model to the data and the efficiency of the methodology was evaluated by Monte Carlo simulations, considering different scenarios for the parametric values. The isotonization of the data yielded an improvement in all the fitting quality parameters evaluated. Among the models used, the logarithmic presented estimates of the parameters more consistent with the values reported in literature.

Latent residual analysis in binary regression with skewed link

Model diagnostics is an integral part of model determination and an important part of the model diagnostics is residual analysis. We adapt and implement residuals considered in the literature for the probit, logistic and skew-probit links under binary regression. New latent residuals for the skew-probit link are proposed here. We have detected the presence of outliers using the residuals proposed here for different models in a simulated dataset and a real medical dataset.

A log-Birnbaum-Saunders regression model with asymmetric errors

This paper introduces a skewed log-Birnbaum-Saunders regression model based on the skewed sinh-normal distribution proposed by Leiva et al. [A skewed sinh-normal distribution and its properties and application to air pollution, Comm. Statist. Theory Methods 39 (2010), pp. 426-443]. Some influence methods, such as the local influence and generalized leverage, are presented. Additionally, we derived the normal curvatures of local influence under some perturbation schemes. An empirical application to a real data set is presented in order to illustrate the usefulness of the proposed model.

Multi-element determination in Brazilian honey samples by inductively coupled plasma mass spectrometry and estimation of geographic origin with data mining techniques

Multi-element analysis of honey samples was carried out with the aim of developing a reliable method of tracing the origin of honey. Forty-two chemical elements were determined (Al, Cu, Pb, Zn, Mn, Cd, Tl, Co, Ni, Rb, Ba, Be, Bi, U, V, Fe, Pt, Pd, Te, Hf, Mo, Sn, Sb, P, La, Mg, I, Sm, Tb, Dy, Sd, Th, Pr, Nd, Tm, Yb, Lu, Gd, Ho, Er, Ce, Cr) by inductively coupled plasma mass spectrometry (ICP-MS). Then, three machine learning tools for classification and two for attribute selection were applied in order to prove that it is possible to use data mining tools to find the region where honey originated. Our results clearly demonstrate the potential of Support Vector Machine (SVM), Multilayer Perceptron (MLP) and Random Forest (RF) chemometric tools for honey origin identification. Moreover, the selection tools allowed a reduction from 42 trace element concentrations to only 5. (C) 2012 Elsevier Ltd. All rights reserved.

Covariates of high-risk human papillomavirus (HPV) infections are distinct for incident CIN1, CIN2 and CIN3 as disclosed by competing-risks regression models

Background: In addition to the oncogenic human papillomavirus (HPV), several cofactors are needed in cervical carcinogenesis, but whether the HPV covariates associated with incident i) CIN1 are different from those of incident ii) CIN2 and iii) CIN3 needs further assessment. Objectives: To gain further insights into the true biological differences between CIN1, CIN2 and CIN3, we assessed HPV covariates associated with incident CIN1, CIN2, and CIN3. Study Design and Methods: HPV covariates associated with progression to CIN1, CIN2 and CIN3 were analysed in the combined cohort of the NIS (n = 3,187) and LAMS study (n = 12,114), using competing-risks regression models (in panel data) for baseline HR-HPV-positive women (n = 1,105), who represent a sub-cohort of all 1,865 women prospectively followed-up in these two studies. Results: Altogether, 90 (4.8%), 39 (2.1%) and 14 (1.4%) cases progressed to CIN1, CIN2, and CIN3, respectively. Among these baseline HR-HPV-positive women, the risk profiles of incident GIN I, CIN2 and CIN3 were unique in that completely different HPV covariates were associated with progression to CIN1, CIN2 and CIN3, irrespective which categories (non-progression, CIN1, CIN2, CIN3 or all) were used as competing-risks events in univariate and multivariate models. Conclusions: These data confirm our previous analysis based on multinomial regression models implicating that distinct covariates of HR-HPV are associated with progression to CIN1, CIN2 and CIN3. This emphasises true biological differences between the three grades of GIN, which revisits the concept of combining CIN2 with CIN3 or with CIN1 in histological classification or used as a common end-point, e.g., in HPV vaccine trials.

Geographic differentiation of polymorphism in the Plasmodium falciparum malaria vaccine candidate gene SERA5

SERA5 is regarded as a promising malaria vaccine candidate of the most virulent human malaria parasite Plasmodium falciparum. SERA5 is a 120 kDa abundantly expressed blood-stage protein containing a papain-like protease. Since substantial polymorphism in blood-stage vaccine candidates may potentially limit their efficacy, it is imperative to fully investigate polymorphism of the SERA5 gene (sera5). In this study, we performed evolutionary and population genetic analysis of sera5. The level of inter-species divergence (kS = 0.076) between P. falciparum and Plasmodium reichenowi, a closely related chimpanzee malaria parasite is comparable to that of housekeeping protein genes. A signature of purifying selection was detected in the proenzyme and enzyme domains. Analysis of 445 near full-length P. falciparum sera5 sequences from nine countries in Africa, Southeast Asia, Oceania and South America revealed extensive variations in the number of octamer repeat (OR) and serine repeat (SR) regions as well as substantial level of single nucleotide polymorphism (SNP) in non-repeat regions (2562 bp). Remarkably, a 14 amino acid sequence of SERA5 (amino acids 59-72) that is known to be the in vitro target of parasite growth inhibitory antibodies was found to be perfectly conserved in all 445 worldwide isolates of P. falciparum evaluated. Unlike other major vaccine target antigen genes such as merozoite surface protein-1, apical membrane antigen-1 or circumsporozoite protein, no strong evidence for positive selection was detected for SNPs in the non-repeat regions of sera5. A biased geographical distribution was observed in SNPs as well as in the haplotypes of the sera5 OR and SR regions. In Africa, OR- and SR-haplotypes with low frequency (<5%) and SNPs with minor allele frequency (<5%) were abundant and were mostly continent-specific. Consistently, significant genetic differentiation, assessed by the Wright's fixation index (FST) of inter-population variance in allele frequencies, was detected for SNPs and both OR- and SR-haplotypes among almost all parasite populations. The exception was parasite populations between Tanzania and Ghana, suggesting frequent gene flow in Africa. The present study points to the importance of investigating whether biased geographical distribution for SNPs and repeat variants in the OR and SR regions affect the reactivity of human serum antibodies to variants. (C) 2011 Elsevier Ltd. All rights reserved.

The log-exponentiated generalized gamma regression model for censored data

For the first time, we introduce a generalized form of the exponentiated generalized gamma distribution [Cordeiro et al. The exponentiated generalized gamma distribution with application to lifetime data, J. Statist. Comput. Simul. 81 (2011), pp. 827-842.] that is the baseline for the log-exponentiated generalized gamma regression model. The new distribution can accommodate increasing, decreasing, bathtub- and unimodal-shaped hazard functions. A second advantage is that it includes classical distributions reported in the lifetime literature as special cases. We obtain explicit expressions for the moments of the baseline distribution of the new regression model. The proposed model can be applied to censored data since it includes as sub-models several widely known regression models. It therefore can be used more effectively in the analysis of survival data. We obtain maximum likelihood estimates for the model parameters by considering censored data. We show that our extended regression model is very useful by means of two applications to real data.

Understanding bias in geographic range size estimates

Aim Estimates of geographic range size derived from natural history museum specimens are probably biased for many species. We aim to determine how bias in these estimates relates to range size. Location We conducted computer simulations based on herbarium specimen records from localities ranging from the southern United States to northern Argentina. Methods We used theory on the sampling distribution of the mean and variance to develop working hypotheses about how range size, defined as area of occupancy (AOO), was related to the inter-specific distribution of: (1) mean collection effort per area across the range of a species (MC); (2) variance in collection effort per area across the range of a species (VC); and (3) proportional bias in AOO estimates (PBias: the difference between the expected value of the estimate of AOO and true AOO, divided by true AOO). We tested predictions from these hypotheses using computer simulations based on a dataset of more than 29,000 herbarium specimen records documenting occurrences of 377 plant species in the tribe Bignonieae (Bignoniaceae). Results The working hypotheses predicted that the mean of the inter-specific distribution of MC, VC and PBias were independent of AOO, but that the respective variance and skewness decreased with increasing AOO. Computer simulations supported all but one prediction: the variance of the inter-specific distribution of VC did not decrease with increasing AOO. Main conclusions Our results suggest that, despite an invariant mean, the dispersion and symmetry of the inter-specific distribution of PBias decreases as AOO increases. As AOO increased, range size was less severely underestimated for a large proportion of simulated species. However, as AOO increased, range size estimates having extremely low bias were less common.

A New Method to Predict the Epidemiology of Fungal Keratitis by Monitoring the Sales Distribution of Antifungal Eye Drops in Brazil

Purpose: Fungi are a major cause of keratitis, although few medications are licensed for their treatment. The aim of this study is to observe the variation in commercialisation of antifungal eye drops, and to predict the seasonal distribution of fungal keratitis in Brazil. Methods: Data from a retrospective study of antifungal eye drops sales from the only pharmaceutical ophthalmologic laboratory, authorized to dispense them in Brazil (Opthalmos) were gathered. These data were correlated with geographic and seasonal distribution of fungal keratitis in Brazil between July 2002 and June 2008. Results: A total of 26,087 antifungal eye drop units were sold, with a mean of 2.3 per patient. There was significant variation in antifungal sales during the year (p < 0.01). A linear regression model displayed a significant association between reduced relative humidity and antifungal drug sales (R-2 = 0.17, p < 0.01). Conclusions: Antifungal eye drops sales suggest that there is a seasonal distribution of fungal keratitis. A possible interpretation is that the third quarter of the year (a period when the climate is drier), when agricultural activity is more intense in Brazil, suggests a correlation with a higher incidence of fungal keratitis. A similar model could be applied to other diseases, that are managed with unique, or few, and monitorable medications to predict epidemiological aspects.

Tráfego veicular e mortalidade por doenças do aparelho circulatório em homens adultos

OBJETIVO: Analisar a associação entre indicadores de exposição à poluição por tráfego veicular e mortalidade por doenças do aparelho circulatório em homens adultos. MÉTODOS: Foram analisadas informações sobre vias e volume de tráfego no ano de 2007 fornecidas pela companhia de engenharia de tráfego local. Mortalidade por doenças do aparelho circulatório no ano de 2005 entre homens ≥ 40 anos foram obtidas do registro de mortalidade do Programa de Aprimoramento de Informações de Mortalidade do Município de São Paulo, SP. Dados socioeconômicos do Censo 2000 e informações sobre a localização dos serviços de saúde também foram coletados. A exposição foi avaliada pela densidade de vias e volume de tráfego para cada distrito administrativo. Foi calculada regressão (α = 5%) entre esses indicadores de exposição e as taxas de mortalidade padronizadas, ajustando os modelos para variáveis socioeconômicas, número de serviços de saúde nos distritos e autocorrelação espacial. RESULTADOS: A correlação entre densidade de vias e volume de tráfego foi modesta (r² = 0,28). Os distritos do centro apresentaram os maiores valores de densidade de vias. O modelo de regressão espacial de densidade de vias indicou associação com mortalidade por doenças do aparelho circulatório (p = 0,017). Não se observou associação no modelo de volume de tráfego. Em ambos os modelos – vias e volume de tráfego (veículos leves/pesados) – a variável socioeconômica foi estatisticamente signifi cante. CONCLUSÕES: A associação entre mortalidade por doenças do aparelho circulatório e densidade de vias converge com a literatura e encoraja a realização de mais estudos epidemiológicos em nível individual e com métodos mais acurados de avaliação da exposição.

AGE INFLUENCE ON THE HEART RATE BEHAVIOR ON THE REST-EXERCICIO TRANSITION: AN ANALYSIS BY DELTAS AND LINEAR REGRESSION

Background: Changes in heart rate during rest-exercise transition can be characterized by the application of mathematical calculations, such as deltas 0-10 and 0-30 seconds to infer on the parasympathetic nervous system and linear regression and delta applied to data range from 60 to 240 seconds to infer on the sympathetic nervous system. The objective of this study was to test the hypothesis that young and middle-aged subjects have different heart rate responses in exercise of moderate and intense intensity, with different mathematical calculations. Methods: Seven middle-aged men and ten young men apparently healthy were subject to constant load tests (intense and moderate) in cycle ergometer. The heart rate data were submitted to analysis of deltas (0-10, 0-30 and 60-240 seconds) and simple linear regression (60-240 seconds). The parameters obtained from simple linear regression analysis were: intercept and slope angle. We used the Shapiro-Wilk test to check the distribution of data and the "t" test for unpaired comparisons between groups. The level of statistical significance was 5%. Results: The value of the intercept and delta 0-10 seconds was lower in middle age in two loads tested and the inclination angle was lower in moderate exercise in middle age. Conclusion: The young subjects present greater magnitude of vagal withdrawal in the initial stage of the HR response during constant load exercise and higher speed of adjustment of sympathetic response in moderate exercise.