On estimation and influence diagnostics for zero-inflated negative binomial regression models

The zero-inflated negative binomial model is used to account for overdispersion detected in data that are initially analyzed under the zero-Inflated Poisson model A frequentist analysis a jackknife estimator and a non-parametric bootstrap for parameter estimation of zero-inflated negative binomial regression models are considered In addition an EM-type algorithm is developed for performing maximum likelihood estimation Then the appropriate matrices for assessing local influence on the parameter estimates under different perturbation schemes and some ways to perform global influence analysis are derived In order to study departures from the error assumption as well as the presence of outliers residual analysis based on the standardized Pearson residuals is discussed The relevance of the approach is illustrated with a real data set where It is shown that zero-inflated negative binomial regression models seems to fit the data better than the Poisson counterpart (C) 2010 Elsevier B V All rights reserved

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 impact of El Nino Southern Oscillation phases on off-season maize yield for a subtropical region of Brazil

In recent years, maize has become one of the main alternative crops for the autumn winter growing season in the central-western and southeastern regions of Brazil. However, water deficits, sub-optimal temperatures and low solar radiation levels are common problems that are experienced during this growing season by local farmers. One methodology to assess the impact of variable weather conditions on crop production is the use of crop simulation models. The goal of this study was to evaluate the effect of climate variability on maize yield for a subtropical region of Brazil. Specific objectives for this study were (1) to analyse the effect of El Nino Southern Oscillation (ENSO) on precipitation and air temperature for four locations in the state of Sao Paulo and (2) to analyse the impact of ENSO on maize grown off-season for the same four locations using a crop simulation model. For each site, historical weather data were categorised as belonging to one of three phases of ENSO: El Nino (warm sea surface temperature anomalies in the Pacific), La Nina (cool sea surface temperature anomalies) or neutral, based on an index derived from observed sea surface temperature anomalies. During El Nino, there is a tendency for an increase in the rainfall amount during May for the four selected locations, and also during April, mainly in three of the locations, resulting in an increase in simulated maize yield planted between February 15 and March 15. In general, there was a decrease in the simulated yield for maize grown off-season during neutral years. This study showed how a crop model can be used to assess the impact of climate variability on the yield of maize grown off-season in a subtropical region of Brazil. The outcomes of this study can be very useful for both policy makers and local farmers for agricultural planning and decision making. Copyright (C) 2009 Royal Meteorological Society

Pricing farm-level agricultural insurance: a Bayesian approach

This paper applies Hierarchical Bayesian Models to price farm-level yield insurance contracts. This methodology considers the temporal effect, the spatial dependence and spatio-temporal models. One of the major advantages of this framework is that an estimate of the premium rate is obtained directly from the posterior distribution. These methods were applied to a farm-level data set of soybean in the State of the Parana (Brazil), for the period between 1994 and 2003. The model selection was based on a posterior predictive criterion. This study improves considerably the estimation of the fair premium rates considering the small number of observations.

Parametric and nonparametric statistical modelling of crop yield: implications for pricing crop insurance contracts

This article considers alternative methods to calculate the fair premium rate of crop insurance contracts based on county yields. The premium rate was calculated using parametric and nonparametric approaches to estimate the conditional agricultural yield density. These methods were applied to a data set of county yield provided by the Statistical and Geography Brazilian Institute (IBGE), for the period of 1990 through 2002, for soybean, corn and wheat, in the State of Paran. In this article, we propose methodological alternatives to pricing crop insurance contracts resulting in more accurate premium rates in a situation of limited data.

Crop rotation in center-pivot for phytonematode control: density variation, pathogenicity and crop loss estimation

Crop rotation in center-pivot for phytonematode control: density variation, pathogenicity and crop loss estimation A field study conducted over three consecutive years, on a farm using crop rotation system under center-pivot and infested with the nematodes Pratylenchus brachyurus, P. zeae, Meloidogyne incognita, Paratrichodorus minor, Helicotylenchus dihystera, Mesocriconema ornata and M. onoense, demonstrated that intensive crop systems provide conditions for the maintenance of high densities of polyphagous phytonematodes. Of the crops established on the farm (cotton, maize, soybean and cowpea), cotton and soybean suffered the most severe crop losses, caused respectively by M. incognita and P. brachyurus. Since maize is a good host for both nematodes, but tolerant of M. incognita, its exclusion from cropping system would be favorable to the performance of cotton, soybean and cowpea. Results from experiments carried out in controlled conditions confirmed the pathogenicity of P. brachyurus on cotton. Additional management with genetic resistance was useful in fields infested with M. incognita, although the soybean performance was affected by low resistance of the cultivars used for P. brachyurus. In conclusion, crop rotation must be carefully planned in areas infested with polyphagous nematodes, specifically in the case of occurrence of two or more major pathogenic nematodes.

Yield loss caused by huanglongbing in different sweet orange cultivars in So Paulo, Brazil

Citrus huanglongbing (HLB) reduces an affected orchard`s economic life. This work aimed to characterize yield loss due to HLB for different sweet orange cultivars and determine the relationship between disease severity and yield. Disease severity and yield were assessed on 949 individual trees distributed in 11 different blocks from sweet orange cultivars Hamlin, Westin, Pera and Valencia. In each block, plants showing a range of HLB severity levels and asymptomatic plants were selected. Total yield (weight of harvested fruit), mean weight of asymptomatic and symptomatic fruit, relative yield (symptomatic tree yield/mean yield of asymptomatic trees from the same block) and relative number of fruits (fruit number from symptomatic tree/mean number of fruits from asymptomatic trees from the same block) were determined. The weight of symptomatic fruit was lower than the weight of asymptomatic fruit, but the weights of asymptomatic and symptomatic fruit were not correlated with disease severity, indicating that the effects of HLB were restricted to symptomatic branches. The relationship of the relative yield with HLB severity can be satisfactorily described by a negative exponential model. The rates of yield decrease as a function of disease severity were similar for all assessed cultivars. A relative yield (up to 19%) was observed even for trees where disease severity was 100%. The strong linear relationship between relative number of fruits per tree and the relative yield per tree suggested that the yield reduction was due primarily to early fruit drop or lack of fruit set on affected branches.

Role of healthy nursery plants in orange yield during eight years of Citrus Variegated Chlorosis epidemics

The production and commercialization of citrus seedlings inspected and produced in protected screen-houses has become mandatory in Sao Paulo State, Brazil since January 2003. This law was intended to avoid the dispersion of Citrus Variegated Chlorosis (CVC), disease caused by Xylella fastidiosa. Our objective was to compare the yield over 8 years of `Natal` sweet orange trees grafted onto Rangpur lime obtained from healthy nursery plants and from plants artificially inoculated with X. fastidiosa. Yield was evaluated in an orchard planted in February 1999 with two treatments: (i) trees from healthy nursery plant, and (ii) trees from plants artificially inoculated with X. fastidiosa. The mean yield was 21% higher in trees from healthy nursery plants, as compared to trees from inoculated nursery plants. This difference represents a gain of approximately 203 boxes of 40.8 kg each, considering a planting density of 550 plants per hectare. (C) 2011 Elsevier B.V. All rights reserved.

Identification of microsatellites markers associated with yield components and quality parameters in sugarcane

Marker assisted selection depends on the identification of tightly linked association between marker and the trait of interest. In the present work, functional (EST-SSRs) and genomic (gSSRs) microsatellite markers were used to detect putative QTLs for sugarcane yield components (stalk number, diameter and height) and as well as for quality parameters (Brix, Pol and fibre) in plant cane. The mapping population (200 individuals) was derived from a bi-parental cross (IACSP95-3018 x IACSP93-3046) from the IAC Sugarcane Breeding Program. As the map is under construction, single marker trait association analysis based on the likelihood ratio test was undertaken to detect the QTLs. Of the 215 single dose markers evaluated (1:1 and 3:1), 90 (42%) were associated with putative QTLs involving 43 microsatellite primers (18 gSSRs and 25 EST-SSRs). For the yield components, 41 marker/trait associations were found: 20 for height, 6 for diameter and 15 for stalk number. An EST-SSRs marker with homology to non-phototropic hypocotyls 4 (NPH4) protein was associated with a putative QTL with positive effect for diameter as also with a negative effect for stalk number. In relation to the quality parameters, 18 marker trait associations were found for Brix, 19 for Pol, and 12 for fibre. For fibre, 58% of the QTLs detected showed a negative effect on this trait. Some makers associated with QTLs with a negative effect for fibre showed a positive effect for Pol, reflecting the negative correlation generally observed between these traits.

Analysis of genomic and functional RFLP derived markers associated with sucrose content, fiber and yield QTLs in a sugarcane (Saccharum spp.) commercial cross

Expressed sequence tags derived markers have a great potential to be used in functional map construction and QTL tagging. In the present work, sugarcane genomic probes and expressed sequence tags having homology to genes, mostly involved in carbohydrate metabolism were used in RFLP assays to identify putative QTLs as well as their epistatic interactions for fiber content, cane yield, pol and tones of sugar per hectare, at two crop cycles in a progeny derived from a bi-parental cross of sugarcane elite materials. A hundred and twenty marker trait associations were found, of which 26 at both crop cycle and 32 only at first ratoon cane. A sucrose synthase derived marker was associated with a putative QTL having a high negative effect on cane yield and also with a QTL having a positive effect on Pol at both crop cycles. Fifty digenic epistatic marker interactions were identified for the four traits evaluated. Of these, only two were observed at both crop cycles.

An integrated molecular map of yellow passion fruit based on simultaneous maximum-likelihood estimation of linkage and linkage phases

The development of genetic maps for auto-incompatible species, such as the yellow passion fruit (Passiflora edulis Sims f.flavicarpa Deg.) is restricted due to the unfeasibility of obtaining traditional mapping populations based on inbred lines. For this reason, yellow passion fruit linkage maps were generally constructed using a strategy known as two-way pseudo-testeross, based on monoparental dominant markers segregating in a 1:1 fashion. Due to the lack of information from these markers in one of the parents, two individual (parental) maps were obtained. However, integration of these maps is essential, and biparental markers can be used for such an operation. The objective of our study was to construct an integrated molecular map for a full-sib population of yellow passion fruit combining different loci configuration generated from amplified fragment length polymorphisms (AFLPs) and microsatellite markers and using a novel approach based on simultaneous maximum-likelihood estimation of linkage and linkage phases, specially designed for outcrossing species. Of the total number of loci, approximate to 76%, 21%, 0.7%, and 2.3% did segregate in 1:1, 3:1, 1:2:1, and 1:1:1:1 ratios, respectively. Ten linkage groups (LGs) were established with a logarithm of the odds (LOD) score >= 5.0 assuming a recombination fraction : <= 0.35. On average, 24 markers were assigned per LG, representing a total map length of 1687 cM, with a marker density of 6.9 cM. No markers were placed as accessories on the map as was done with previously constructed individual maps.

Tree performance and fruit yield and quality of `Okitsu` Satsuma mandarin grafted on 12 rootstocks

The citriculture in Brazil, as well as in other important regions in the world, is based on very few mandarin cultivars. This fact leads to a short harvest period and higher prices for off-season fruit. The `Okitsu` Satsuma (Citrus unshiu Marc.) is among the earliest ripening mandarin cultivars and it is considered to be tolerant to, citrus canker (Xanthomonas citri subsp. citri Schaad et al.) and to citrus variegated chlorosis (Xylella fastidiosa Wells et al.). Despite having regular fruit quality under hot climate conditions, the early fruit maturation and absence of seeds of `Okitsu` fruits are well suited for the local market in the summer(December through March), when the availability of citrus fruits for fresh consumption is limited. Yet, only a few studies have been conducted in Brazil on rootstocks for `Okitsu`. Consequently, a field trial was carried out in Bebeclouro, Sao Paulo State, to evaluate the horticultural performance of `Okitsu` Satsuma mandarin budded onto 12 rootstocks: the citrandarin `Changsha` mandarin (Citrus reticulata Blanco) x Poncirus trifoliata `English Small`: the hybrid Rangpur lime (Citrus limonia Osbeck) x `Swingle` citrumelo (P. trifoliata (L.) Raf. x Citrus paradisi Macfad.); the trifoliates (P. trifoliata (L) Raf)`Rubidoux`,`FCAV` and `Flying Dragon`(P. trifoliata var. monstrosa); the mandarins `Sun Chu Sha Kat`(C. reticulata Blanco) and `Sunki`(Citrus sunki (Hayata) Hort. ex. Tanaka); the Rangpur limes (C. limonia Osbeck) `Cravo Limeira` and `Cravo FCAV`;`Carrizo` citrange (Citrus sinensis x P. trifoliata), `Swingle` citrumelo (P. trifoliata x C. paradisi), and `Orlando` tangelo (C. paradisi x Citrus tangerina cv. `Dancy`). The experimental grove was planted in 2001, using a 6 m x 3 m spacing, in a randomized block design. No supplementary irrigation was applied. Fruit yield, canopy volume, and fruit quality were assessed for each rootstock. A cluster multivariate analysis identified three different rootstock pairs with similar effects on plant growth, yield and fruit quality of `Okitsu` mandarin. The `Flying Dragon `trifoliate had a unique effect over the `Okitsu` trees performance, inducing lower canopy volume and higher yield efficiency and fruit quality, and might be suitable for high-density plantings. The `Cravo Limeira` and `Cravo FCAV` Rangpur limes induced early-ripening of fruits, with low fruit quality. `Sun Chu Sha Kat` and `Sunki` mandarins and the `Orlando` tangelo conferred lower yield efficiency and less content of soluble solids for the latter rootstock. (C) 2009 Elsevier B.V. All rights reserved.

Seed development, yield and quality of two palm species growing in different tropical forest types in SE Brazil: implications for ecological restoration

Natural forest remnants have been set as seed production fields to supply seeds of native tree species for tropical forest restoration, but the effect of different forest types on seed production has not been accessed to date for palm species. In this work, we studied seed development, yield, and quality of two palm species in different tropical forest types in SE Brazil. Seed production of palmiteiro (Euterpe edulis) and queen-palm (Syagrus romanzoffiana), which are largely used in restoration efforts due to their importance for vertebrate frugivores, were studied in natural remnants of Atlantic Rainforest, Restinga Forest, Seasonally Dry Forest, and Cerrado Forest. We studied seed development, yield, size, and germination of seed lots produced in some of these forest types, including seeds harvested in 2008, 2009, and both years. Seed yield and quality, as well as seed dry mass in 2009, were higher for palmiteiro seeds produced in the Atlantic Rainforest, while queen-palm seeds produced at the Restinga Forest showed the higher mass and yield, but the lowest physiological potential. Consequently, these natural differences of seed yield and quality have to be taken into account for establishing standards for seed commercialization and analysis, seed pricing, and seedling production in forest nurseries.

Parametric correlation functions to model the structure of permanent environmental (co)variances in milk yield random regression models

The objective of the present study was to estimate milk yield genetic parameters applying random regression models and parametric correlation functions combined with a variance function to model animal permanent environmental effects. A total of 152,145 test-day milk yields from 7,317 first lactations of Holstein cows belonging to herds located in the southeastern region of Brazil were analyzed. Test-day milk yields were divided into 44 weekly classes of days in milk. Contemporary groups were defined by herd-test-day comprising a total of 2,539 classes. The model included direct additive genetic, permanent environmental, and residual random effects. The following fixed effects were considered: contemporary group, age of cow at calving (linear and quadratic regressions), and the population average lactation curve modeled by fourth-order orthogonal Legendre polynomial. Additive genetic effects were modeled by random regression on orthogonal Legendre polynomials of days in milk, whereas permanent environmental effects were estimated using a stationary or nonstationary parametric correlation function combined with a variance function of different orders. The structure of residual variances was modeled using a step function containing 6 variance classes. The genetic parameter estimates obtained with the model using a stationary correlation function associated with a variance function to model permanent environmental effects were similar to those obtained with models employing orthogonal Legendre polynomials for the same effect. A model using a sixth-order polynomial for additive effects and a stationary parametric correlation function associated with a seventh-order variance function to model permanent environmental effects would be sufficient for data fitting.

Random regression models to estimate test-day milk yield genetic parameters Holstein cows in Southeastern Brazil

A total of 152,145 weekly test-day milk yield records from 7317 first lactations of Holstein cows distributed in 93 herds in southeastern Brazil were analyzed. Test-day milk yields were classified into 44 weekly classes of DIM. The contemporary groups were defined as herd-year-week of test-day. The model included direct additive genetic, permanent environmental and residual effects as random and fixed effects of contemporary group and age of cow at calving as covariable, linear and quadratic effects. Mean trends were modeled by a cubic regression on orthogonal polynomials of DIM. Additive genetic and permanent environmental random effects were estimated by random regression on orthogonal Legendre polynomials. Residual variances were modeled using third to seventh-order variance functions or a step function with 1, 6,13,17 and 44 variance classes. Results from Akaike`s and Schwarz`s Bayesian information criterion suggested that a model considering a 7th-order Legendre polynomial for additive effect, a 12th-order polynomial for permanent environment effect and a step function with 6 classes for residual variances, fitted best. However, a parsimonious model, with a 6th-order Legendre polynomial for additive effects and a 7th-order polynomial for permanent environmental effects, yielded very similar genetic parameter estimates. (C) 2008 Elsevier B.V. All rights reserved.