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.     

On Fitting Transformation Model to Survey Data

In this article, we consider fitting a semiparametric linear model to survey data with censored observations. The specific goal of the paper is to extend the methods of Cheng et al. (1995 Cheng, S.C., Wei, L.J., Ying, Z. (1995). Analysis of transformation models with censored data. Biometrika 82(4):835–845.[Crossref], [Web of Science ®] , [Google Scholar]) and Chen et al. (2002 Chen, K., Jin, Z. Ying, Z. (2002). Semiparametric analysis of transformation models with censored data. Biometrika 89:659–668.[Crossref], [Web of Science ®] , [Google Scholar]) to the case when the sample has been drawn from a population using a complex sampling design. Similar to the approach of Lin (2000 Lin, D.Y. (2000). On fitting Cox’s proportional hazards models to survey data. Biometrika 87:37–47.[Crossref], [Web of Science ®] , [Google Scholar]), we regard the survey population as a random sample from an infinite universe and accounts for this randomness in the statistical inference. A simulation study is conducted to investigate the performance of the proposed estimators.     

An objective Bayesian estimation in a two-period crossover design

Two-period crossover design is one of the commonly used designs in clinical trials. But, the estimation of treatment effect is complicated by the possible presence of carryover effect. It is known that ignoring the carryover effect when it exists can lead to poor estimates of the treatment effect. The classical approach by Grizzle (1965 Grizzle, J.E. (1965). The two-period change-over design and its use in clinical trials. Biometrics 21:467–480. See Grizzle (1974) for corrections.[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]) consists of two stages. First, a preliminary test is conducted on carryover effect. If the carryover effect is significant, analysis is based only on data from period one; otherwise, analysis is based on data from both periods. A Bayesian approach with improper priors was proposed by Grieve (1985 Grieve, A.P. (1985). A Bayesian analysis of the two-period crossover design for clinical trials. Biometrics 41:979–990.[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]) which uses a mixture of two models: a model with carryover effect and another without. The indeterminacy of the Bayes factor due to the arbitrary constant in the improper prior was addressed by assigning a minimally discriminatory value to the constant. In this article, we present an objective Bayesian estimation approach to the two-period crossover design which is also based on a mixture model, but using the commonly recommended Zellner–Siow g-prior. We provide simulation studies and a real data example and compare the numerical results with Grizzle (1965 Grizzle, J.E. (1965). The two-period change-over design and its use in clinical trials. Biometrics 21:467–480. See Grizzle (1974) for corrections.[Crossref], [PubMed], [Web of Science ®] , [Google Scholar])’s and Grieve (1985 Grieve, A.P. (1985). A Bayesian analysis of the two-period crossover design for clinical trials. Biometrics 41:979–990.[Crossref], [PubMed], [Web of Science ®] , [Google Scholar])’s approaches.     

An adaptive multivariate CUSUM control chart for signaling a range of location shifts

ABSTRACT

In this work, we proposed an adaptive multivariate cumulative sum (CUSUM) statistical process control chart for signaling a range of location shifts. This method was based on the multivariate CUSUM control chart proposed by Pignatiello and Runger (1990 Pignatiello, J.J., Runger, G.C. (1990). Comparisons of multivariate CUSUM charts. J. Qual. Technol. 22(3):173–186.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]), but we adopted the adaptive approach similar to that discussed by Dai et al. (2011 Dai, Y., Luo, Y., Li, Z., Wang, Z. (2011). A new adaptive CUSUM control chart for detecting the multivariate process mean. Qual. Reliab. Eng. Int. 27(7):877–884.[Crossref], [Web of Science ®] , [Google Scholar]), which was based on a different CUSUM method introduced by Crosier (1988 Crosier, R.B. (1988). Multivariate generalizations of cumulative sum quality-control schemes. Technometrics 30(3):291–303.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]). The reference value in this proposed procedure was changed adaptively in each run, with the current mean shift estimated by exponentially weighted moving average (EWMA) statistic. By specifying the minimal magnitude of the mean shift, our proposed control chart achieved a good overall performance for detecting a range of shifts rather than a single value. We compared our adaptive multivariate CUSUM method with that of Dai et al. (2001 Dai, Y., Luo, Y., Li, Z., Wang, Z. (2011). A new adaptive CUSUM control chart for detecting the multivariate process mean. Qual. Reliab. Eng. Int. 27(7):877–884.[Crossref], [Web of Science ®] , [Google Scholar]) and the non adaptive versions of these two methods, by evaluating both the steady state and zero state average run length (ARL) values. The detection efficiency of our method showed improvements over the comparative methods when the location shift is unknown but falls within an expected range.     

On goodness-of-fit tests for Aalen's additive model based on left-truncated right-censored or doubly censored data

ABSTRACT

Gandy and Jensen (2005 Gandy, A., Jensen, U. (2005). On goodness-of-fit tests for Aalen's additive risk model. Scan. J. Stat. 32:425–445.[Crossref], [Web of Science ®] , [Google Scholar]) proposed goodness-of-fit tests for Aalen's additive risk model. In this article, we demonstrate that the approach of Gandy and Jensen (2005 Gandy, A., Jensen, U. (2005). On goodness-of-fit tests for Aalen's additive risk model. Scan. J. Stat. 32:425–445.[Crossref], [Web of Science ®] , [Google Scholar]) can be applied to left-truncated right-censored (LTRC) data and doubly censored data. A simulation study is conducted to investigate the performance of the proposed tests. The proposed tests are illustrated using heart transplant data.     

Alternative Approximations to Value-At-Risk: A Comparison

This article compares three value-at-risk (VaR) approximation methods suggested in the literature: Cornish and Fisher (1937 Cornish, E.A., Fisher, R.A. (1937). Moments and cumulants in the specification of distributions. Revue de l’Institut International de Statistique 5:307–320.[Crossref] , [Google Scholar]), Sillitto (1969 Sillitto, G.P. (1969). Derivation of approximants to the inverse distribution function of a continuous univariate population from the order statistics of a sample. Biometrika 56:641–650.[Crossref], [Web of Science ®] , [Google Scholar]), and Liu (2010 Liu, W.-H. (2010). Estimation and testing of portfolio value-at-risk based on L-comoment matrices. Journal of Futures Markets 30:897–908.[Crossref], [Web of Science ®] , [Google Scholar]). Simulation results are obtained for three families of distributions: student-t, skewed-normal, and skewed-t. We recommend the Sillitto approximation as the best method to evaluate the VaR when the financial return has an unknown, skewed, and heavy-tailed distribution.     

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.     

A Strong Limit Theorem for Weighted Sums of Negatively Dependent Random Variables

In this article, we establish a new complete convergence theorem for weighted sums of negatively dependent random variables. As corollaries, many results on the almost sure convergence and complete convergence for weighted sums of negatively dependent random variables are obtained. In particular, the results of Jing and Liang (2008 Jing, B.Y., Liang, H.Y. (2008). Strong limit theorems for weighted sums of negatively associated random variables. J. Theor. Probab. 21:890–909.[Crossref], [Web of Science ®] , [Google Scholar]), Sung (2012 Sung, S.H. (2012). Complete convergence for weighted sums of negatively dependent random variables. Stat. Pap. 53:73–82.[Crossref], [Web of Science ®] , [Google Scholar]), and Wu (2010) can be obtained.     

Quantile-based Chernoff distance for truncated random variables

Several probability distributions such as power-Pareto distribution (see Gilchrist 2000 Gilchrist, W. 2000. Statistical modelling with quantile functions. Boca Raton, FL: Chapman and Hall/CRC.[Crossref] , [Google Scholar] and Hankin and Lee 2006 Hankin, R. K. S., and A. Lee. 2006. A new family of non-negative distributions. Australian and New Zealand Journal of Statistics 48:67–78.[Crossref], [Web of Science ®] , [Google Scholar]), various forms of lambda distributions (see Ramberg and Schmeiser 1974 Ramberg, J. S., and B. W. Schmeiser. 1974. An appropriate method for generating asymmetric random variables. Communications of the ACM 17:78–82.[Crossref], [Web of Science ®] , [Google Scholar] and Freimer et al. 1988 Freimer, M., S. Mudholkar, G. Kollia, and C. T. Lin. 1988. A study of the generalized lambda family. Communications in Statistics - Theory and Methods 17:3547–67.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]), Govindarajulu distribution (see Nair, Sankaran, and Vineshkumar 2012 Nair, U. N., P. G. Sankaran, and B. Vineshkumar. 2012. The Govindarajulu distribution: some properties and applications. Communications in Statistics—Theory and Methods 41:4391–406.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]), etc., do not have manageable distribution functions, though they have tractable quantile functions. Hence, analytical study of the properties of Chernoff distance of two random variables associated with these distributions via traditional distribution function-based tool becomes difficult. To make this simple, in this paper, we introduce quantile-based Chernoff distance for (left or right) truncated random variables and study its various properties. Some useful bounds as well as characterization results are obtained.     

Further Results for a General Family of Bivariate Copulas

A bivariate family of copulas has been initiated by Cuadras-Augé (1981 Cuadras, C.M., Augé, J. (1981). A continuous general multivariate distribution and its properties. Commun. Statist. (A) Theor. Meth. 10:339–353.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) and Marshall (1996 Marshall, A.W. (1996). Copulas, marginals, and joint distributions. In: Distributions with fixed marginals and related topics. IMS Lecture Notes Monogr. Ser. 28:213–222.[Crossref] , [Google Scholar]). Recently, Durante (2007 Durante, F. (2007). A new family of symmetric bivariate copulas. C. R. Math. Acad. Sci. Paris 344:195–198.[Crossref], [Web of Science ®] , [Google Scholar]) considered this family as a general family of symmetric bivariate copulas indexed by a generator function and studied some of its dependence properties. In this article, we obtain and describe further aspects of dependence for this family. For example, we have proved that the family has positive likelihood ratio dependence structure if and only if the family reduces to some well-known copulas. We also derive several proper forms for the generator function of this family. Considering a multivariate extension of the bivariate family of copulas provided by Durante et al. (2007 Durante, F., Quesada-Molina, J.J., Flores, M. (2007). On a family of multivariate copulas for aggregation processes. Inform. Sci. 177(24):5715–5724.[Crossref], [Web of Science ®] , [Google Scholar]), some dependence properties are studied. Finally, some positive dependence stochastic orderings for two random vectors having a copula from the proposed families, are discussed.     

The Identification of Potential Longitudinal Biomarkers and Measurements of Effectiveness for Biomarkers as Surrogates in Multivariate Survival Data

A proposed method based on frailty models is used to identify longitudinal biomarkers or surrogates for a multivariate survival. This method is an extention of earlier models by Wulfsohn and Tsiatis (1997 Wulfsohn , M. S. , Tsiatis , A. A. ( 1997 ). A joint model for survival and longitudinal data measured with error . Biometrics 53 ( 1 ): 330 – 339 .[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]) and Song et al. (2002 Song , X. , Davidian , M. , Tsiatis , A. A. ( 2002 ). A Semiparametric likelihood approach to joint modeling of longitudinal and time-to-event data . Biometrics 58 ( 4 ): 742 – 753 .[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]). In this article, similar to Henderson et al. (2002 Henderson , R. , Diggle , P. J. , Dobson , A. ( 2002 ). Identification and efficacy of longitudinal markers for survival . Biostatistics 3 ( 1 ): 33 – 50 .[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]), a joint likelihood function combines the likelihood functions of the longitudinal biomarkers and the multivariate survival times. We use simulations to explore how the number of individuals, the number of time points per individual and the functional form of the random effects from the longitudianl biomarkers influence the power to detect the association of a longitudinal biomarker and the multivariate survival time. The proposed method is illustrate by using the gastric cancer data.     

A Simple Corrected Score for Logistic Regression with Errors-in-Covariates

We develop a simple corrected score for logistic regression with errors-in-covariates. The new method is an extension of the consistent functional methods proposed by Huang and Wang (2001) and is closely related to the corrected score method by Nakamura (1990 Nakamura, T. (1990). Corrected score function for errors-in-variables models: Methodology and application to generalized linear models. Biometrika. 77:127–137.[Crossref], [Web of Science ®] , [Google Scholar]) and Stefanski (1989 Stefanski, L.A. (1989). Unbiased estimation of a nonlinear function a normal mean with application to measurement error models. Commun. Stat. Ser. A - Theory Methods. 18:4335–4358.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]). The new method requires that the measurement error distribution is known, but does not require normality. The new method yields a consistent and asymptotically normal estimator under regularity conditions. We examine the finite-sample performance of the new estimator through simulation studies. Finally, we illustrate the new method by applying it to an AIDS study.     

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.     

Bayesian Capture-Recapture Analysis: An Application in Modeling Semelparity of a Neotropical Didelphid Marsupial

We develop a series of Bayesian statistical models for estimating survival of a neotropic didelphid marsupial, the Brazilian gracile mouse opossum (Gracilinanus microtarsus). These models are based on the Cormack–Jolly–Seber model (Cormack, 1964 Cormack , R. M. ( 1964 ). Estimates of survival from the sighting of marked animals . Biometrika 51 : 429 – 438 .[Crossref], [Web of Science ®] , [Google Scholar]; Jolly 1965 Jolly , G. M. ( 1965 ). Explicit estimates from capture-recapture data with both death and immigration stochastic model . Biometrika 52 : 225 – 247 .[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]; Seber 1965 Seber , G. A. F. ( 1965 ). A note on the multiple recapture census . Biometrika 52 : 249 – 259 .[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]) with both survival and recapture rates expressed as a function of covariates using a logit link. The proposed models allow taking into account heterogeneity in capture probability caused by the existence of different groups of individuals in the population. The models were applied to two cohorts (Cohort, 2000, 2001) with the first one including 14 and the second one 15 sampling occasions. The best models for each of the cohorts indicate that G. microtarsus is best described as partially semelparous, a condition in which mortality after the first mating is high but graded over time, with a fraction of males surviving for a second breeding season (Boonstra, 2005 Boonstra , R. ( 2005 ). Equipped for life: the adaptive role of the stress axis in male mammals . Journal of Mammalogy 86 : 236 – 247 .[Crossref], [Web of Science ®] , [Google Scholar]).     

Restricted estimation and testing of hypothesis in linear measurement errors models

ABSTRACT

In this article, the linear models with measurement error both in the response and in the covariates are considered. Following Shalabh et al. (2007 Shalabh, Garg, G., Misra, N. (2007). Restricted regression estimation in measurement error models. Comput. Stat. Data Anal. 52:1149–1166.[Crossref], [Web of Science ®] , [Google Scholar], 2009 Shalabh, Garg, G., Misra, N. (2009). Use of prior information in the consistent estimation of regression coefficients in measurement error models. J. Multivariate Anal. 100:1498–1520.[Crossref], [Web of Science ®] , [Google Scholar]), we propose several restricted estimators for the regression coefficients. The consistency and asymptotic normality of the restricted estimators are established. Furthermore, we also discuss the superiority of the restricted estimators to unrestricted estimators under Pitman closeness criterion. We also develop several variance estimators and establish their asymptotic distributions. Wald-type statistics are constructed for testing the linear restrictions. Finally, Monte Carlo simulations are conducted to illustrate the finite-sample properties of the proposed estimators.     

Hazards Regression for Cancer Studies

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 length-biased 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 proportional hazards model. In this article, by modeling growth function as a function of covariates, we demonstrate 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. A simulation study is conducted to demonstrate the potential usefulness of the proposed estimators for the regression parameters in the proportional and additive hazards model.     

Inferences on partial linear models with right censored data by empirical likelihood

ABSTRACT

This article considers inference for partial linear models with right censored data. We use empirical likelihood based on the Buckley and James (1979 Buckley, J., James, I. (1979). Linear regression with censored data. Biometrika 66:429–436.[Crossref], [Web of Science ®] , [Google Scholar]) estimating equation to derive the confidence region for the regression parameter. We introduce an adjusted empirical likelihood ratio statistic for the parameter of interest and show that its limiting distribution is standard chi-square. A simulation is carried out to compare our method with the synthetic data approach in Wang and Li (2002 Wang, Q.-H., Li, G. (2002). Empirical Likelihood Semiparametric Regression Analysis under Random Censorship. J. Multivariate Anal. 83:469–486.[Crossref], [Web of Science ®] , [Google Scholar]).     

On Augmented Franses Tests for Seasonal Unit Roots

This article extends the results reported in del Barrio Castro, Osborn and Taylor (2012 del Barrio Castro, T., Osborn, D.R., Taylor, A. M.R. (2012). On augmented HEGY tests for seasonal unit roots. Econometric Theor. 18:1121–1143.[Crossref], [Web of Science ®] , [Google Scholar]) to the approach followed by Franses (1991a Franses, P. H. (1991a). Model selection and seasonality in time series. Tibergen Institute Series, 18. [Google Scholar],b Franses, P.H. (1991b). Seasonality, non-stationarity and the forecasting of monthly time series. Int. J. Forecast. 7:199–208.[Crossref], [Web of Science ®] , [Google Scholar]) to test for seasonal unit roots, providing the asymptotic representation to the seasonal unit roots tests proposed by Franses for a general number of seasons S.     

Convergence of weighted sums for sequences of pairwise NQD random variables

ABSTRACT

The authors discuss the convergence for weighted sums of pairwise negatively quadrant dependent (NQD) random variables and obtain some new results which extend and improve the result of Bai and Cheng (2000) Bai, Z.D., Cheng, P.E. (2000). Marcinkiewicz strong laws for linear statistics. Stat. Probab. Lett. 46:105–112.[Crossref], [Web of Science ®] , [Google Scholar]. In addition, we relax some restrictions of the conditions in their result. Some new methods are used in this article which differ from that of Bai and Cheng (2000) Bai, Z.D., Cheng, P.E. (2000). Marcinkiewicz strong laws for linear statistics. Stat. Probab. Lett. 46:105–112.[Crossref], [Web of Science ®] , [Google Scholar].     

On General Continuous Triangular and Two-sided Power Distributions

In this article, we directly introduce the continuous version of the general discrete triangular distributions (Kokonendji and Zocchi, 2010 Kokonendji, C.C., Zocchi, S.S. (2010). Extensions of discrete triangular distribution and boundary bias in kernel estimation for discrete functions. Statist. Probab. Lett. 80:1655–1662.[Crossref], [Web of Science ®] , [Google Scholar]). It is bounded and, in general, unimodal with pike. It contains thus a very useful class of two-sided power distributions (van Dorp and Kotz, 2002a Van Dorp, J.R., Kotz, S. (2002a). A novel extension of the triangular distribution and its parameter estimation. Statistician 51:1–17. [Google Scholar],b Van Dorp, J.R., Kotz, S. (2002b). The standard two-sided power distribution and its properties; with applications in financial engineering. Amer. Statistician 56:90–99.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar], 2003 Van Dorp, J.R., Kotz, S. (2003). Generalization of two-sided power distributions and their convolution. Commun. Statist. Theor. Meth. 32:1703–1723.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]). Moments, particular cases, limit distributions, and relations between parameters are straightforwardly derived.