Modeling bivariate survival data using shared inverse Gaussian frailty model

ABSTRACT

The shared frailty models are often used to model heterogeneity in survival analysis. The most common shared frailty model is a model in which hazard function is a product of a random factor (frailty) and the baseline hazard function which is common to all individuals. There are certain assumptions about the baseline distribution and the distribution of frailty. In this paper, we consider inverse Gaussian distribution as frailty distribution and three different baseline distributions, namely the generalized Rayleigh, the weighted exponential, and the extended Weibull distributions. With these three baseline distributions, we propose three different inverse Gaussian shared frailty models. We also compare these models with the models where the above-mentioned distributions are considered without frailty. We develop 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. A search of the literature suggests that currently no work has been done for these three baseline distributions with a shared inverse Gaussian frailty so far. We also apply these three models by using a real-life bivariate survival data set 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 using the Bayesian model selection criteria.     

Modeling of Inverse Gaussian Frailty Model for Bivariate Survival Data

Shared frailty models are often 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 article, we consider inverse Gaussian distribution as frailty distribution and three different baseline distributions namely, Weibull, generalized exponential, and exponential power distribution. With these three baseline distributions, we propose three different inverse Gaussian shared frailty models. To estimate the parameters involved in these models we adopt Markov Chain Monte Carlo (MCMC) approach. We present a simulation study to compare the true values of the parameters with the estimated values. Also, we apply these three models to a real life bivariate survival data set 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 kidney infection and a better model is suggested for the data.     

Bayesian Estimation of Parameters and Comparison of Shared Gamma Frailty Models

In this article, we introduce shared gamma frailty models with three different baseline distributions namely, Weibull, generalized exponential and exponential power distributions. We develop 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 three models to a real life bivariate survival dataset of McGilchrist and Aisbett (1991 McGilchrist, C. A. and Aisbett, C. W. 1991. Regression with frailty in survival analysis. Biometrics, 47: 461–466. [Crossref], [PubMed], [Web of Science ®] , [Google Scholar]) related to kidney infection data and a better model is suggested for the data.     

A Modification of a Percentile Estimation Procedure Based on Generalized Polya Urns

Adaptive designs find an important application in the estimation of unknown percentiles for an underlying dose-response curve. A nonparametric adaptive design was suggested by Mugno et al. (2004 Mugno, R.A., Zhus, W., Rosenberger, W.F. (2004). Adaptive urn designs for estimating several percentiles of a dose-response curve. Statist. Med. 23(13):2137–2150.[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]) to simultaneously estimate multiple percentiles of an unknown dose-response curve via generalized Polya urns. In this article, we examine the properties of the design proposed by Mugno et al. (2004 Mugno, R.A., Zhus, W., Rosenberger, W.F. (2004). Adaptive urn designs for estimating several percentiles of a dose-response curve. Statist. Med. 23(13):2137–2150.[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]) when delays in observing responses are encountered. Using simulations, we evaluate a modification of the design under varying group sizes. Our results demonstrate unbiased estimation with minimal loss in efficiency when compared to the original compound urn design.     

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.     

Tuning Parameter Selection in Penalized Frailty Models

The penalized likelihood approach of Fan and Li (2001 Fan, J., Li, R. (2001). Variable selection via nonconcave penalized likelihood and its oracle properties. Journal of the American Association 96:1348–1360.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar], 2002 Fan, J., Li, R. (2002). Variable selection for Cox’s proportional hazards model and frailty model. The Annals of Statistics 30:74–99.[Crossref], [Web of Science ®] , [Google Scholar]) differs from the traditional variable selection procedures in that it deletes the non-significant variables by estimating their coefficients as zero. Nevertheless, the desirable performance of this shrinkage methodology relies heavily on an appropriate selection of the tuning parameter which is involved in the penalty functions. In this work, new estimates of the norm of the error are firstly proposed through the use of Kantorovich inequalities and, subsequently, applied to the frailty models framework. These estimates are used in order to derive a tuning parameter selection procedure for penalized frailty models and clustered data. In contrast with the standard methods, the proposed approach does not depend on resampling and therefore results in a considerable gain in computational time. Moreover, it produces improved results. Simulation studies are presented to support theoretical findings and two real medical data sets are analyzed.     

Discrete Beta-Exponential Distribution

Consider the estimation of the regression parameters in the usual linear model. For design densities with infinite support, it has been shown by Faraldo Roca and González Manteiga [1] Faraldo Roca, P. and González Manteiga, W. 1987. “Efficiency of a new class of linear regression estimates obtained by preliminary nonparametric estimation”. In New Perspectives in Theoretical and Applied Statistics Edited by: Puri, M. L., Vilaplana, J. P. and Wertz, W. 229–242. New York: John Wiley.  [Google Scholar] that it is possible to modify the classical least squares procedure and to obtain estimators for the regression parameters whose MSE's (mean squared errors) are smaller than those of the usual least squares estimators. The modification consists of presmoothing the response variables by a kernel estimator of the regression function. These authors also show that the gain in efficiency is not possible for a design density with compact support. We show that in this case local linear presmoothing does not fix this inefficiency problem, in spite of the well known fact that local linear fitting automatically corrects the bias in the endpoints of the (design density) support. We demonstrate on a theoretical basis how this inefficiency problem can be rectified in the compact design case: we prove that presmoothing with boundary kernels studied in Müller [2] Müller, H.-G. 1991. Smooth optimum kernel estimators near endpoints. Biometrika, 78: 521–530. [Crossref], [Web of Science ®] , [Google Scholar] and Müller and Wang [3] Müller, H.-G. and Wang, J.-L. 1994. Hazard rate estimation under random censoring with varying kernels and bandwidths. Biometrics, 50: 61–76. [Crossref], [PubMed], [Web of Science ®] , [Google Scholar] leads to regression estimators which are superior over the least squares estimators. A very careful analytic treatment is needed to arrive at these asymptotic results.     

On the Properties of the SAVE Directions

In quadratic discriminant analysis, the use of SAVE (Cook and Weisberg, 1991 Cook, R.D., Weisberg, S. (1991). Discussion of Li (1991). J. Amer. Statist. Assoc. 86:328–32.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]; Pardoe et al., 2007 Pardoe, I., Yin, X., Cook, R. (2007). Graphical tools for quadratic discriminant analysis. Technometrics 49:172–183.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]) is often recommended for dimension-reduction purposes. However, the associated directions tend to over-emphasize the differences of the groups in dispersion, ignoring at the same time those in location. This behavior makes often the plots of the corresponding canonical coordinates difficult to interpret. In this article, the properties of SAVE are investigated and related to those of the SIR and SIRII components. Applications with real data are presented. Comparisons with previous work in this area are also discussed.     

Finite-sample distributions of the Wald,likelihood ratio,and Lagrange multiplier test statistics in the classical linear model

In this note, it is shown that the finite-sample distributions of the Wald, likelihood ratio, and Lagrange multiplier statistics in the classical linear regression model are members of the generalized beta model introduced by McDonald and Xu (1995a McDonald, J.B., Xu, Y.J. (1995a). A generalization of the beta distribution with applications. J. Econom. 66:133–152.[Crossref], [Web of Science ®] , [Google Scholar]). This is useful for examining the properties of these test statistics. For example, this characterization makes it easy to find distribution, quantile, and density functions for each test statistic, makes it clear why Wald tests may overreject the null hypothesis using asymptotic critical values, and formalizes the fact that the Lagrange multiplier statistic follows a distribution with bounded support.     

Linear transformation models for survival analysis with tumor growth information in cancer screening study

The complication in analyzing tumor data is that the tumors detected in a screening program tend to be slowly progressive tumors, which is the so-called left-truncated sampling that is inherent in screening studies. Under the assumption that all subjects have the same tumor growth function, Ghosh (2008 Ghosh, D. (2008). Proportional hazards regression for cancer studies. Biometrics 64:141–148.[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]) developed estimation procedures for the Cox proportional hazards model. Shen (2011a Shen, P.-S. (2011a). Proportional hazards regression for cancer screening data. J. Stat. Comput. Simul. 18:367–377.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]) demonstrated that Ghosh (2008 Ghosh, D. (2008). Proportional hazards regression for cancer studies. Biometrics 64:141–148.[Crossref], [PubMed], [Web of Science ®] , [Google Scholar])'s approach can be extended to the case when each subject has a specific growth function. In this article, under linear transformation model, we present a general framework to the analysis of data from cancer screening studies. We developed estimation procedures under linear transformation model, which includes Cox's model as a special case. A simulation study is conducted to demonstrate the potential usefulness of the proposed estimators.     

On extreme order statistics from heterogeneous beta distributions with applications

This paper treats the problem of stochastic comparisons for the extreme order statistics arising from heterogeneous beta distributions. Some sufficient conditions involved in majorization-type partial orders are provided for comparing the extreme order statistics in the sense of various magnitude orderings including the likelihood ratio order, the reversed hazard rate order, the usual stochastic order, and the usual multivariate stochastic order. The results established here strengthen and extend those including Kochar and Xu (2007 Kochar, S.C., Xu, M. (2007). Stochastic comparisons of parallel systems when components have proportional hazard rates. Probab. Eng. Inf. Sci. 21:597–609.[Crossref], [Web of Science ®] , [Google Scholar]), Mao and Hu (2010 Mao, T., Hu, T. (2010). Equivalent characterizations on orderings of order statistics and sample ranges. Probab. Eng. Inf. Sci. 24:245–262.[Crossref], [Web of Science ®] , [Google Scholar]), Balakrishnan et al. (2014 Balakrishnan, N., Barmalzan, G., Haidari, A. (2014). On usual multivariate stochastic ordering of order statistics from heterogeneous beta variables. J. Multivariate Anal. 127:147–150.[Crossref], [Web of Science ®] , [Google Scholar]), and Torrado (2015 Torrado, N. (2015). On magnitude orderings between smallest order statistics from heterogeneous beta distributions. J. Math. Anal. Appl. 426:824–838.[Crossref], [Web of Science ®] , [Google Scholar]). A real application in system assembly and some numerical examples are also presented to illustrate the theoretical results.     

Universally optimal designs under interference models with and without block effects

The concept of neighbor designs was introduced and defined by Rees (1967 Rees, D.H. (1967). Some designs of use in serology. Biometrics 23:779–791.[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]) along with giving some methods of their construction. Henceforth, many methods of construction of neighbor designs as well as of their generalizations are available in the literature. However, there are only few results on their optimality. Therefore, the purpose of this article is to give an overview of study on this problem. Recent results on optimality of specified neighbor balanced designs under various interference models with block effects are presented and then these results are compared with respective models where block effects are not significant.     

Generalized linear regression models incorporating original outcome distributions

ABSTRACT

In this article, we propose an approach for incorporating continuous and discrete original outcome distributions into the usual exponential family regression models. The new approach is an extension of the works of Suissa (1991 Suissa, S. (1991). Binary methods for continuous outcomes: A parametric alternative. J. Clin. Epidemiol. 44:241–248.[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]) and Suissa and Blais (1995 Suissa, S., Blais, L. (1995). Binary regression with continuous outcomes. Stat. Med. 14:247–255.[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]), which present methods to estimate the risk of an event defined in a sample subspace of an original continuous outcome variable. Simulation studies are presented in order to illustrate the performance of the developed methodology. Real data sets are analyzed by using the proposed models.     

Bimodal symmetric-asymmetric power-normal families

This article proposes new symmetric and asymmetric distributions applying methods analogous as the ones in Kim (2005 Kim, H.J. (2005). On a class of two-piece skew-normal distributions. Statist.: J. Theoret. Appl. Statist. 39:537–553.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]) and Arnold et al. (2009 Arnold, B.C., H.W. Gómez, and H.S. Salinas. (2009). On multiple constraint skewed models. Statist. J. Theoret. Appl. Statist. 43: 279–293.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]) to the exponentiated normal distribution studied in Durrans (1992 Durrans, S.R. (1992). Distributions of fractional order statistics in hydrology. Water Resour. Res. 28:1649–1655.[Crossref], [Web of Science ®] , [Google Scholar]), that we call the power-normal (PN) distribution. The proposed bimodal extension, the main focus of the paper, is called the bimodal power-normal model and is denoted by BPN(α) model, where α is the asymmetry parameter. The authors give some properties including moments and maximum likelihood estimation. Two important features of the model proposed is that its normalizing constant has closed and simple form and that the Fisher information matrix is nonsingular, guaranteeing large sample properties of the maximum likelihood estimators. Finally, simulation studies and real applications reveal that the proposed model can perform well in both situations.     

An efficient use of moment's ratios of scrambling variables in a randomized response technique

This paper aimed at providing an efficient new unbiased estimator for estimating the proportion of a potentially sensitive attribute in survey sampling. The suggested randomization device makes use of the means, variances of scrambling variables, and the two scalars lie between “zero” and “one.” Thus, the same amount of information has been used at the estimation stage. The variance formula of the suggested estimator has been obtained. We have compared the proposed unbiased estimator with that of Kuk (1990 Kuk, A.Y.C. (1990). Asking sensitive questions inderectely. Biometrika 77:436–438.[Crossref], [Web of Science ®] , [Google Scholar]) and Franklin (1989 Franklin, L.A. (1989). A comparision of estimators for randomized response sampling with continuous distribution s from a dichotomous population. Commun. Stat. Theor. Methods 18:489–505.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]), and Singh and Chen (2009 Singh, S., Chen, C.C. (2009). Utilization of higher order moments of scrambling variables in randomized response sampling. J. Stat. Plann. Inference. 139:3377–3380.[Crossref], [Web of Science ®] , [Google Scholar]) estimators. Relevant conditions are obtained in which the proposed estimator is more efficient than Kuk (1990 Kuk, A.Y.C. (1990). Asking sensitive questions inderectely. Biometrika 77:436–438.[Crossref], [Web of Science ®] , [Google Scholar]) and Franklin (1989 Franklin, L.A. (1989). A comparision of estimators for randomized response sampling with continuous distribution s from a dichotomous population. Commun. Stat. Theor. Methods 18:489–505.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]) and Singh and Chen (2009 Singh, S., Chen, C.C. (2009). Utilization of higher order moments of scrambling variables in randomized response sampling. J. Stat. Plann. Inference. 139:3377–3380.[Crossref], [Web of Science ®] , [Google Scholar]) estimators. The optimum estimator (OE) in the proposed class of estimators has been identified which finally depends on moments ratios of the scrambling variables. The variance of the optimum estimator has been obtained and compared with that of the Kuk (1990 Kuk, A.Y.C. (1990). Asking sensitive questions inderectely. Biometrika 77:436–438.[Crossref], [Web of Science ®] , [Google Scholar]) and Franklin (1989 Franklin, L.A. (1989). A comparision of estimators for randomized response sampling with continuous distribution s from a dichotomous population. Commun. Stat. Theor. Methods 18:489–505.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]) estimator and Singh and Chen (2009 Singh, S., Chen, C.C. (2009). Utilization of higher order moments of scrambling variables in randomized response sampling. J. Stat. Plann. Inference. 139:3377–3380.[Crossref], [Web of Science ®] , [Google Scholar]) estimator. It is interesting to mention that the “optimum estimator” of the class of estimators due to Singh and Chen (2009 Singh, S., Chen, C.C. (2009). Utilization of higher order moments of scrambling variables in randomized response sampling. J. Stat. Plann. Inference. 139:3377–3380.[Crossref], [Web of Science ®] , [Google Scholar]) depends on the parameter π under investigation which limits the use of Singh and Chen (2009 Singh, S., Chen, C.C. (2009). Utilization of higher order moments of scrambling variables in randomized response sampling. J. Stat. Plann. Inference. 139:3377–3380.[Crossref], [Web of Science ®] , [Google Scholar]) OE in practice while the proposed OE in this paper is free from such a constraint. The proposed OE depends only on the moments ratios of scrambling variables. This is an advantage over the Singh and Chen (2009 Singh, S., Chen, C.C. (2009). Utilization of higher order moments of scrambling variables in randomized response sampling. J. Stat. Plann. Inference. 139:3377–3380.[Crossref], [Web of Science ®] , [Google Scholar]) estimator. Numerical illustrations are given in the support of the present study when the scrambling variables follow normal distribution. Theoretical and empirical results are very sound and quite illuminating in the favor of the present study.     

The Marshall–Olkin extended generalized Rayleigh distribution: Properties and applications

In this paper, a new extension for the generalized Rayleigh distribution is introduced. The proposed model, called Marshall–Olkin extended generalized Rayleigh distribution, arises based on the scheme introduced by Marshall and Olkin (1997) Marshall, A.W., Olkin, I. (1997). A new method for adding a parameter to a family of distributions with application to the exponential and Weibull families. Biometrika 84:641–652.[Crossref], [Web of Science ®] , [Google Scholar]. A comprehensive account of the mathematical properties of the new distribution is provided. We discuss about the estimation of the model parameters based on two estimation methods. Empirical applications of the new model to real data are presented for illustrative purposes.     

Inequalities for sums of squares of reranked differences involving ties and missing data

An inequality for the sum of squares of rank differences associated with Spearman’s rank correlation coefficient, when ties and missing data are present in both rankings, was established numerically in Loukas and Papaioannou (1991 Loukas, S., Papaioannou, T. (1991). Rank correlation inequalities with ties and missing data. Stat. Probab. Lett. 11:53–56.[Crossref], [Web of Science ®] , [Google Scholar]). That inequality is improved and generalized.     

Using the Inverse Maximum Ratio-Λ as a Technique to Quantify Surface Uniformity

Techniques used in variability assessment are subsequently used to draw conclusions regarding the “spread”/uniformity of data curves. Due to the limitations of these techniques, they are not adequate for circumstances where data manifest with multiple peaks. Examples of these manifestations (in three-dimensional space) include under-foot pressure distributions recorded for different types of footwear (Becerro-de-Bengoa-Vallejo et al., 2014 Biau, D.J. (2011). In brief: Standard deviation and standard error. Clinical Orthopaedics and Related Research 469(9):2661–2664.[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]; Cibulka et al., 1994 Cibulka, M.T., Sinacore, D.R., Mueller, M.J. (1994). Shin splints and forefoot contact running: A case report. Journal of Orthopaedic &; Sports Physical Therapy 20(2):98–102.[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]; Davies et al., 2003 Davies, M.B., Betts, R.P., Scott, I.R. (2003). Optical plantar pressure analysis following internal fixation for displaced intra-articular os calcis fractures. Foot &; Ankle International 24(11):851–856.[PubMed], [Web of Science ®] , [Google Scholar]), surface textures and interfaces designed to impact friction, and and and molecular surface structures such as viral epitopes (Torras and Garcia-Valls, 2004 Torras, C., Garcia-Valls, R. (2004). Quantification of membrane morphology by interpretation of scanning electron microscopy images. Journal of Membrane Science 233(1–2):119–127.[Crossref], [Web of Science ®] , [Google Scholar]; Pacejka, 1997; Fustaffson, 1997). This article proposes a technique for generating a single variable – Λ that will quantify the uniformity of such surfaces. We define and validate this technique using several mathematical and graphical models.     

Bootstrap methods for multivariate hypothesis testing

The nonparametric and parametric bootstrap methods for multivariate hypothesis testing are developed. They are used to approximate the null distribution of the test statistics proposed by Duchesne and Francq (2015 Duchesne, P., Francq, C. (2015). Multivariate hypothesis testing using generalized and {2}-inverses—with applications. Statistics 49:475–496.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]), resulting in bootstrap testing procedures. In the problem of testing for the mean vector of a multivariate distribution, the asymptotic validity of the bootstrap methods is proved. The finite sample performance of the new solutions is demonstrated by means of Monte Carlo simulation studies. They indicate that for small-sample size, the bootstrap tests provide a better finite sample properties than the asymptotic tests considered by Duchesne and Francq (2015 Duchesne, P., Francq, C. (2015). Multivariate hypothesis testing using generalized and {2}-inverses—with applications. Statistics 49:475–496.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]).     

The odd log-logistic generalized half-normal lifetime distribution: Properties and applications

Cooray and Ananda (2008 Cooray, K., Ananda, M.M.A. (2008). A Generalization of the half-normal distribution with applications to lifetime data. Commun. Stat. - Theory Methods 37:1323–1337.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]) pioneered a lifetime model commonly used in reliability studies. Based on this distribution, we propose a new model called the odd log-logistic generalized half-normal distribution for describing fatigue lifetime data. Various of its structural properties are derived. We discuss the method of maximum likelihood to fit the model parameters. For different parameter settings and sample sizes, some simulation studies compare the performance of the new lifetime model. It can be very useful, and its superiority is illustrated by means of a real dataset.