Kernel Estimation of Rate Function for Recurrent Event Data

Abstract.  Recurrent event data are largely characterized by the rate function but smoothing techniques for estimating the rate function have never been rigorously developed or studied in statistical literature. This paper considers the moment and least squares methods for estimating the rate function from recurrent event data. With an independent censoring assumption on the recurrent event process, we study statistical properties of the proposed estimators and propose bootstrap procedures for the bandwidth selection and for the approximation of confidence intervals in the estimation of the occurrence rate function. It is identified that the moment method without resmoothing via a smaller bandwidth will produce a curve with nicks occurring at the censoring times, whereas there is no such problem with the least squares method. Furthermore, the asymptotic variance of the least squares estimator is shown to be smaller under regularity conditions. However, in the implementation of the bootstrap procedures, the moment method is computationally more efficient than the least squares method because the former approach uses condensed bootstrap data. The performance of the proposed procedures is studied through Monte Carlo simulations and an epidemiological example on intravenous drug users.     

GMM tests for the Katz family of distributions

Generalized method of moments (GMM) is used to develop tests for discriminating discrete distributions among the two-parameter family of Katz distributions. Relationships involving moments are exploited to obtain identifying and over-identifying restrictions. The asymptotic relative efficiencies of tests based on GMM are analyzed using the local power approach and the approximate Bahadur efficiency. The paper also gives results of Monte Carlo experiments designed to check the validity of the theoretical findings and to shed light on the small sample properties of the proposed tests. Extensions of the results to compound Poisson alternative hypotheses are discussed.     

Creating a reference to a complex emergency situation using time series methods: war in Guinea-Bissau 1998-1999

Impacts of complex emergencies or relief interventions have often been evaluated by absolute mortality compared to international standardized mortality rates. A better evaluation would be to compare with local baseline mortality of the affected populations. A projection of population-based survival data into time of emergency or intervention based on information from before the emergency may create a local baseline reference. We find a log-transformed Gaussian time series model where standard errors of the estimated rates are included in the variance to have the best forecasting capacity. However, if time-at-risk during the forecasted period is known then forecasting might be done using a Poisson time series model with overdispersion. Whatever, the standard error of the estimated rates must be included in the variance of the model either in an additive form in a Gaussian model or in a multiplicative form by overdispersion in a Poisson model. Data on which the forecasting is based must be modelled carefully concerning not only calendar-time trends but also periods with excessive frequency of events (epidemics) and seasonal variations to eliminate residual autocorrelation and to make a proper reference for comparison, reflecting changes over time during the emergency. Hence, when modelled properly it is possible to predict a reference to an emergency-affected population based on local conditions. We predicted childhood mortality during the war in Guinea-Bissau 1998-1999. We found an increased mortality in the first half-year of the war and a mortality corresponding to the expected one in the last half-year of the war.     

Statistik für bivariate gemischte Poisson–Prozesse am Beispiel der Kraftfahrthaftpflichtversicherung

Zusammenfassung: In diesem Artikel wird der Weg von einem univariaten gemischten Poisson–Prozess, der in vielen Bereichen zum Z?hlen von Ereignissen benutzt wird, zu einem bivariaten gemischten Poisson–Prozess aufgezeigt. Dazu werden einige Eigenschaften des bivariaten Prozesses angegeben. Im zweiten Teil der Arbeit wird gezeigt, wie mit Hilfe dieses Prozesses der übergang von einem herk?mmlichen Bonus–Malus–System in der Kraftfahrthaftpflichtversicherung zu einem Bonus–Malus–System mit Berücksichtigung der Schadenart beschritten werden kann. Dazu wird zuerst eine Modellprüfung der gegebenen Daten vorgenommen und sodann werden für verschiedene mischende Verteilungen die Verteilungsparameter gesch?tzt und Nettopr?mien angegeben sowie die Prognosegenauigkeit getestet.

---

Summary: In this paper we show that the model of the bivariate mixed Poisson process arises in a natural way from the univariate mixed Poisson process, which is used in several areas for counting certain events. Furthermore we state some properties of the bivariate process. In the second part of the paper we illustrate how by means of the bivariate mixed Poisson process a bonus–malus system handling different types of accidents can be derived from the classical bonus–malus system in third–party liability insurance. To this end we first check the model on the given data and then estimate distribution parameters and compute net premiums for different mixing distributions as well as test the prediction probabilities.

* Vortrag am Dresdner Forum zur Versicherungsmathematik: Tarifierung in Erst- und Rückversicherung am 25. Juni 2004. Für die Unterstützung zu dieser Arbeit m?chte der Autor Lothar Partzsch, Klaus D. Schmidt (beide Dresden) und Friedemann Spies (München) recht herzlich danken.     

Contre-exemple au théorème de P.C. Consul sur la factorisation des distributions de Poisson généralisées

This note exhibits two independent random variables on integers, X₁ and X₂, such that neither X₁ nor X₂ has a generalized Poisson distribution, but X₁ + X₂ has. This contradicts statements made by Professor Consul in his recent book.     

A note on LeCam’s bound for the distance between the Poisson binomial and the Poisson distribution

n possibly different success probabilities p ₁, p ₂, ..., p _n is frequently approximated by a Poisson distribution with parameter λ = p ₁ + p ₂ + ... + p _n . LeCam's bound p ² ₁ + p ² ₂ + ... + p _n ² for the total variation distance between both distributions is particularly useful provided the success probabilities are small. The paper presents an improved version of LeCam's bound if a generalized d-dimensional Poisson binomial distribution is to be approximated by a compound Poisson distribution. Received: May 10, 2000; revised version: January 15, 2001     

CreditRisk+模型采用Poisson分布所产生的经济资本计量误差分析

本文对CreditRisk+模型采用Poisson分布近似债务人违约事件分布这一关键步骤进行系统研究。首先分析了Poisson分布在CreditRisk+模型中的作用,并从理论上证明了采用Poisson分布作为债务人违约事件分布的近似,会导致CreditRisk+模型计算出来的经济资本高估贷款组合的实际风险水平;然后以债务人违约事件服从两点分布并采用蒙特卡罗模拟计算出来的经济资本为参照值,对债务人违约概率的大小与这一近似所引起的经济资本计量误差率进行了敏感性试验,发现为将这一近似所引起的误差率控制在10%的范围内,债务人违约概率的取值不应超过0.2。     

Hostage taking: Understanding terrorism event dynamics

This paper employs advanced time series methods to identify the dynamic properties of three hostage taking series. The immediate and long run multipliers of three covariates—successful past negotiations, violent ends, and deaths—are identified. Each hostage series responds differently to the covariates. Past concessions have the strongest impact on generating future kidnapping events, supporting the conventional wisdom to abide by a stated no-concession policy. Each hostage series has different changepoints caused by a variety of circumstances. Skyjackings and kidnappings are negatively correlated, while skyjackings and other hostage events are positively correlated. Policy recommendations are offered.     

The empirical Bayes estimators of the parameter of the Poisson distribution with a conjugate gamma prior under Stein's loss function

For the hierarchical Poisson and gamma model, we calculate the Bayes posterior estimator of the parameter of the Poisson distribution under Stein's loss function which penalizes gross overestimation and gross underestimation equally and the corresponding Posterior Expected Stein's Loss (PESL). We also obtain the Bayes posterior estimator of the parameter under the squared error loss and the corresponding PESL. Moreover, we obtain the empirical Bayes estimators of the parameter of the Poisson distribution with a conjugate gamma prior by two methods. In numerical simulations, we have illustrated: The two inequalities of the Bayes posterior estimators and the PESLs; the moment estimators and the Maximum Likelihood Estimators (MLEs) are consistent estimators of the hyperparameters; the goodness-of-fit of the model to the simulated data. The numerical results indicate that the MLEs are better than the moment estimators when estimating the hyperparameters. Finally, we exploit the attendance data on 314 high school juniors from two urban high schools to illustrate our theoretical studies.     

A new stratified three-stage unrelated randomized response model for estimating a rare sensitive attribute based on the Poisson distribution

This article suggests an efficient method of estimating a rare sensitive attribute which is assumed following Poisson distribution by using three-stage unrelated randomized response model instead of the Land et al. model (2011 Land, M., S. Singh, and S. A. Sedory. 2011. Estimation of a rare sensitive attribute using poisson distribution. Statistics 46 (3):351–60. doi:10.1080/02331888.2010.524300.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]) when the population consists of some different sized clusters and clusters selected by probability proportional to size(:pps) sampling. A rare sensitive parameter is estimated by using pps sampling and equal probability two-stage sampling when the parameter of a rare unrelated attribute is assumed to be known and unknown.

We extend this method to the case of stratified population by applying stratified pps sampling and stratified equal probability two-stage sampling. An empirical study is carried out to show the efficiency of the two proposed methods when the parameter of a rare unrelated attribute is assumed to be known and unknown.