On a bounded bimodal two-sided distribution fitted to the Old-Faithful geyser data

In this paper, we shall develop a novel family of bimodal univariate distributions (also allowing for unimodal shapes) and demonstrate its use utilizing the well-known and almost classical data set involving durations and waiting times of eruptions of the Old-Faithful geyser in Yellowstone park. Specifically, we shall analyze the Old-Faithful data set with 272 data points provided in Dekking et al. [3]. In the process, we develop a bivariate distribution using a copula technique and compare its fit to a mixture of bivariate normal distributions also fitted to the same bivariate data set. We believe the fit-analysis and comparison is primarily illustrative from an educational perspective for distribution theory modelers, since in the process a variety of statistical techniques are demonstrated. We do not claim one model as preferred over the other.     

A strategy for constructing multivariate distributions

Given a random vector (X₁,…, X_n) for which the univariate and bivariate marginal distributions belong to some specified families of distributions, we present a procedure for constructing families of multivariate distributions with the specified univariate and bivariate margins. Some general properties of the resulting families of multivariate distributions are reviewed. This procedure is illustrated by generalizing the bivariate Plackett (1965) and Clayton (1978) distributions to three dimensions. In addition to providing rich families of models for data analysis, this method of construction provides a convenient way of simulating observations from multivariate distributions with specific types of univariate and bivariate marginal distributions. A general algorithm for simulating random observations from these families of multivariate distributions is presented     

A bivariate extension of the beta generated distribution derived from copulas

In this paper, we introduce a new class of bivariate distributions whose marginals are beta-generated distributions. Copulas are employed to construct this bivariate extension of the beta-generated distributions. It is shown that when Archimedean copulas and convex beta generators are used in generating bivariate distributions, the copulas of the resulting distributions also belong to the Archimedean family. The dependence of the proposed bivariate distributions is examined. Simulation results for beta generators and an application to financial risk management are presented.     

On the probability of ruin in a continuous risk model with two types of delayed claims

ABSTRACT

This article studies a risk model involving one type of main claims and two types of by-claims, which is an extension of the general risk model with delayed claims. We suppose that every main claim may not induce any by-claims or may induce one by-claim belonging to one of the two types of by-claims with a certain probability. In addition, assume that the by-claim and its associated main claim may occur at the same time and that the occurrence of the by-claim may be delayed. An integro-differential equation system for survival probabilities is derived by using two auxiliary risk models. The expression of the survival probability is obtained by applying Laplace transforms and Rouché theorem. Furthermore, we provide a method for solving the survival probability when the two by-claim amounts satisfy different exponential distributions. As a special case, an explicit expression of survival probability is given when all the claim amounts obey the same exponential distribution. Finally, numerical results are provided to examine the proposed method.     

A bivariate geometric distribution allowing for positive or negative correlation

In this paper, we propose a new bivariate geometric model, derived by linking two univariate geometric distributions through a specific copula function, allowing for positive and negative correlations. Some properties of this joint distribution are presented and discussed, with particular reference to attainable correlations, conditional distributions, reliability concepts, and parameter estimation. A Monte Carlo simulation study empirically evaluates and compares the performance of the proposed estimators in terms of bias and standard error. Finally, in order to demonstrate its usefulness, the model is applied to a real data set.     

On the Gerber–Shiu discounted penalty function in a risk model with delayed claims

In this paper, we consider an extension to the continuous time risk model for which the occurrence of the claim may be delayed and the time of delay for the claim is assumed to be random. Two types of dependent claims, main claims and by-claims, are defined, where every by-claim is induced by the main claim. The time of occurrence of a by-claim is later than that of its associate main claim and the time of delay for the occurrence of a by-claim is random. An integro-differential equations system for the Gerber–Shiu discounted penalty function is established using the auxiliary risk models. Both the system of Laplace transforms of the Gerber–Shiu discounted penalty functions and the Gerber–Shiu discounted penalty functions with zero initial surplus are obtained. From Lagrange interpolating theorem, we prove that the Gerber–Shiu discounted penalty function satisfies a defective renewal equation. Exact representation for the solution of this equation is derived through an associated compound geometric distribution. Finally, examples are given with claim sizes that have exponential and a mixture of exponential distributions.     

Flexible Bivariate INAR(1) Processes Using Copulas

Multivariate count time series data occur in many different disciplines. The class of INteger-valued AutoRegressive (INAR) processes has the great advantage to consider explicitly both the discreteness and autocorrelation characterizing this type of data. Moreover, extensions of the simple INAR(1) model to the multi-dimensional space make it possible to model more than one series simultaneously. However, existing models do not offer great flexibility for dependence modelling, allowing only for positive correlation. In this work, we consider a bivariate INAR(1) (BINAR(1)) process where cross-correlation is introduced through the use of copulas for the specification of the joint distribution of the innovations. We mainly emphasize on the parametric case that arises under the assumption of Poisson marginals. Other marginal distributions are also considered. A short application on a bivariate financial count series illustrates the model.     

Non‐parametric Copula Estimation Under Bivariate Censoring

In this paper, we consider non‐parametric copula inference under bivariate censoring. Based on an estimator of the joint cumulative distribution function, we define a discrete and two smooth estimators of the copula. The construction that we propose is valid for a large range of estimators of the distribution function and therefore for a large range of bivariate censoring frameworks. Under some conditions on the tails of the distributions, the weak convergence of the corresponding copula processes is obtained in l^∞([0,1]²). We derive the uniform convergence rates of the copula density estimators deduced from our smooth copula estimators. Investigation of the practical behaviour of these estimators is performed through a simulation study and two real data applications, corresponding to different censoring settings. We use our non‐parametric estimators to define a goodness‐of‐fit procedure for parametric copula models. A new bootstrap scheme is proposed to compute the critical values.     

An absolutely continuous bivariate Topp–Leone distribution: A useful model on a bounded domain

We introduce an absolutely continuous bivariate generalization of the Topp–Leone distribution, which is a special member of the proportional reversed hazard family using a one-parameter bivariate exchangeable distribution. We show that a copula approach could also be used in defining the bivariate Topp–Leone distribution. The marginal distributions of the new bivariate distribution have also Topp–Leone distributions. We study its distributional and dependence properties. We estimate the parameters by maximum-likelihood procedure, perform a simulation study on the estimators, and apply them to a real data set. Furthermore, we give a way of generating bivariate distributions using the proposed distribution.     

Bayesian inference for a flexible class of bivariate beta distributions

Several bivariate beta distributions have been proposed in the literature. In particular, Olkin and Liu [A bivariate beta distribution. Statist Probab Lett. 2003;62(4):407–412] proposed a 3 parameter bivariate beta model which Arnold and Ng [Flexible bivariate beta distributions. J Multivariate Anal. 2011;102(8):1194–1202] extend to 5 and 8 parameter models. The 3 parameter model allows for only positive correlation, while the latter models can accommodate both positive and negative correlation. However, these come at the expense of a density that is mathematically intractable. The focus of this research is on Bayesian estimation for the 5 and 8 parameter models. Since the likelihood does not exist in closed form, we apply approximate Bayesian computation, a likelihood free approach. Simulation studies have been carried out for the 5 and 8 parameter cases under various priors and tolerance levels. We apply the 5 parameter model to a real data set by allowing the model to serve as a prior to correlated proportions of a bivariate beta binomial model. Results and comparisons are then discussed.     

A bivariate regression model with cure fraction

The use of bivariate distributions plays a fundamental role in survival and reliability studies. In this paper, we introduce a location-scale model for bivariate survival times based on the copula to model the dependence of bivariate survival data with cure fraction. We create the correlation structure between the failure times using the Clayton family of copulas, which is assumed to have any distribution. It turns out that the model becomes very flexible with respect to the choice of the marginal distributions. For the proposed model, we consider inferential procedures based on constrained parameters under maximum likelihood. We derive the appropriate matrices for assessing local influence under different perturbation schemes and present some ways to perform global influence analysis. The relevance of the approach is illustrated using a real data set and a diagnostic analysis is performed to select an appropriate model.     

On a class of bivariate mixed Sarmanov distributions

Multivariate distributions are more and more used to model the dependence encountered in many fields. However, classical multivariate distributions can be restrictive by their nature, while Sarmanov's multivariate distribution, by joining different marginals in a flexible and tractable dependence structure, often provides a valuable alternative. In this paper, we introduce some bivariate mixed Sarmanov distributions with the purpose to extend the class of bivariate Sarmanov distributions and to obtain new dependency structures. Special attention is paid to the bivariate mixed Sarmanov distribution with Poisson marginals and, in particular, to the resulting bivariate Sarmanov distributions with negative binomial and with Poisson‐inverse Gaussian marginals; these particular types of mixed distributions have possible applications in, for example modelling bivariate count data. The extension to higher dimensions is also discussed. Moreover, concerning the dependency structure, we also present some correlation formulas.     

Investigating dependence between frequency and severity via simple generalized linear models

Recently, a body of literature proposed new models relaxing a widely-used but controversial assumption of independence between claim frequency and severity in non-life insurance rate making. This paper critically reviews a generalized linear model approach, where a dependence between claim frequency and severity is introduced by treating frequency as a covariate in a regression model for severity. As an extension of this approach, we propose a dispersion model for severity. For this model, the information loss caused by using average severity rather than individual severity is examined in detail and the parameter estimators suffering from low efficiency are identified. We also provide analytical solutions for the aggregate sum to help rate making. We show that the simple functional form used in current research may not properly reflect the real underlying dependence structure. A real data analysis is given to explain our analytical findings.     

Semiparametric Global Cross-ratio Models for Bivariate Censored Data

Abstract.  Multivariate correlated failure time data arise in many medical and scientific settings. In the analysis of such data, it is important to use models where the parameters have simple interpretations. In this paper, we formulate a model for bivariate survival data based on the Plackett distribution. The model is an alternative to the Gamma frailty model proposed by Clayton and Oakes. The parameter in this distribution has a very appealing odds ratio interpretation for dependence between the two failure times; in addition, it allows for negative dependence. We develop novel semiparametric estimation and inference procedures for the model. The asymptotic results of the estimator are developed. The performance of the proposed techniques in finite samples is examined using simulation studies; in addition, the proposed methods are applied to data from an observational study in cancer.     

Accelerated life regression modelling of dependent bivariate time‐to‐event data

To analyze bivariate time‐to‐event data from matched or naturally paired study designs, researchers frequently use a random effect called frailty to model the dependence between within‐pair response measurements. The authors propose a computational framework for fitting dependent bivariate time‐to‐event data that combines frailty distributions and accelerated life regression models. In this framework users can choose from several parametric options for frailties, as well as the conditional distributions for within‐pair responses. The authors illustrate the flexibility that their framework represents using paired data from a study of laser photocoagulation therapy for retinopathy in diabetic patients.     

Power-normal distribution

Recently, Gupta and Gupta [Analyzing skewed data by power-normal model, Test 17 (2008), pp. 197–210] proposed the power-normal distribution for which normal distribution is a special case. The power-normal distribution is a skewed distribution, whose support is the whole real line. Our main aim of this paper is to consider bivariate power-normal distribution, whose marginals are power-normal distributions. We obtain the proposed bivariate power-normal distribution from Clayton copula, and by making a suitable transformation in both the marginals. Lindley–Singpurwalla distribution also can be used to obtain the same distribution. Different properties of this new distribution have been investigated in detail. Two different estimators are proposed. One data analysis has been performed for illustrative purposes. Finally, we propose some generalizations to multivariate case also along the same line and discuss some of its properties.     

A new test for single against competing risks models in duration analysis

This paper shows that the single-risk duration model with two event types is a limiting case of bivariate dependent competing risks model, where the joint distribution of event times are degenerate. Then a new test is proposed for the null hypothesis of single risk against dependent competing risks model under the proportional hazard model assumption.     

Modified Sarhan–Balakrishnan singular bivariate distribution

Recently Sarhan and Balakrishnan [2007. A new class of bivariate distribution and its mixture. Journal of Multivariate Analysis 98, 1508–1527] introduced a new bivariate distribution using generalized exponential and exponential distributions. They discussed several interesting properties of this new distribution. Unfortunately, they did not discuss any estimation procedure of the unknown parameters. In this paper using the similar idea as of Sarhan and Balakrishnan [2007. A new class of bivariate distribution and its mixture. Journal of Multivariate Analysis 98, 1508–1527], we have proposed a singular bivariate distribution, which has an extra shape parameter. It is observed that the marginal distributions of the proposed bivariate distribution are more flexible than the corresponding marginal distributions of the Marshall–Olkin bivariate exponential distribution, Sarhan–Balakrishnan's bivariate distribution or the bivariate generalized exponential distribution. Different properties of this new distribution have been discussed. We provide the maximum likelihood estimators of the unknown parameters using EM algorithm. We reported some simulation results and performed two data analysis for illustrative purposes. Finally we propose some generalizations of this bivariate model.     

Claim Dependence Induced by Common Effects in Hierarchical Credibility Models

In the usual credibility model, observations are made of a risk or group of risks selected from a population, and claims are assumed to be independent among different risks. However, there are some problems in practical applications and this assumption may be violated in some situations. Some credibility models allow for one source of claim dependence only, that is, across time for an individual insured risk or a group of homogeneous insured risks. Some other credibility models have been developed on a two-level common effects model that allows for two possible sources of dependence, namely, across time for the same individual risk and between risks. In this paper, we argue for the notion of modeling claim dependence on a three-level common effects model that allows for three possible sources of dependence, namely, across portfolios, across individuals and simultaneously across time within individuals. We also obtain the corresponding credibility premiums hierarchically using the projection method. Then we derive the general hierarchical structure or multi-level credibility premiums for the models with h-level of common effects.     

Gamma frailty models for bivariate survival data

In this paper, we introduce the shared gamma frailty models with two different baseline distributions namely, the generalized log-logistic and the generalized Weibull. We introduce the Bayesian estimation procedure to estimate the parameters involved in these models. We present a simulation study to compare the true values of the parameters with the estimated values. We apply these models to a real-life bivariate survival data set of McGilchrist and Aisbett related to the kidney infection data and a better model is suggested for the data.