Correlated gamma frailty models for bivariate survival data

Frailty models are used in the survival analysis to account for the unobserved heterogeneity in individual risks to disease and death. To analyze the bivariate data on related survival times (e.g., matched pairs experiments, twin or family data) the shared frailty models were suggested. Shared frailty models are used despite their limitations. To overcome their disadvantages correlated frailty models may be used. In this article, we introduce the gamma correlated frailty models with two different baseline distributions namely, the generalized log logistic, and the generalized Weibull. We introduce the Bayesian estimation procedure using Markov chain Monte Carlo (MCMC) technique 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. Also we apply these models to a real life bivariate survival dataset related to the kidney infection data and a better model is suggested for the data.     

Shared frailty models with baseline generalized Pareto distribution

Abstract

In this article, we have considered three different shared frailty models under the assumption of generalized Pareto Distribution as baseline distribution. Frailty models have been used in the survival analysis to account for the unobserved heterogeneity in an individual risks to disease and death. These three frailty models are with gamma frailty, inverse Gaussian frailty and positive stable frailty. Then we introduce the Bayesian estimation procedure using Markov chain Monte Carlo (MCMC) technique to estimate the parameters. We applied these three models to a kidney infection data and find the best fitted model for kidney infection data. We present a simulation study to compare true value of the parameters with the estimated values. Model comparison is made using Bayesian model selection criterion and a well-fitted model is suggested for the kidney infection data.     

Shared gamma frailty models based on reversed hazard rate for modeling Australian twin data

In this paper, we consider shared gamma frailty model with the reversed hazard rate (RHR) with two different baseline distributions, namely the generalized inverse Rayleigh and the exponentiated Gumbel distributions. With these two baseline distributions we propose two different shared frailty models. We develop the Bayesian estimation procedure using Markov Chain Monte Carlo technique 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. A search of the literature suggests that currently no work has been done for these two baseline distributions with a shared gamma frailty with the RHR so far. We also apply these two models by using a real life bivariate survival data set of Australian twin data given by Duffy et a1. (1990) and a better model is suggested for the data.     

Shared frailty models based on reversed hazard rate for modified inverse Weibull distribution as baseline distribution

The unknown or unobservable risk factors in the survival analysis cause heterogeneity between individuals. Frailty models are used in the survival analysis to account for the unobserved heterogeneity in individual risks to disease and death. To analyze the bivariate data on related survival times, the shared frailty models were suggested. The most common shared frailty model is a model in which frailty act multiplicatively on the hazard function. In this paper, we introduce the shared gamma frailty model and the inverse Gaussian frailty model with the reversed hazard rate. We introduce the Bayesian estimation procedure using Markov chain Monte Carlo (MCMC) technique to estimate the parameters involved in the model. We present a simulation study to compare the true values of the parameters with the estimated values. We also apply the proposed models to the Australian twin data set and a better model is suggested.     

Shared frailty model based on reversed hazard rate for left censored data

Frailty models are used in the survival analysis to account for the unobserved heterogeneity in the individual risks to disease and death. To analyze the bivariate data on related survival times (e.g., matched pairs experiments, twin or family data), the shared frailty models were suggested. In this article, we introduce the shared gamma frailty models with the reversed hazard rate. We develop the Bayesian estimation procedure using the Markov chain Monte Carlo (MCMC) technique to estimate the parameters involved in the model. We present a simulation study to compare the true values of the parameters with the estimated values. We apply the model to a real life bivariate survival dataset.     

Gamma shared frailty model based on reversed hazard rate

Abstract

Frailty models are used in survival analysis to account for unobserved heterogeneity in individual risks to disease and death. To analyze bivariate data on related survival times (e.g., matched pairs experiments, twin, or family data), shared frailty models were suggested. Shared frailty models are frequently used to model heterogeneity in survival analysis. The most common shared frailty model is a model in which hazard function is a product of random factor(frailty) and baseline hazard function which is common to all individuals. There are certain assumptions about the baseline distribution and distribution of frailty. In this paper, we introduce shared gamma frailty models with reversed hazard rate. We introduce Bayesian estimation procedure using Markov Chain Monte Carlo (MCMC) technique to estimate the parameters involved in the model. We present a simulation study to compare the true values of the parameters with the estimated values. Also, we apply the proposed model to the Australian twin data set.     

The log-exponentiated generalized gamma regression model for censored data

For the first time, we introduce a generalized form of the exponentiated generalized gamma distribution [Cordeiro et al. The exponentiated generalized gamma distribution with application to lifetime data, J. Statist. Comput. Simul. 81 (2011), pp. 827–842.] that is the baseline for the log-exponentiated generalized gamma regression model. The new distribution can accommodate increasing, decreasing, bathtub- and unimodal-shaped hazard functions. A second advantage is that it includes classical distributions reported in the lifetime literature as special cases. We obtain explicit expressions for the moments of the baseline distribution of the new regression model. The proposed model can be applied to censored data since it includes as sub-models several widely known regression models. It therefore can be used more effectively in the analysis of survival data. We obtain maximum likelihood estimates for the model parameters by considering censored data. We show that our extended regression model is very useful by means of two applications to real data.     

Shared inverse Gaussian frailty models based on additive hazards

Frailty models are used in the survival analysis to account for the unobserved heterogeneity in individual risks to disease and death. To analyze the bivariate data on related survival times (e.g., matched pairs experiments, twin, or family data), the shared frailty models were suggested. These models are based on the assumption that frailty acts multiplicatively to hazard rate. In this article, we assume that frailty acts additively to hazard rate. We introduce the shared inverse Gaussian frailty models with three different baseline distributions, namely the generalized log-logistic, the generalized Weibull, and exponential power distribution. We introduce the Bayesian estimation procedure using Markov chain Monte Carlo technique to estimate the parameters involved in these models. We apply these models to a real-life bivariate survival dataset of McGilchrist and Aisbett (1991 McGilchrist, C.A., Aisbett, C.W. (1991). Regression with frailty in survival analysis. Biometrics 47:461–466.[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]) related to the kidney infection data, and a better model is suggested for the data.     

Modeling Australian twin data using shared positive stable frailty models based on reversed hazard rate

Unobserved heterogeneity, also called frailty, is a major concern in the application of survival analysis. The shared frailty models allow for the statistical dependence between the observed survival data. In this paper, we consider shared positive stable frailty model with the reversed hazard rate (RHR) with three different baseline distributions, namely the exponentiated Gumbel, the generalized Rayleigh, and the generalized inverse Rayleigh distributions. With these three baseline distributions we propose three different shared frailty models. We develop the Bayesian estimation procedure using Markov Chain Monte Carlo technique 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. A search of the literature suggests that currently no work has been done for these three baseline distributions with a shared positive stable frailty with the RHR so far. We also apply these three models by using a real-life bivariate survival data set of Australian twin data given by Duffy et a1. (1990 Duffy, D.L., Martin, N.G., Mathews, J.D. (1990). Appendectomy in Australian twins. Aust. J. Hum. Genet. 47(3):590–592.[PubMed], [Web of Science ®] , [Google Scholar]) and a better model is suggested for the data.     

Generalized Poisson–Lindley linear model for count data

The purpose of this paper is to develop a new linear regression model for count data, namely generalized-Poisson Lindley (GPL) linear model. The GPL linear model is performed by applying generalized linear model to GPL distribution. The model parameters are estimated by the maximum likelihood estimation. We utilize the GPL linear model to fit two real data sets and compare it with the Poisson, negative binomial (NB) and Poisson-weighted exponential (P-WE) models for count data. It is found that the GPL linear model can fit over-dispersed count data, and it shows the highest log-likelihood, the smallest AIC and BIC values. As a consequence, the linear regression model from the GPL distribution is a valuable alternative model to the Poisson, NB, and P-WE models.     

The gamma-exponentiated exponential distribution

In this paper, we introduce a new distribution generated by gamma random variables. We show that this distribution includes as a special case the distribution of the lower record value from a sequence of i.i.d. random variables from a population with the exponentiated (generalized) exponential distribution. The properties of this distribution are derived and the estimation of the model parameters is discussed. Some applications to real data sets are finally presented for illustration.     

A new exponentiated class of distributions: Properties and applications

Motivated by series-parallel-series systems, we introduce a new generator of continuous distributions with three extra parameters called the generalized exponentiated class of distributions. We study its mathematical properties and introduce a bivariate extension of the class. We discuss the estimation of its parameters by maximum likelihood and illustrate the potentiality of the class by applications to two real datasets.     

Bayesian Conditional Mean Estimation in Log‐Normal Linear Regression Models with Finite Quadratic Expected Loss

Log‐normal linear regression models are popular in many fields of research. Bayesian estimation of the conditional mean of the dependent variable is problematic as many choices of the prior for the variance (on the log‐scale) lead to posterior distributions with no finite moments. We propose a generalized inverse Gaussian prior for this variance and derive the conditions on the prior parameters that yield posterior distributions of the conditional mean of the dependent variable with finite moments up to a pre‐specified order. The conditions depend on one of the three parameters of the suggested prior; the other two have an influence on inferences for small and medium sample sizes. A second goal of this paper is to discuss how to choose these parameters according to different criteria including the optimization of frequentist properties of posterior means.     

Generalized Laplacian Distributions and Autoregressive Processes

Generalized Laplacian distribution is considered. A new distribution called geometric generalized Laplacian distribution is introduced and its properties are studied. First- and higher-order autoregressive processes with these stationary marginal distributions are developed and studied. Simulation studies are conducted and trajectories of the process are obtained for selected values of the parameters. Various areas of application of these models are discussed.     

GQL estimation in linear dynamic models for panel data

In the econometrics literature, it is standard practice to use the existing instrumental variables as well as generalized method of moments approaches for the estimation of the parameters of a linear dynamic mixed model for panel data. In this paper, we introduce a generalized quasi-likelihood estimation approach that produces estimates with smaller mean squared errors when compared with the aforementioned and other existing approaches.     

Modified slashed generalized exponential distribution

Abstract

In this article, we introduce a new distribution for modeling positive data sets with high kurtosis, the modified slashed generalized exponential distribution. The new model can be seen as a modified version of the slashed generalized exponential distribution. It arises as a quotient of two independent random variables, one being a generalized exponential distribution in the numerator and a power of the exponential distribution in the denominator. We studied various structural properties (such as the stochastic representation, density function, hazard rate function and moments) and discuss moment and maximum likelihood estimating approaches. Two real data sets are considered in which the utility of the new model in the analysis with high kurtosis is illustrated.     

The skew generalized t distribution as the scale mixture of a skew exponential power distribution and its applications in robust estimation

In this paper, we consider the family of skew generalized t (SGT) distributions originally introduced by Theodossiou [P. Theodossiou, Financial data and the skewed generalized t distribution, Manage. Sci. Part 1 44 (12) ( 1998), pp. 1650–1661] as a skew extension of the generalized t (GT) distribution. The SGT distribution family warrants special attention, because it encompasses distributions having both heavy tails and skewness, and many of the widely used distributions such as Student's t, normal, Hansen's skew t, exponential power, and skew exponential power (SEP) distributions are included as limiting or special cases in the SGT family. We show that the SGT distribution can be obtained as the scale mixture of the SEP and generalized gamma distributions. We investigate several properties of the SGT distribution and consider the maximum likelihood estimation of the location, scale, and skewness parameters under the assumption that the shape parameters are known. We show that if the shape parameters are estimated along with the location, scale, and skewness parameters, the influence function for the maximum likelihood estimators becomes unbounded. We obtain the necessary conditions to ensure the uniqueness of the maximum likelihood estimators for the location, scale, and skewness parameters, with known shape parameters. We provide a simple iterative re-weighting algorithm to compute the maximum likelihood estimates for the location, scale, and skewness parameters and show that this simple algorithm can be identified as an EM-type algorithm. We finally present two applications of the SGT distributions in robust estimation.     

A new class of generalized Lindley distributions with applications

A new four-parameter class of generalized Lindley (GL) distribution called the beta-generalized Lindley (BGL) distribution is proposed. This class of distributions contains the beta-Lindley, GL and Lindley distributions as special cases. Expansion of the density of the BGL distribution is obtained. The properties of these distributions, including hazard function, reverse hazard function, monotonicity property, shapes, moments, reliability, mean deviations, Bonferroni and Lorenz curves are derived. Measures of uncertainty such as Renyi entropy and s-entropy as well as Fisher information are presented. Method of maximum likelihood is used to estimate the parameters of the BGL and related distributions. Finally, real data examples are discussed to illustrate the applicability of this class of models.     

A new generalization of the Weibull-geometric distribution with bathtub failure rate

In this paper, the researchers attempt to introduce a new generalization of the Weibull-geometric distribution. The failure rate function of the new model is found to be increasing, decreasing, upside-down bathtub, and bathtub-shaped. The researchers obtained the new model by compounding Weibull distribution and discrete generalized exponential distribution of a second type, which is a generalization of the geometric distribution. The new introduced model contains some previously known lifetime distributions as well as a new one. Some basic distributional properties and moments of the new model are discussed. Estimation of the parameters is illustrated and the model with two known real data sets is examined.     

The shape of the hazard function: Does the generalized gamma have the last word?

This paper compares the five-parameter beta generalized gamma (BGG) distribution to the three-parameter generalized gamma (GG). Both distributions include the four standard hazard shapes that we believe is an important property for any parametric family. For several BGG distributions, we select matching GGs and compute the Kullback-Liebler distance, observing remarkable agreement. We explore the beta parameters' influence on the matched GG parameters, detecting a strong connection between the distributions. Lastly, we compare the distributions using two real-data examples. We conclude from these comparisons that the BGG is not likely to be more useful for analytical purposes than the simpler GG.