Finite time ruin probability and structural density properties in the presence of dependence in insurance risk model

In the present paper, we consider the classical compound Poisson risk model with dependence between claim sizes and claim inter-arrival time. We attempt to analyze the approximation of finite time ruin probability. The finite time ruin probabilities are plotted for fixed threshold value associated to the claim inter-arrival time and also for fixed dependence parameter in Nelsen (2006) copula separately. Additionally, a general form for joint density of the interclaim times and claim sizes is considered. With respect to the classical Gerber-Shiu's (1998) function, first some structural density properties of dependent collective risk model is obtained. Then the ladder height probability density function of claim sizes is computed and the dependency structure investigated for Erlang interclaim time. As the application, some dependent models of the interclaim times and claim sizes are studied.     

Uniform Asymptotic Estimates for Ruin Probabilities with Exponential Lévy Process Investment Returns and Two-sided Linear Heavy-tailed Claims

This paper investigates the ruin probabilities of a renewal risk model with stochastic investment returns and dependent claim sizes. The investment is described as a portfolio of one risk-free asset and one risky asset whose price process is an exponential Lévy process. The claim sizes are assumed to follow a two-sided linear process with independent and identically distributed step sizes. When the step-size distribution is heavy tailed, the paper establishes some uniform asymptotic formulas of ruin probabilities.     

On Erlang(2) Risk Process Perturbed by Diffusion

ABSTRACT

In this article, we consider an Erlang(2) risk process perturbed by diffusion. From the extreme value distribution of Brownian motion with drift and the renewal theory, we show that the survival probability satisfies an integral equation. We then give the bounds for the ultimate ruin probability and the ruin probability caused by claim. By introducing a random walk associated with the proposed risk process, we define an adjustment-coefficient. The relation between the adjustment-coefficient and the bound is given and the Lundberg-type inequality for the bound is obtained. Also, a formula of Pollaczek–Khinchin type for the bound is derived. Using these results, the bound can be calculated when claim sizes are exponentially distributed.     

A surplus process involving a compound Poisson counting process and applications

Abstract

In this paper, we introduce a surplus process involving a compound Poisson counting process, which is a generalization of the classical ruin model where the claim-counting process is a homogeneous Poisson process. The incentive is to model batch arrival of claims using a counting process that is based on a compound distribution. This reduces the difficulty of modeling claim amounts and is consistent with industrial data. Recursive formula, some properties and relevant main ruin theory results are provided. Further, we consider applications involving zero-truncated negative binomial and zero-truncated binomial batch arrivals when the claim amounts follow exponential or Erlang distribution.     

Bayesian sample size determination for phase IIA clinical trials using historical data and semi‐parametric prior's elicitation

The Simon's two‐stage design is the most commonly applied among multi‐stage designs in phase IIA clinical trials. It combines the sample sizes at the two stages in order to minimize either the expected or the maximum sample size. When the uncertainty about pre‐trial beliefs on the expected or desired response rate is high, a Bayesian alternative should be considered since it allows to deal with the entire distribution of the parameter of interest in a more natural way. In this setting, a crucial issue is how to construct a distribution from the available summaries to use as a clinical prior in a Bayesian design. In this work, we explore the Bayesian counterparts of the Simon's two‐stage design based on the predictive version of the single threshold design. This design requires specifying two prior distributions: the analysis prior, which is used to compute the posterior probabilities, and the design prior, which is employed to obtain the prior predictive distribution. While the usual approach is to build beta priors for carrying out a conjugate analysis, we derived both the analysis and the design distributions through linear combinations of B‐splines. The motivating example is the planning of the phase IIA two‐stage trial on anti‐HER2 DNA vaccine in breast cancer, where initial beliefs formed from elicited experts' opinions and historical data showed a high level of uncertainty. In a sample size determination problem, the impact of different priors is evaluated.     

The ruin probabilities of a discrete time risk model with one-sided linear claim sizes and dependent risks

This article investigates the ruin probabilities of a discrete time risk model with dependent claim sizes and dependent relation between insurance risks and financial risks. The risk-free and risky investments of an insurer lead to stochastic discount factors {θ_n}_{n ? 1}. The claim sizes are assumed to follow a one-sided linear process with independent and identically distributed (i.i.d.) innovations {?_n}_{n ? 1}. The i.i.d. random pairs {(?_n, θ_n)}_{n ? 1} follow a common bivariate Sarmanov-dependent distribution. When the common distribution of the innovations is heavy tailed, we establish some asymptotic estimates for the ruin probabilities of this discrete time risk model.     

Bayesiam simultaneous estimation for several multinomial distributions

A Bayesian method is proposed for estimating the cell probabilities of several multinomial distributions. Parameters of different distributions are taken to be a priori exchangeable. The prior specification is based upon mixtures of a hierarchical distribution, referred to as the multivariate “Dirichlet-Dirichlet” distribution. The analysis is facilitated by a multinomial approximation relating to the multinomial-Dirichlet distribution. The posterior estimates depend upon measures of entropy for the various distributions and shrink the individual observed proportions towards values obtained by pooling the data across the distributions. As well as incorporating prior information they are particularly useful when some of the cell frequencies are zero. We use them to investigate a numerical classification of males of various vocations, according to cause of death.     

Asymptotic ruin probabilities of a two-dimensional renewal risk model with dependent inter-arrival times

Abstract

This article mainly considers the uniform asymptotics for the finite-time ruin probabilities of a two-dimensional renewal risk model with heavy-tailed claims. In this model, the two claim-number processes are arbitrarily dependent and each of them is generated by widely orthant dependent claim inter-arrival times. Two types of ruin are studied and for each type of ruin, an asymptotic formula for the finite-time ruin probability is established. These formulae possess a certain uniformity feature in the time horizon.     

Bayesian analysis of mixture modelling using the multivariate t distribution

A finite mixture model using the multivariate t distribution has been shown as a robust extension of normal mixtures. In this paper, we present a Bayesian approach for inference about parameters of t-mixture models. The specifications of prior distributions are weakly informative to avoid causing nonintegrable posterior distributions. We present two efficient EM-type algorithms for computing the joint posterior mode with the observed data and an incomplete future vector as the sample. Markov chain Monte Carlo sampling schemes are also developed to obtain the target posterior distribution of parameters. The advantages of Bayesian approach over the maximum likelihood method are demonstrated via a set of real data.     

Upper bounds for ruin probabilities under model uncertainty

Abstract

In this paper, we investigate some ruin problems for risk models that contain uncertainties on both claim frequency and claim size distribution. The problems naturally lead to the evaluation of ruin probabilities under the so-called G-expectation framework. We assume that the risk process is described as a class of G-compound Poisson process, a special case of the G-Lévy process. By using the exponential martingale approach, we obtain the upper bounds for the two-sided ruin probability as well as the ruin probability involving investment. Furthermore, we derive the optimal investment strategy under the criterion of minimizing this upper bound. Finally, we conclude that the upper bound in the case with investment is less than or equal to the case without investment.     

Bayesian estimation for the M/G/1 queue using a phase-type approximation

This article deals with Bayesian inference and prediction for M/G/1 queueing systems. The general service time density is approximated with a class of Erlang mixtures which are phase-type distributions. Given this phase-type approximation, an explicit evaluation of measures such as the stationary queue size, waiting time and busy period distributions can be obtained. Given arrival and service data, a Bayesian procedure based on reversible jump Markov Chain Monte Carlo methods is proposed to estimate system parameters and predictive distributions.     

A BAYESIAN SIMULATION APPROACH TO INFERENCE ON A MULTI‐STATE LATENT FACTOR INTENSITY MODEL

The influence of economic conditions on the movement of a variable between states (for example a change in credit rating from A to B) can be modelled using a multi‐state latent factor intensity framework. Estimation of this type of model is, however, not straightforward, as transition probabilities are involved and the model contains a few highly analytically intractable distributions. In this paper, a Bayesian approach is adopted to manage the distributions. The innovation in the sampling algorithm used to obtain the posterior distributions of the model parameters includes a particle filter step and a Metropolis–Hastings step within a Gibbs sampler. The feasibility and accuracy of the proposed sampling algorithm is supported with a few simulated examples. The paper contains an application concerning what caused 1049 firms to change their credit ratings over a span of ten years.     

A Bayesian Hierarchical Model for Multiple Comparisons in Mixed Models

We propose a Bayesian hierarchical model for multiple comparisons in mixed models where the repeated measures on subjects are described with the subject random effects. The model facilitates inferences in parameterizing the successive differences of the population means, and for them, we choose independent prior distributions that are mixtures of a normal distribution and a discrete distribution with its entire mass at zero. For the other parameters, we choose conjugate or vague priors. The performance of the proposed hierarchical model is investigated in the simulated and two real data sets, and the results illustrate that the proposed hierarchical model can effectively conduct a global test and pairwise comparisons using the posterior probability that any two means are equal. A simulation study is performed to analyze the type I error rate, the familywise error rate, and the test power. The Gibbs sampler procedure is used to estimate the parameters and to calculate the posterior probabilities.     

Scale mixtures log-Birnbaum–Saunders regression models with censored data: a Bayesian approach

The main objective of this paper is to develop a full Bayesian analysis for the Birnbaum–Saunders (BS) regression model based on scale mixtures of the normal (SMN) distribution with right-censored survival data. The BS distributions based on SMN models are a very general approach for analysing lifetime data, which has as special cases the Student-t-BS, slash-BS and the contaminated normal-BS distributions, being a flexible alternative to the use of the corresponding BS distribution or any other well-known compatible model, such as the log-normal distribution. A Gibbs sample algorithm with Metropolis–Hastings algorithm is used to obtain the Bayesian estimates of the parameters. Moreover, some discussions on the model selection to compare the fitted models are given and case-deletion influence diagnostics are developed for the joint posterior distribution based on the Kullback–Leibler divergence. The newly developed procedures are illustrated on a real data set previously analysed under BS regression models.     

Diffusion approximations for insurance risk processes

We consider a risk-reserve process for an insurance company where premium income and the claim sum process are modeled as a renewal reward processes. Moreover, dividends are paid out according to a barrier rule. The aim of the article is to establish a diffusion approximation of this model and to compute ruin probabilities (in finite and in infinite time) and other relevant statistics approximately using the limiting diffusion process. We also demonstrate that, under special circumstances, there exists a stationary distribution for the limiting diffusion.     

Subexponential Distributions — Large Deviations with Applications to Insurance and Queueing Models

This paper presents a fine large‐deviations theory for heavy‐tailed distributions whose tails are heavier than exp(?√t and have finite second moment. Asymptotics for first passage times are derived. The results are applied to estimate the finite time ruin probabilities in insurance as well as the busy period in a GI/G/1 queueing model.     

Bayesian analysis of a multiple-recapture model

In the Bayesian analysis of a multiple-recapture census, different diffuse prior distributions can lead to markedly different inferences about the population size N. Through consideration of the Fisher information matrix it is shown that the number of captures in each sample typically provides little information about N. This suggests that if there is no prior information about capture probabilities, then knowledge of just the sample sizes and not the number of recaptures should leave the distribution of Nunchanged. A prior model that has this property is identified and the posterior distribution is examined. In particular, asymptotic estimates of the posterior mean and variance are derived. Differences between Bayesian and classical point and interval estimators are illustrated through examples.     

Forecasting future values of changing sequences

Sequences of independent random variables are observed and on the basis of these observations future values of the process are forecast. The Bayesian predictive density of k future observations for normal, exponential, and binomial sequences which change exactly once are analyzed for several cases. It is seen that the Bayesian predictive densities are mixtures of standard probability distributions. For example, with normal sequences the Bayesian predictive density is a mixture of either normal or t-distributions, depending on whether or not the common variance is known. The mixing probabilities are the same as those occurring in the corresponding posterior distribution of the mean(s) of the sequence. The predictive mass function of the number of future successes that will occur in a changing Bernoulli sequence is computed and point and interval predictors are illustrated.     

Skew-probit measurement error models

In this paper we extend the structural probit measurement error model by considering the unobserved covariate to follow a skew-normal distribution. The new model is termed the structural skew-normal probit model. As in the normal case, the likelihood function is obtained analytically, and can be maximized by using existing statistical software. A Bayesian approach using Markov chain Monte Carlo techniques for generating from the posterior distributions is also developed. A simulation study demonstrates the usefulness of the approach in avoiding attenuation which arises with the naive procedure. Moreover, a comparison of predicted and true success probabilities indicates that it seems to be more efficient to use the skew probit model when the distribution of the covariate (predictor) is skew. An application to a real data set is also provided.     

On a discrete interaction risk model with delayed claims and stochastic incomes under random discount rates

In this paper, we study a discrete interaction risk model with delayed claims and stochastic incomes in the framework of the compound binomial model. A generalized Gerber-Shiu discounted penalty function is proposed to analyse this risk model in which the interest rates follow a Markov chain with finite state space. We derive an explicit expression for the generating function of this Gerber-Shiu discounted penalty function. Furthermore, we derive a recursive formula and a defective renewal equation for the original Gerber-Shiu discounted penalty function. As an application, the joint distributions of the surplus one period prior to ruin and the deficit at ruin, as well as the probabilities of ruin are obtained. Finally, some numerical illustrations from a specific example are also given.