Selecting statistical models to study the relationship between soybean yield and soil physical properties

Statistical models allow the representation of data sets and the estimation and/or prediction of the behavior of a given variable through its interaction with the other variables involved in a phenomenon. Among other different statistical models, are the autoregressive state-space models (ARSS) and the linear regression models (LR), which allow the quantification of the relationships among soil-plant-atmosphere system variables. To compare the quality of the ARSS and LR models for the modeling of the relationships between soybean yield and soil physical properties, Akaike's Information Criterion, which provides a coefficient for the selection of the best model, was used in this study. The data sets were sampled in a Rhodic Acrudox soil, along a spatial transect with 84 points spaced 3 m apart. At each sampling point, soybean samples were collected for yield quantification. At the same site, soil penetration resistance was also measured and soil samples were collected to measure soil bulk density in the 0-0.10 m and 0.10-0.20 m layers. Results showed autocorrelation and a cross correlation structure of soybean yield and soil penetration resistance data. Soil bulk density data, however, were only autocorrelated in the 0-0.10 m layer and not cross correlated with soybean yield. The results showed the higher efficiency of the autoregressive space-state models in relation to the equivalent simple and multiple linear regression models using Akaike's Information Criterion. The resulting values were comparatively lower than the values obtained by the regression models, for all combinations of explanatory variables.

The three phases of the ensemble forecasting of niche models: geographic range and shifts in climatically suitable areas of Utetheisa ornatrix (Lepidoptera, Arctiidae)

Species' geographic ranges are usually considered as basic units in macroecology and biogeography, yet it is still difficult to measure them accurately for many reasons. About 20 years ago, researchers started using local data on species' occurrences to estimate broad scale ranges, thereby establishing the niche modeling approach. However, there are still many problems in model evaluation and application, and one of the solutions is to find a consensus solution among models derived from different mathematical and statistical models for niche modeling, climatic projections and variable combination, all of which are sources of uncertainty during niche modeling. In this paper, we discuss this approach of ensemble forecasting and propose that it can be divided into three phases with increasing levels of complexity. Phase I is the simple combination of maps to achieve a consensual and hopefully conservative solution. In Phase II, differences among the maps used are described by multivariate analyses, and Phase III consists of the quantitative evaluation of the relative magnitude of uncertainties from different sources and their mapping. To illustrate these developments, we analyzed the occurrence data of the tiger moth, Utetheisa ornatrix (Lepidoptera, Arctiidae), a Neotropical moth species, and modeled its geographic range in current and future climates.

Impact of livestock on a mosquito community (Diptera: Culicidae) in a Brazilian tropical dry forest

AbstractINTRODUCTION: This study evaluated the effects of cattle removal on the Culicidae mosquito community structure in a tropical dry forest in Brazil.METHODS: Culicidae were collected during dry and wet seasons in cattle presence and absence between August 2008 and October 2010 and assessed using multivariate statistical models.RESULTS: Cattle removal did not significantly alter Culicidae species richness and abundance. However, alterations were noted in Culicidae community composition.CONCLUSIONS: This is the first study to evaluate the impact of cattle removal on Culicidae community structure in Brazil and demonstrates the importance of assessing ecological parameters such as community species composition.

Prevalence of self-reported constipation in adults from the general population

OBJECTIVE To estimate the prevalence of self-reported constipation and associated factors in the general population of a Brazilian city. METHOD Secondary analysis of an epidemiological study, population-based, cross-sectional study, about bowel habits of Brazilian population. A total of 2,162 individuals were interviewed using two instruments: sociodemographic data and the adapted and validated Brazilian version of the "Bowel Function in the Community" tool. RESULTS There was a prevalence of 25.2% for the self-reported constipation, 37.2% among women and 10.2% among men. Stroke and old age were associated with constipation in the three statistical models used. CONCLUSION The prevalence found showed to be similar to the findings in the literature, although some associated factors obtained here have never been investigated.

Beaufils ranges to assess the cotton nutrient status in the southern region of Mato Grosso

The relationships between nutrient contents and indices of the Diagnosis and Recommendation Integrated System (DRIS) are a useful basis to determine appropriate ranges for the interpretation of leaf nutrient contents. The purpose of this study was to establish Beaufils ranges from statistical models of the relationship between foliar concentrations and DRIS indices, generated by two systems of DRIS norms - the F value and natural logarithm transformation - and assess the nutritional status of cotton plants, based on these Beaufils ranges. Yield data from plots (average acreage 100 ha) and foliar concentrations of macro and micronutrients of cotton (Gossypium hirsutum r. latifolium) plants, in the growing season 2004/2005, were stored in a database. The criterion to define the reference population consisted of plots with above-average yields + 0.5 standard deviation (over 4,575 kg ha-1 seed cotton yield). The best-fitting statistical model of the relationship between foliar nutrient concentrations and DRIS indices was linear, with R² > 0.8090, p < 0.01, except for N, with R² = 0.5987, p < 0.01. The two criteria were effective to diagnose the plant nutritional status. The diagnoses were not random, but based on the effectiveness of the chi-square-tested method. The agreement between the methods to assess the nutritional status was 92.59-100 %, except for S, with 74.07 % agreement.

Genetic associations in the detection of QTLs for wheat spike-related traits

The objective of this work was to assess the genetic diversity and population structure of wheat genotypes, to detect significant and stable genetic associations, as well as to evaluate the efficiency of statistical models to identify chromosome regions responsible for the expression of spike-related traits. Eight important spike characteristics were measured during five growing seasons in Serbia. A set of 30 microsatellite markers positioned near important agronomic loci was used to evaluate genetic diversity, resulting in a total of 349 alleles. The marker-trait associations were analyzed using the general linear and mixed linear models. The results obtained for number of allelic variants per locus (11.5), average polymorphic information content value (0.68), and average gene diversity (0.722) showed that the exceptional level of polymorphism in the genotypes is the main requirement for association studies. The population structure estimated by model-based clustering distributed the genotypes into six subpopulations according to log probability of data. Significant and stable associations were detected on chromosomes 1B, 2A, 2B, 2D, and 6D, which explained from 4.7 to 40.7% of total phenotypic variations. The general linear model identified a significantly larger number of marker-trait associations (192) than the mixed linear model (76). The mixed linear model identified nine markers associated to six traits.

Annonaceae substitution rates: a codon model perspective

The Annonaceae includes cultivated species of economic interest and represents an important source of information for better understanding the evolution of tropical rainforests. In phylogenetic analyses of DNA sequence data that are used to address evolutionary questions, it is imperative to use appropriate statistical models. Annonaceae are cases in point: Two sister clades, the subfamilies Annonoideae and Malmeoideae, contain the majority of Annonaceae species diversity. The Annonoideae generally show a greater degree of sequence divergence compared to the Malmeoideae, resulting in stark differences in branch lengths in phylogenetic trees. Uncertainty in how to interpret and analyse these differences has led to inconsistent results when estimating the ages of clades in Annonaceae using molecular dating techniques. We ask whether these differences may be attributed to inappropriate modelling assumptions in the phylogenetic analyses. Specifically, we test for (clade-specific) differences in rates of non-synonymous and synonymous substitutions. A high ratio of nonsynonymous to synonymous substitutions may lead to similarity of DNA sequences due to convergence instead of common ancestry, and as a result confound phylogenetic analyses. We use a dataset of three chloroplast genes (rbcL, matK, ndhF) for 129 species representative of the family. We find that differences in branch lengths between major clades are not attributable to different rates of non-synonymous and synonymous substitutions. The differences in evolutionary rate between the major clades of Annonaceae pose a challenge for current molecular dating techniques that should be seen as a warning for the interpretation of such results in other organisms.

Rainfall erosivity for the State of Rio de Janeiro estimated by artificial neural network

The Artificial Neural Networks (ANNs) are mathematical models method capable of estimating non-linear response plans. The advantage of these models is to present different responses of the statistical models. Thus, the objective of this study was to develop and to test ANNs for estimating rainfall erosivity index (EI30) as a function of the geographical location for the state of Rio de Janeiro, Brazil and generating a thematic visualization map. The characteristics of latitude, longitude e altitude using ANNs were acceptable to estimating EI30 and allowing visualization of the space variability of EI30. Thus, ANN is a potential option for the estimate of climatic variables in substitution to the traditional methods of interpolation.

Economic analysis of cashew early dwarf Crop BRS - 189 depending on water levels and doses of potassium fertilization

Water and fertilizer among the production factors are the elements that most restrict the production of cashew. The precise amount of these factors is essential to the success of the crop yield. This research aimed to determine the best factor-product ratio and analyze technical and economic indicators, of productivity of the cashew clone BRS 189 (Anacardium occidentale) to production factors water and potassium. The experiment was conducted from May 2009 to December 2009 in an experimental area of 56.0 m x 112.0 m in the irrigated Curu - Pentecoste, located in the municipality of Pentecoste, Ceará, Brazil. Production factors water (W) and potassium (K) were the independent variables and productivity (Y), the dependent variable. Ten statistical models that have proven satisfactory for obtaining production function were tested. The marginal rate of substitution was obtained through the ratio of the potassium marginal physical product and the water marginal physical product. The most suited model to the conditions of the experiment was the quadratic polynomial without intercept and interaction. Considering that the price of the water was 0.10 R$ mm -1, the price of the potassium 2.19 R$ kg -1 and the price of the cashew 0.60 R$ kg-1, the amounts of water and K2O to obtain the maximum net income were 6,349.1 L plant-1 of water and 128.7 g plant -1year, -1 respectively. Substituting the values obtained in the production function, the maximum net income was achieved with a yield of 7,496.8 kg ha-1 of cashew.

The effect of performing corrections on reported uterine cancer mortality data in the city of São Paulo

Reports of uterine cancer deaths that do not specify the subsite of the tumor threaten the quality of the epidemiologic appraisal of corpus and cervix uteri cancer mortality. The present study assessed the impact of correcting the estimated corpus and cervix uteri cancer mortality in the city of São Paulo, Brazil. The epidemiologic assessment of death rates comprised the estimation of magnitudes, trends (1980-2003), and area-level distribution based on three strategies: i) using uncorrected death certificate information; ii) correcting estimates of corpus and cervix uteri mortality by fully reallocating unspecified deaths to either one of these categories, and iii) partially correcting specified estimates by maintaining as unspecified a fraction of deaths certified as due to cancer of "uterus not otherwise specified". The proportion of uterine cancer deaths without subsite specification decreased from 42.9% in 1984 to 20.8% in 2003. Partial and full corrections resulted in considerable increases of cervix (31.3 and 48.8%, respectively) and corpus uteri (34.4 and 55.2%) cancer mortality. Partial correction did not change trends for subsite-specific uterine cancer mortality, whereas full correction did, thus representing an early indication of decrease for cervical neoplasms and stability for tumors of the corpus uteri in this population. Ecologic correlations between mortality and socioeconomic indices were unchanged for both strategies of correcting estimates. Reallocating unspecified uterine cancer mortality in contexts with a high proportion of these deaths has a considerable impact on the epidemiologic profile of mortality and provides more reliable estimates of cervix and corpus uteri cancer death rates and trends.

Evaluation of statistical and geostatistical models of digital soil properties mapping in tropical mountain regions

Soil properties have an enormous impact on economic and environmental aspects of agricultural production. Quantitative relationships between soil properties and the factors that influence their variability are the basis of digital soil mapping. The predictive models of soil properties evaluated in this work are statistical (multiple linear regression-MLR) and geostatistical (ordinary kriging and co-kriging). The study was conducted in the municipality of Bom Jardim, RJ, using a soil database with 208 sampling points. Predictive models were evaluated for sand, silt and clay fractions, pH in water and organic carbon at six depths according to the specifications of the consortium of digital soil mapping at the global level (GlobalSoilMap). Continuous covariates and categorical predictors were used and their contributions to the model assessed. Only the environmental covariates elevation, aspect, stream power index (SPI), soil wetness index (SWI), normalized difference vegetation index (NDVI), and b3/b2 band ratio were significantly correlated with soil properties. The predictive models had a mean coefficient of determination of 0.21. Best results were obtained with the geostatistical predictive models, where the highest coefficient of determination 0.43 was associated with sand properties between 60 to 100 cm deep. The use of a sparse data set of soil properties for digital mapping can explain only part of the spatial variation of these properties. The results may be related to the sampling density and the quantity and quality of the environmental covariates and predictive models used.

Optimization of solvent mixtures for extraction from bark of Schinus terebinthifolius by a statistical mixture-design technique and development of a UV-Vis spectrophotometric method for analysis of total polyphenols in the extract

A statistical mixture-design technique was used to study the effects of different solvents and their mixtures on the yield, total polyphenol content, and antioxidant capacity of the crude extracts from the bark of Schinus terebinthifolius Raddi (Anacardiaceae). The experimental results and their response-surface models showed that ternary mixtures with equal portions of all the three solvents (water, ethanol and acetone) were better than the binary mixtures in generating crude extracts with the highest yield (22.04 ± 0.48%), total polyphenol content (29.39 ± 0.39%), and antioxidant capacity (6.38 ± 0.21). An analytical method was developed and validated for the determination of total polyphenols in the extracts. Optimal conditions for the various parameters in this analytical method, namely, the time for the chromophoric reaction to stabilize, wavelength of the absorption maxima to be monitored, the reference standard and the concentration of sodium carbonate were determined to be 5 min, 780 nm, pyrogallol, and 14.06% w v-1, respectively. UV-Vis spectrophotometric monitoring of the reaction under these conditions proved the method to be linear, specific, precise, accurate, reproducible, robust, and easy to perform.

Models and applications for risk assessment and prediction of Asian soybean rust epidemics

Asian rust of soybean [Glycine max (L.) Merril] is one of the most important fungal diseases of this crop worldwide. The recent introduction of Phakopsora pachyrhizi Syd. & P. Syd in the Americas represents a major threat to soybean production in the main growing regions, and significant losses have already been reported. P. pachyrhizi is extremely aggressive under favorable weather conditions, causing rapid plant defoliation. Epidemiological studies, under both controlled and natural environmental conditions, have been done for several decades with the aim of elucidating factors that affect the disease cycle as a basis for disease modeling. The recent spread of Asian soybean rust to major production regions in the world has promoted new development, testing and application of mathematical models to assess the risk and predict the disease. These efforts have included the integration of new data, epidemiological knowledge, statistical methods, and advances in computer simulation to develop models and systems with different spatial and temporal scales, objectives and audience. In this review, we present a comprehensive discussion on the models and systems that have been tested to predict and assess the risk of Asian soybean rust. Limitations, uncertainties and challenges for modelers are also discussed.

A statistical approach to define some tofu processing conditions

The aim of this work was to make tofu from soybean cultivar BRS 267 under different processing conditions in order to evaluate the influence of each treatment on the product quality. A fractional factorial 2(5-1) design was used, in which independent variables (thermal treatment, coagulant concentration, coagulation time, curd cutting, and draining time) were tested at two different levels. The response variables studied were hardness, yield, total solids, and protein content of tofu. Polynomial models were generated for each response. To obtain tofu with desirable characteristics (hardness ~4 N, yield 306 g tofu.100 g-1 soybeans, 12 g proteins.100 g-1 tofu and 22 g solids.100 g-1 tofu), the following processing conditions were selected: heating until boiling plus 10 minutes in water bath, 2% dihydrated CaSO4 w/w, 10 minutes coagulation, curd cutting, and 30 minutes draining time.