Extending Zelterman’s approach for robust estimation of population size to zero-truncated clustered data

Estimation of population size with missing zero-class is an important problem that is encountered in epidemiological assessment studies. Fitting a Poisson model to the observed data by the method of maximum likelihood and estimation of the population size based on this fit is an approach that has been widely used for this purpose. In practice, however, the Poisson assumption is seldom satisfied. Zelterman (1988) has proposed a robust estimator for unclustered data that works well in a wide class of distributions applicable for count data. In the work presented here, we extend this estimator to clustered data. The estimator requires fitting a zero-truncated homogeneous Poisson model by maximum likelihood and thereby using a Horvitz-Thompson estimator of population size. This was found to work well, when the data follow the hypothesized homogeneous Poisson model. However, when the true distribution deviates from the hypothesized model, the population size was found to be underestimated. In the search of a more robust estimator, we focused on three models that use all clusters with exactly one case, those clusters with exactly two cases and those with exactly three cases to estimate the probability of the zero-class and thereby use data collected on all the clusters in the Horvitz-Thompson estimator of population size. Loss in efficiency associated with gain in robustness was examined based on a simulation study. As a trade-off between gain in robustness and loss in efficiency, the model that uses data collected on clusters with at most three cases to estimate the probability of the zero-class was found to be preferred in general. In applications, we recommend obtaining estimates from all three models and making a choice considering the estimates from the three models, robustness and the loss in efficiency. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)

On the clustering of winter storm loss events over Germany

During the last decades, several windstorm series hit Europe leading to large aggregated losses. Such storm series are examples of serial clustering of extreme cyclones, presenting a considerable risk for the insurance industry. Clustering of events and return periods of storm series for Germany are quantified based on potential losses using empirical models. Two reanalysis data sets and observations from German weather stations are considered for 30 winters. Histograms of events exceeding selected return levels (1-, 2- and 5-year) are derived. Return periods of historical storm series are estimated based on the Poisson and the negative binomial distributions. Over 4000 years of general circulation model (GCM) simulations forced with current climate conditions are analysed to provide a better assessment of historical return periods. Estimations differ between distributions, for example 40 to 65 years for the 1990 series. For such less frequent series, estimates obtained with the Poisson distribution clearly deviate from empirical data. The negative binomial distribution provides better estimates, even though a sensitivity to return level and data set is identified. The consideration of GCM data permits a strong reduction of uncertainties. The present results support the importance of considering explicitly clustering of losses for an adequate risk assessment for economical applications.

Destructive weighted Poisson cure rate models

In this paper, we develop a flexible cure rate survival model by assuming the number of competing causes of the event of interest to follow a compound weighted Poisson distribution. This model is more flexible in terms of dispersion than the promotion time cure model. Moreover, it gives an interesting and realistic interpretation of the biological mechanism of the occurrence of event of interest as it includes a destructive process of the initial risk factors in a competitive scenario. In other words, what is recorded is only from the undamaged portion of the original number of risk factors.

A hands-on approach for fitting long-term survival models under the GAMLSS framework

In many data sets from clinical studies there are patients insusceptible to the occurrence of the event of interest. Survival models which ignore this fact are generally inadequate. The main goal of this paper is to describe an application of the generalized additive models for location, scale, and shape (GAMLSS) framework to the fitting of long-term survival models. in this work the number of competing causes of the event of interest follows the negative binomial distribution. In this way, some well known models found in the literature are characterized as particular cases of our proposal. The model is conveniently parameterized in terms of the cured fraction, which is then linked to covariates. We explore the use of the gamlss package in R as a powerful tool for inference in long-term survival models. The procedure is illustrated with a numerical example. (C) 2009 Elsevier Ireland Ltd. All rights reserved.

COM-Poisson cure rate survival models and an application to a cutaneous melanoma data

In this paper, we develop a flexible cure rate survival model by assuming the number of competing causes of the event of interest to follow the Conway-Maxwell Poisson distribution. This model includes as special cases some of the well-known cure rate models discussed in the literature. Next, we discuss the maximum likelihood estimation of the parameters of this cure rate survival model. Finally, we illustrate the usefulness of this model by applying it to a real cutaneous melanoma data. (C) 2009 Elsevier B.V. All rights reserved.

Spatial distribution of Diaphorina citri Kuwayama (Hemiptera: Psyllidae) in citrus orchards

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Distribuição Espacial da Lagarta-do-Cartucho, Spodoptera frugiperda (J.E. Smith) (Lepidoptera: Noctuidae), na Cultura do Milho

O estudo da distribuição espacial da lagarta-do-cartucho, Spodoptera frugiperda (J.E. Smith), na cultura do milho é fundamental para garantir a utilização de estratégias de controle, otimização de técnicas de amostragem, determinação de danos econômicos e incorporação da dinâmica espacial dentro do modelo populacional. O experimento foi conduzido em três campos contendo 100 parcelas cada, sendo amostradas 10 plantas ao acaso por parcela, num total de 1000 plantas por campo em cinco datas de amostragem. Foram contados o número de lagartas pequenas (menor que 10 mm) e grandes (maior que 10 mm) por planta. As lagartas pequenas apresentaram um ajuste muito bom à distribuição binomial negativa e não à Poisson, indicando que esta categoria larval encontra-se agregada no campo. Os números de lagartas grandes por planta ajustaram-se razoavelmente à distribuição binomial negativa, com algumas datas ajustando-se à distribuição de Poisson. Portanto, as lagartas pequenas (alta densidade populacional) têm distribuição agregada no campo, enquanto que as lagartas grandes (baixa densidade populacional) podem ser mais dispersas no campo, tendendo à aleatoriedade.

DISTRIBUTION OF THE NUMBER OF EURYTREMA SP EGGS PER GRAM OF FECES IN NATURALLY INFECTED CATTLE

Eurytrema sp. egg counts (epg) in the feces of naturally infected cattle were performed and the technique employed showed 94.2% probability of detecting positive cases of the infection with a single examination independently of the host parasite burden. It was also demonstrated that the epg of Eurytrema sp. follows a negative binomial distribution model and is characterized by its small magnitude.

Density and spatial distribution of Atta sexdens rubropilosa and Atta laevigata colonies (Hym., Formicidae) in Eucalyptus spp. forests

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

The Poisson-exponential lifetime distribution

In this paper we proposed a new two-parameters lifetime distribution with increasing failure rate. The new distribution arises on a latent complementary risk problem base. The properties of the proposed distribution are discussed, including a formal proof of its probability density function and explicit algebraic formulae for its reliability and failure rate functions, quantiles and moments, including the mean and variance. A simple EM-type algorithm for iteratively computing maximum likelihood estimates is presented. The Fisher information matrix is derived analytically in order to obtaining the asymptotic covariance matrix. The methodology is illustrated on a real data set. © 2010 Elsevier B.V. All rights reserved.

Spatial distribution and architecture of acromyrmex landolti forel (Hymenoptera, Formicidae) nests in pastures of Southwestern Bahia, Brazil

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Distribuição espacial e plano de amostragem sequencial para Euschistus heros (Fabricius, 1794) (Heteroptera: Pentatomidae) na cultura da soja convencional e transgênica

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Distribuição espacial e plano de amostragem sequencial para Piezodorus guildinii (Westwood, 1837) (Hemiptera: Pentatomidae) na cultura de soja transgênica RR®

Pós-graduação em Agronomia (Entomologia Agrícola) - FCAV

Dinâmica populacional, distribuição espacial e plano de amostragem sequencial de Euschistus heros (Fabricius, 1794) (Heteroptera: Pentatomidae) em cultivares de soja

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Stochastic modeling of bacterial protein domains distribution to predict horizontal gene transfer

In questo elaborato, abbiamo tentato di modellizzare i processi che regolano la presenza dei domini proteici. I domini proteici studiati in questa tesi sono stati ottenuti dai genomi batterici disponibili nei data base pubblici (principalmente dal National Centre for Biotechnology Information: NCBI) tramite una procedura di simulazione computazionale. Ci siamo concentrati su organismi batterici in quanto in essi la presenza di geni trasmessi orizzontalmente, ossia che parte del materiale genetico non provenga dai genitori, e assodato che sia presente in una maggiore percentuale rispetto agli organismi più evoluti. Il modello usato si basa sui processi stocastici di nascita e morte, con l'aggiunta di un parametro di migrazione, usato anche nella descrizione dell'abbondanza relativa delle specie in ambito delle biodiversità ecologiche. Le relazioni tra i parametri, calcolati come migliori stime di una distribuzione binomiale negativa rinormalizzata e adattata agli istogrammi sperimentali, ci induce ad ipotizzare che le famiglie batteriche caratterizzate da un basso valore numerico del parametro di immigrazione abbiano contrastato questo deficit con un elevato valore del tasso di nascita. Al contrario, ipotizziamo che le famiglie con un tasso di nascita relativamente basso si siano adattate, e in conseguenza, mostrano un elevato valore del parametro di migrazione. Inoltre riteniamo che il parametro di migrazione sia direttamente proporzionale alla quantità di trasferimento genico orizzontale effettuato dalla famiglia batterica.