Generalized linear models for dependent frequency and severity of insurance claims

Traditionally, claim counts and amounts are assumed to be independent in non-life insurance. This paper explores how this often unwarranted assumption can be relaxed in a simple way while incorporating rating factors into the model. The approach consists of fitting generalized linear models to the marginal frequency and the conditional severity components of the total claim cost; dependence between them is induced by treating the number of claims as a covariate in the model for the average claim size. In addition to being easy to implement, this modeling strategy has the advantage that when Poisson counts are assumed together with a log-link for the conditional severity model, the resulting pure premium is the product of a marginal mean frequency, a modified marginal mean severity, and an easily interpreted correction term that reflects the dependence. The approach is illustrated through simulations and applied to a Canadian automobile insurance dataset.     

Bayesian cohort models for general cohort table analyses

Summary New Bayesian cohort models designed to resolve the identification problem in cohort analysis are proposed in this paper. At first, the basic cohort model which represents the statistical structure of time-series social survey data in terms of age, period and cohort effects is explained. The logit cohort model for qualitative data from a binomial distribution and the normal-type cohort model for quantitative data from a normal distribution are considered as two special cases of the basic model. In order to overcome the identification problem in cohort analysis, a Bayesian approach is adopted, based on the assumption that the effect parameters change gradually. A Bayesian information criterion ABIC is introduced for the selection of the optimal model. This approach is so flexible that both the logit and the normal-type cohort models can be made applicable, not only to standard cohort tables but also to general cohort tables in which the range of age group is not equal to the interval between periods. The practical utility of the proposed models is demonstrated by analysing two data sets from the literature on cohort analysis. The Institute of Statistical Mathematics     

The Poisson random effect model for experience ratemaking: Limitations and alternative solutions

Poisson random effect models with a shared random effect have been widely used in actuarial science for analyzing the number of claims. In particular, the random effect is a key factor in a posteriori risk classification. However, the necessity of the random effect may not be properly assessed due to the dual role of the random effect; it affects both the marginal distribution of the number of claims and the dependence among the numbers of claims obtained from an individual over time. We first show that the score test for the nullity of the variance of the shared random effect can falsely indicate significant dependence among the numbers of claims even though they are independent. To mitigate this problem, we propose to separate the dual role of the random effect by introducing additional random effects to capture the overdispersion part, which are called saturated random effects. In order to circumvent heavy computational issues by the saturated random effects, we choose a gamma distribution for the saturated random effects because it gives the closed form of marginal distribution. In fact, this choice leads to the negative binomial random effect model that has been widely used for the analysis of frequency data. We show that safer conclusions about the a posteriori risk classification can be made based on the negative binomial mixed model under various situations. We also derive the score test as a sufficient condition for the existence of the a posteriori risk classification based on the proposed model.     

On the Type I multivariate zero-truncated hurdle model with applications in health insurance

In the general insurance modeling literature, there has been a lot of work based on univariate zero-truncated models, but little has been done in the multivariate zero-truncation cases, for instance a line of insurance business with various classes of policies. There are three types of zero-truncation in the multivariate setting: only records with all zeros are missing, zero counts for one or some classes are missing, or zeros are completely missing for all classes. In this paper, we focus on the first case, the so-called Type I zero-truncation, and a new multivariate zero-truncated hurdle model is developed to study it. The key idea of developing such a model is to identify a stochastic representation for the underlying random variables, which enables us to use the EM algorithm to simplify the estimation procedure. This model is used to analyze a health insurance claims dataset that contains claim counts from different categories of claims without common zero observations.     

Drawing inferences from logit models for panel data

Logit models have been widely used in marketing to predict brand choice and to make inference about the impact of marketing mix variables on these choices. Most researchers have followed the pioneering example of Guadagni and Little, building choice models and drawing inference conditional on the assumption that the logit model is the correct specification for household purchase behaviour. To the extent that logit models fail to adequately describe household purchase behaviour, statistical inferences from them may be flawed. More importantly, marketing decisions based on these models may be incorrect. This research applies White's robust inference method to logit brand choice models. The method does not impose the restrictive assumption that the assumed logit model specification be true. A sandwich estimator of the covariance ‘corrected’ for possible mis‐specification is the basis for inference about logit model parameters. An important feature of this method is that it yields correct standard errors for the marketing mix parameter estimates even if the assumed logit model specification is not correct. Empirical examples include using household panel data sets from three different product categories to estimate logit models of brand choice. The standard errors obtained using traditional methods are compared with those obtained by White's robust method. The findings illustrate that incorrectly assuming the logit model to be true typically yields standard errors which are biased downward by 10–40 per cent. Conditions under which the bias is particularly severe are explored. Under these conditions, the robust approach is recommended. Copyright © 2000 John Wiley & Sons, Ltd.     

Identification of graphical models for nonignorable nonresponse of binary outcomes in longitudinal studies

In this paper, we use directed acyclic graphs (DAGs) with temporal structure to describe models of nonignorable nonresponse mechanisms for binary outcomes in longitudinal studies, and we discuss identification of these models under an assumption that the sequence of variables has the first-order Markov dependence, that is, the future variables are independent of the past variables conditional on the present variables. We give a stepwise approach for checking identifiability of DAG models. For an unidentifiable model, we propose adding completely observed variables such that this model becomes identifiable.     

Multiway Dependence in Exponential Family Conditional Distributions

Conditionally specified statistical models are frequently constructed from one-parameter exponential family conditional distributions. One way to formulate such a model is to specify the dependence structure among random variables through the use of a Markov random field (MRF). A common assumption on the Gibbsian form of the MRF model is that dependence is expressed only through pairs of random variables, which we refer to as the “pairwise-only dependence” assumption. Based on this assumption, J. Besag (1974, J. Roy. Statist. Soc. Ser. B36, 192–225) formulated exponential family “auto-models” and showed the form that one-parameter exponential family conditional densities must take in such models. We extend these results by relaxing the pairwise-only dependence assumption, and we give a necessary form that one-parameter exponential family conditional densities must take under more general conditions of multiway dependence. Data on the spatial distribution of the European corn borer larvae are fitted using a model with Bernoulli conditional distributions and several dependence structures, including pairwise-only, three-way, and four-way dependencies.     

Modeling special-day effects for forecasting intraday electricity demand

We propose and apply a novel approach for modeling special-day effects to predict electricity demand in Korea. Notably, we model special-day effects on an hourly rather than a daily basis. Hourly specified predictor variables are implemented in the regression model with a seasonal autoregressive moving average (SARMA) type error structure in order to efficiently reflect the special-day effects. The interaction terms between the hour-of-day effects and the hourly based special-day effects are also included to capture the unique intraday patterns of special days more accurately. The multiplicative SARMA mechanism is employed in order to identify the double seasonal cycles, namely, the intraday effect and the intraweek effect. The forecast results of the suggested model are evaluated by comparing them with those of various benchmark models for the following year. The empirical results indicate that the suggested model outperforms the benchmark models for both special- and non-special day predictions.     

Tree-Structured Clustering in Fixed Effects Models

Fixed effects models are very flexible because they do not make assumptions on the distribution of effects and can also be used if the heterogeneity component is correlated with explanatory variables. A disadvantage is the large number of effects that have to be estimated. A recursive partitioning (or tree based) method is proposed that identifies clusters of units that share the same effect. The approach reduces the number of parameters to be estimated and is useful in particular if one is interested in identifying clusters with the same effect on a response variable. It is shown that the method performs well and outperforms competitors like the finite mixture model in particular if the heterogeneity component is correlated with explanatory variables. In two applications the usefulness of the approach to identify clusters that share the same effect is illustrated. Supplementary materials for this article are available online.     

A revisit of spinning disk models. Part II: linear transverse vibrations

In this paper, two recently derived linear models for the transverse vibrations of a spinning plate are considered. The disk is modelled as a pure plate with no membrane effects. Furthermore, the effect of the rotary inertia of the plate is taken into account. The first model is based on the assumption of linear (Kirchhoff) strains. The second model is based on the assumption of non-linear (von Karman) strains. The merits of both models are considered and their predictions are compared with those of the traditional linear model.     

An application of a hybrid approach to modeling a flexible manufacturing system

Both analytic and simulation models were used to analyze the capabilities and requirements of an automated circuit card manufacturing system. Analytic models were used to determine the sensitivity of the measures of effectiveness (MOEs) to various design parameters. This analysis gave approximate results and bounded the range of input parameters for the simulation model. A detailed simulation model was required for use during both the design and production phases of the project. This simulation model incorporated only those variables to which the MOEs are most sensitive, and provided additional features to observe system behavior. The benefits and appropriate uses for each class of models are discussed.     

复合索赔次数模型与混合索赔次数模型中的若干结果

本文利用完整描述方法研究复合索赔次数模型与混合索赔次数模型中总索赔次数的概率分布 ,得到了十余例典型索赔次数模型的相关结果 ,这些结果推广了文献 [1]、[2 ]、[6 ]的有关结论。     

Discovering model structure for partially linear models

Partially linear models (PLMs) have been widely used in statistical modeling, where prior knowledge is often required on which variables have linear or nonlinear effects in the PLMs. In this paper, we propose a model-free structure selection method for the PLMs, which aims to discover the model structure in the PLMs through automatically identifying variables that have linear or nonlinear effects on the response. The proposed method is formulated in a framework of gradient learning, equipped with a flexible reproducing kernel Hilbert space. The resultant optimization task is solved by an efficient proximal gradient descent algorithm. More importantly, the asymptotic estimation and selection consistencies of the proposed method are established without specifying any explicit model assumption, which assure that the true model structure in the PLMs can be correctly identified with high probability. The effectiveness of the proposed method is also supported by a variety of simulated and real-life examples.

     

Use of the Probability Integral Transformation to Fit Nonlinear Mixed-Effects Models With Nonnormal Random Effects

This article describes a simple computational method for obtaining the maximum likelihood estimates (MLE) in nonlinear mixed-effects models when the random effects are assumed to have a nonnormal distribution. Many computer programs for fitting nonlinear mixed-effects models, such as PROC NLMIXED in SAS, require that the random effects have a normal distribution. However, there is often interest in either fitting models with nonnormal random effects or assessing the sensitivity of inferences to departures from the normality assumption for the random effects. When the random effects are assumed to have a nonnormal distribution, we show how the probability integral transform can be used, in conjunction with standard statistical software for fitting nonlinear mixed-effects models (e.g., PROC NLMIXED in SAS), to obtain the MLEs. Specifically, the probability integral transform is used to transform a normal random effect to a nonnormal random effect. The method is illustrated using a gamma frailty model for clustered survival data and a beta-binomial model for clustered binary data. Finally, the results of a simulation study, examining the impact of misspecification of the distribution of the random effects, are presented.     

高维异方差高斯混合罚模型聚类

罚模型聚类实现了在聚类过程中精简变量的目标,同时如何识别聚类的有效变量成了一个新的问题.在这个问题上,已有的研究有成对罚模型,模型处理了各类数据同方差的情况.考察了异方差情况下的变量选择问题,针对异方差数据提出了两种新的模型,并给出模型的解和算法.模拟数据分析结果表明,异方差数据上两个新模型都有更好的表现.     

基于一般学习效应下的单机排序问题

近来具有学习效应的机器排序问题收到广泛的关注.对于机器排序中工件的实际加工来说,与工件加工位置有关的学习模型更具有现实性.本文研究了工件加工位置和与已经加工过的工件之和有关的一般学习效应模型.首先证明文献中与位置和已经加工过的工件加工时间之和有关的学习模型是本模型的特殊情形.其次对于单机排序问题我们提出一般解法.     

On modeling evaporation

We develop a theoretical model for evaporation of a pure liquid drop on a thermally conductive solid substrate. We discuss a variety of effects regarding evaporation regime, the state of the liquid/gas interface and the content of gas phase. Then, we further consider two models: the one resulting from the one-sided non-equilibrium assumption and the other that assumes diffusion-limited regime and equilibrium at the liquid/gas interface. A single governing equation for the evolution of drop thickness is derived for both models. We show that although the model predicts qualitatively different temperature along liquid/gas and liquid/solid interface, the dynamics of the drops is almost the same.      

Understanding and interpreting generalized ordered logit models

When outcome variables are ordinal rather than continuous, the ordered logit model, aka the proportional odds model (ologit/po), is a popular analytical method. However, generalized ordered logit/partial proportional odds models (gologit/ppo) are often a superior alternative. Gologit/ppo models can be less restrictive than proportional odds models and more parsimonious than methods that ignore the ordering of categories altogether. However, the use of gologit/ppo models has itself been problematic or at least sub-optimal. Researchers typically note that such models fit better but fail to explain why the ordered logit model was inadequate or the substantive insights gained by using the gologit alternative. This paper uses both hypothetical examples and data from the 2012 European Social Survey to address these shortcomings.     

Modeling frequency and severity of claims with the zero-inflated generalized cluster-weighted models

To facilitate applications in general insurance, some extensions are proposed to cluster-weighted models (CWMs). First, we extend CWMs to have generalized cluster-weighted models (GCWMs) by allowing modeling of non-Gaussian distribution of the continuous covariates, as they frequently occur in insurance practice. Secondly, we introduce a zero-inflated extension of GCWM (ZI-GCWM) for modeling insurance claims data with excess zeros coming from heterogeneous sources. Additionally, we give two expectation–optimization (EM) algorithms for parameter estimation given in the proposed models. An appropriate simulation study shows that, for various settings and in contrast to the existing mixture-based approaches, both extended models perform well. Finally, a real data set based on French auto-mobile policies is used to illustrate the application of the proposed extensions.     

The age pattern of transitory mortality jumps and its impact on the pricing of catastrophic mortality bonds

To value catastrophic mortality bonds, a number of stochastic mortality models with transitory jump effects have been proposed. Rather than modeling the age pattern of jump effects explicitly, most of the existing models assume that the distributions of jump effects and general mortality improvements across ages are identical. Nevertheless, this assumption does not seem to be in line with what we observe from historical data. In this paper, we address this problem by introducing a Lee–Carter variant that captures the age pattern of mortality jumps by a distinct collection of parameters. The model variant is then further generalized to permit the age pattern of jump effects to vary randomly. We illustrate the two proposed models with mortality data from the United States and English and Welsh populations, and use them to value hypothetical mortality bonds with similar specifications to the Atlas IX Capital Class B note that was launched in 2013. It is found that the features we consider have a significant impact on the estimated prices.