COMPUTING MARGINAL EFFECTS IN THE BOX–COX MODEL

ABSTRACT

This paper considers computation of fitted values and marginal effects in the Box–Cox regression model. Two methods, 1 the “smearing” technique suggested by Duan (see Ref. [10] Duan, N. 1983. Smearing Estimate: A Nonparametric Retransformation Method. J. Amer. Statistical Assoc., 78: 605–610. [Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) and 2 direct numerical integration, are examined and compared with the “naive” method often used in econometrics.     

THE MEDIAN ESTIMATE OF THE AUTOREGRESSIVE LOCATION PARAMETER

ABSTRACT

The concept of generalized order statistics was introduced by Kamps (1995 Kamps , U. ( 1995 ) A Concept of Generalized Order Statistics . Germany : B. G. Teubner Stuttgart [Crossref] , [Google Scholar]) to unify several concepts that have been used in statistics such as order statistics, record values, and sequential order statistics. Estimation of the parameters of the Burr type XII distribution are obtained based on generalized order statistics. The maximum likelihood and Bayes methods of estimation are used for this purposes. The Bayes estimates are derived by using the approximation form of Lindley (1980 Lindley , D. V. ( 1980 ). Approximate Bayesian methods . J. Trabajos de Estadistica 31 : 223 – 237 .[Crossref] , [Google Scholar]). Estimation based on upper records from the Burr model is obtained and compared by using Monte Carlo simulation study. Our results are specialized to the results of AL-Hussaini and Jaheen (1992 AL-Hussaini , E. K. , Jaheen , Z. F. ( 1992 ). Bayesian estimation of the parameters, reliability and failure rate functions of the Burr type XII failure model . J. Statist. Comput. Simul. 41 : 31 – 40 .[Taylor &; Francis Online] , [Google Scholar]) which are based on ordinary order statistics.     

SECOND-ORDER ASYMPTOTICS FOR SCORE TESTS IN HETEROSKEDASTIC t REGRESSION MODELS

ABSTRACT

This paper develops corrected score tests for heteroskedastic t regression models, thus generalizing results by Cordeiro, Ferrari and Paula[1] Cordeiro, G.M., Ferrari, S.L.P. and Paula, G.A. 1993. Improved Score Tests for Generalized Linear Models. Journal of the Royal Statistical Society B, 55: 661–674.  [Google Scholar] and Cribari-Neto and Ferrari[2] Cribari-Neto, F. and Ferrari, S.L.P. 1995. Second-order Asymptotics for Score Tests in Generalised Linear Models. Biometrika, 82: 426–432. [Crossref], [Web of Science ®] , [Google Scholar] for normal regression models and by Ferrari and Arellano-Valle[3] Ferrari, S.L.P. and Arellano-Valle, R. 1996. Modified Likelihood Ratio and Score Tests in Linear Regression Models Using the t Distribution. Brazilian Journal of Probability and Statistics, 10: 15–33.  [Google Scholar] for homoskedastic t regression models. We present, in matrix notation, Bartlett-type correction formulae to improve score tests in this class of models. The corrected score statistics have a chi-squared distribution to order n ^?1, where n is the sample size. We apply our main result to a few special models and present simulation results comparing the performance of the usual score tests and their corrected versions.     

Flexible parametric modelling of the hazard function in breast cancer studies

In cancer research, study of the hazard function provides useful insights into disease dynamics, as it describes the way in which the (conditional) probability of death changes with time. The widely utilized Cox proportional hazard model uses a stepwise nonparametric estimator for the baseline hazard function, and therefore has a limited utility. The use of parametric models and/or other approaches that enables direct estimation of the hazard function is often invoked. A recent work by Cox et al. [6 Cox, C., Chu, H., Schneider, M. F. and Munoz, A. 2007. Parametric survival analysis and taxonomy of hazard functions for the generalized gamma distribution. Stat. Med., 26: 4352–4374. [Crossref], [PubMed], [Web of Science ®] , [Google Scholar]] has stimulated the use of the flexible parametric model based on the Generalized Gamma (GG) distribution, supported by the development of optimization software. The GG distribution allows estimation of different hazard shapes in a single framework. We use the GG model to investigate the shape of the hazard function in early breast cancer patients. The flexible approach based on a piecewise exponential model and the nonparametric additive hazards model are also considered.     

Random Effects,Fixed Effects and Hausman's Test for the Generalized Mixed Regressive Spatial Autoregressive Panel Data Model

This article suggests random and fixed effects spatial two-stage least squares estimators for the generalized mixed regressive spatial autoregressive panel data model. This extends the generalized spatial panel model of Baltagi et al. (2013 Baltagi, B. H., Egger, P., Pfaffermayr, M. (2013). A generalized spatial panel data model with random effects. Econometric Reviews 32:650–685.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]) by the inclusion of a spatial lag term. The estimation method utilizes the Generalized Moments method suggested by Kapoor et al. (2007 Kapoor, M., Kelejian, H. H., Prucha, I. R. (2007). Panel data models with spatially correlated error components. Journal of Econometrics 127(1):97–130.[Crossref], [Web of Science ®] , [Google Scholar]) for a spatial autoregressive panel data model. We derive the asymptotic distributions of these estimators and suggest a Hausman test a la Mutl and Pfaffermayr (2011 Mutl, J., Pfaffermayr, M. (2011). The Hausman test in a Cliff and Ord panel model. Econometrics Journal 14:48–76.[Crossref], [Web of Science ®] , [Google Scholar]) based on the difference between these estimators. Monte Carlo experiments are performed to investigate the performance of these estimators as well as the corresponding Hausman test.     

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.     

An extended Birnbaum–Saunders distribution: Theory,estimation, and applications

Abstract

Birnbaum and Saunders (1969a Birnbaum, Z.W., Saunders, S.C. (1969a). A new family of life distributions. J. Appl. Probab. 6:319–327.[Crossref], [Web of Science ®] , [Google Scholar]) pioneered a lifetime model which is commonly used in reliability studies. Based on this distribution, a new model called the gamma Birnbaum–Saunders distribution is proposed for describing fatigue life data. Several properties of the new distribution including explicit expressions for the ordinary and incomplete moments, generating and quantile functions, mean deviations, density function of the order statistics, and their moments are derived. We discuss the method of maximum likelihood and a Bayesian approach to estimate the model parameters. The superiority of the new model is illustrated by means of three failure real data sets. We also propose a new extended regression model based on the logarithm of the new distribution. The last model can be very useful to the analysis of real data and provide more realistic fits than other special regression models.     

On some dynamic generalized information measures for bivariate lifetimes

Several generalizations to the concept of Kullback-Leibler divergence measure and Kerridge inaccuracy measure are available in the literature. In a recent paper Kundu (Metrika, 78:415–435, 2015 Kundu, C. 2015. Generalized measures of information for truncated random variables. Metrika 78:415–35.[Crossref], [Web of Science ®] , [Google Scholar]) considered a generalized K-L divergence measure of order (α, β). Nath (Metrika, 13:123–135, 1968 Nath, P. 1968. Inaccuracy and coding theory. Metrika 13:123–35.[Crossref] , [Google Scholar]) has also proposed generalized inaccuracy measure of order α. Here we address the question of extending these measures to higher dimensions with reference to residual lifetimes. In the present work, the generalized divergence and inaccuracy measures are extended for conditional lifetimes of two components having possibly different ages. Several properties, including monotonicity, and bounds of these measures are obtained for conditional random variables. Moreover, we study the effect of (increasing) monotone transformation on these generalized measures.     

A Simple Method for Testing Variance Components in Unbalanced Nested Model

In this article, we consider the unbalanced case of the three fold nested random effects model under partial balance. The distributions of unweighted sums of squares are obtained first. Using the method of generalized p value introduced in Tsui and Weerahandi (1989 Tsui , K. , Weerahandi , S. ( 1989 ). Generalized p-values in significance testing of hypotheses in the presence of nuisance parameters . Journal of the American Statistical Association 84 : 602 – 607 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]), a new method is proposed for hypothesis tests involving functions of variance components. To evaluate the sizes of the generalized p value, a simulation study is conducted. The results indicate that the proposed method performs well under all examined conditions.     

On the Additive Weibull Distribution

We derive explicit expressions for the moments, incomplete moments, quantile function and generating function of the additive Weibull model pioneered by Xie and Lai (1995 Xie, M., Lai, C.D. (1995). Reliability analysis using an additive Weibull model with bathtub-shaped failure rate function. Reliab. Eng. Syst. Safety 52:87–93.[Crossref], [Web of Science ®] , [Google Scholar]), which is a quite flexible distribution for fitting lifetime data with bathtub-shaped failure rate function. In addition, we estimate the model parameters by maximum likelihood and determine the observed information matrix. The flexibility of the additive Weibull distribution is illustrated by means of one application to real data.     

A Generalized Time-Dependent Conditional Linear Model with Left-Truncated and Right-Censored Data

Consider the model φ(S(y | X)) = β(y) ^T X, where φ is a known link function, S(· | X) is the survival function of a response Y given a covariate X = (1, X, X ²,…, X ^p ), and β(y) is an unknown vector of time-dependent regression coefficients. The response Y is subject to left truncation and right censoring. We assume that given X, Y is independent of (C, T) where C and T are censoring and truncation variables with P(C ≥ T) = 1. In this article, with some modification of the assumptions in Lemmas 5 and 6 of Iglesias-Pérez and González-Manteiga (1999 Iglesias-Pérez , C. J. , González-Manteiga , W. G. ( 1999 ). Strong representation of a generalized product-limit estimator for truncated and censored data with some application . J. Nonparametric Statist. 10 : 213 – 244 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]), we present an almost sure representation for the generalized product-limit estimator (GPL) of S(y | X). Based on the GPL and the approach of Teodorescu et al. (2010 Teodorescu , B. , Keilegom , I. V. , Cao , R. ( 2010 ). Generalized time-dependent conditional linear models under left truncation and right censoring . Ann. Instit. Statist. Math. 62 : 465 – 485 .[Crossref], [Web of Science ®] , [Google Scholar]), a least squares estimator of β(y) is obtained and a bootstrap procedure is proposed to choose the optimum bandwidth.     

An Estimator for the Survival Function on the Basis of Observable Censoring Times

Abstract

Chiu [Chiu, S. N. (1999 Chiu, S. N. 1999. An unbiased estimator for the survival function of censored data. Commun. Statist. - Theory Meth., 28(9): 2249–2260. [Taylor & Francis Online], [Web of Science ®] , [Google Scholar]). An unbiased estimator for the survival function of censored data. Commun. Statist. - Theory Meth. 28(9):2249–2260.] proposed a nonparametric estimator for the survival function which is based on observable censoring times in the general censoring model. His estimator is less efficient than the Product-Limit estimator. Considering an informative censoring model this drawback can partially be overcome. This is shown by a nonparametric, uniformly consistent estimator based on observable censoring times within the simple Koziol–Green model. Some asymptotic properties of the new estimator are investigated and it is compared with the well-known ACL-estimator.     

Confidence Intervals for Linear Function of Variance Components in Two-Fold Nested Random Effects Model

In the article, we consider the unbalanced case of the two-way nested random effects model under partial balance. Using the method of generalized confidence intervals (GCIs) introduced in Weeranhandi (1993 Weeranhandi , S. ( 1993 ). Generalized confidence intervals . J. Amer. Statist. Assoc. 88 : 899 – 905 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar] 1995 Weeranhandi , S. ( 1995 ). Exact Statistical Methods for Data Analysis . New York : Springer-Verlag . [Google Scholar]), a new method is proposed for constructing confidence intervals on linear function of variance components. To compare the resulted intervals with the Modified Large Sample (MLS) intervals by Hernandez and Burdick (1993 Hernandez , R. P. , Burdick , R. K. ( 1993 ). Confidence intervals on the total variance in unbalanced two-fold nested designs . Biom. J. 35 : 515 – 522 .[Crossref] , [Google Scholar]), a simulation study is conducted. The results indicate that the proposed method performs better than the MLS method, especially for very unbalanced designs.     

On testing proportionality in the Cox regression model by Andersen's plot

Andersen's plot, a graphical method for testing the proportionality assumption in the Cox Regression Model (Cox, 1972 Cox, D.R. (1972). Regression models and life tables (with discussion). J. Royal Stat. Soc. Ser. B 34:187–220. [Google Scholar]), first proposed by Kay (1977 Kay, R. (1977). Proportional hazard regression models and the analysis of censored survival data. Appl. Stat. 26(3):227–237.[Crossref] , [Google Scholar]) and popularized by Andersen (1982 Andersen, P.K. (1982). Testing goodness of fit of Cox's regression and life model. Biometrics 38:67–77.[Crossref], [Web of Science ®] , [Google Scholar]), has been used widely in biomedical research to check the validity of applying this popular regression model in survival analysis. Our theoretical derivation and examples show that the theoretical basis of this method is flawed. The graphical method should not be used in testing the proportionality. Instead, formal analytical methods based on residuals such as Cox–Snell residual and martingale residual should be used in practice.     

A Bayesian approach to analyse overdispersed longitudinal count data

In this paper, we consider a model for repeated count data, with within-subject correlation and/or overdispersion. It extends both the generalized linear mixed model and the negative-binomial model. This model, proposed in a likelihood context [17 G. Molenberghs, G. Verbeke, and C.G.B. Demétrio, An extended random-effects approach to modeling repeated, overdispersion count data, Lifetime Data Anal. 13 (2007), pp. 457–511.[Web of Science ®] , [Google Scholar],18 G. Molenberghs, G. Verbeke, C.G.B. Demétrio, and A. Vieira, A family of generalized linear models for repeated measures with normal and conjugate random effects, Statist. Sci. 25 (2010), pp. 325–347. doi: 10.1214/10-STS328[Crossref], [Web of Science ®] , [Google Scholar]] is placed in a Bayesian inferential framework. An important contribution takes the form of Bayesian model assessment based on pivotal quantities, rather than the often less adequate DIC. By means of a real biological data set, we also discuss some Bayesian model selection aspects, using a pivotal quantity proposed by Johnson [12 V.E. Johnson, Bayesian model assessment using pivotal quantities, Bayesian Anal. 2 (2007), pp. 719–734. doi: 10.1214/07-BA229[Crossref], [Web of Science ®] , [Google Scholar]].     

Goodness-of-Fit Test and Parameter Estimation for a Proportional Odds Model of Random Censorship

Censored data arise naturally in a number of fields, particularly in problems of reliability and survival analysis. There are several types of censoring, in this article, we will confine ourselves to the right randomly censoring type. Recently, Ahmadi et al. (2010 Ahmadi , J. , Doostparast , M. , Parsian , A. ( 2010 ). Bayes estimation based on random censored data for some life time models under symmetric and asymmetric loss functions . Communcations in Statistics-Theory and Methods , 39 : 3058 – 3071 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) considered the problem of estimating unknown parameters in a general framework based on the right randomly censored data. They assumed that the survival function of the censoring time is free of the unknown parameter. This assumption is sometimes inappropriate. In such cases, a proportional odds (PO) model may be more appropriate (Lam and Leung, 2001 Lam , K. F. , Leung , T. L. ( 2001 ). Marginal likelihood estimation for proportional odds models with right censored data . Lifetime Data Analysis 7 : 39 – 54 .[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]). Under this model, in this article, point and interval estimations for the unknown parameters are obtained. Since it is important to check the adequacy of models upon which inferences are based (Lawless, 2003 Lawless , J. F. (2003). Statistical Models and Methods for Lifetime Data. , 2nd ed. New York : John Wiley & Sons. [Google Scholar], p. 465), two new goodness-of-fit tests for PO model based on right randomly censored data are proposed. The proposed procedures are applied to two real data sets due to Smith (2002 Smith , P. J. ( 2002 ). Analysis of Failure and Survival Data . London : Chapman & Hall, CRC . [Google Scholar]). A Monte Carlo simulation study is conducted to carry out the behavior of the estimators obtained.     

Local influence diagnostics for generalized linear mixed models with overdispersion

Since the seminal paper by Cook and Weisberg [9 R.D. Cook and S. Weisberg, Residuals and Influence in Regression, Chapman &; Hall, London, 1982. [Google Scholar]], local influence, next to case deletion, has gained popularity as a tool to detect influential subjects and measurements for a variety of statistical models. For the linear mixed model the approach leads to easily interpretable and computationally convenient expressions, not only highlighting influential subjects, but also which aspect of their profile leads to undue influence on the model's fit [17 E. Lesaffre and G. Verbeke, Local influence in linear mixed models, Biometrics 54 (1998), pp. 570–582. doi: 10.2307/3109764[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]]. Ouwens et al. [24 M.J.N.M. Ouwens, F.E.S. Tan, and M.P.F. Berger, Local influence to detect influential data structures for generalized linear mixed models, Biometrics 57 (2001), pp. 1166–1172. doi: 10.1111/j.0006-341X.2001.01166.x[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]] applied the method to the Poisson-normal generalized linear mixed model (GLMM). Given the model's nonlinear structure, these authors did not derive interpretable components but rather focused on a graphical depiction of influence. In this paper, we consider GLMMs for binary, count, and time-to-event data, with the additional feature of accommodating overdispersion whenever necessary. For each situation, three approaches are considered, based on: (1) purely numerical derivations; (2) using a closed-form expression of the marginal likelihood function; and (3) using an integral representation of this likelihood. Unlike when case deletion is used, this leads to interpretable components, allowing not only to identify influential subjects, but also to study the cause thereof. The methodology is illustrated in case studies that range over the three data types mentioned.     

Survival Weibull regression model for mismeasured outcomes

In some survival studies, the exact time of the event of interest is unknown, but the event is known to have occurred during a particular period of time (interval-censored data). If the diagnostic tool used to detect the event of interest is not perfectly sensitive and specific, outcomes may be mismeasured; a healthy subject may be diagnosed as sick and a sick one may be diagnosed as healthy. In such cases, traditional survival analysis methods produce biased estimates for the time-to-failure distribution parameters (Paggiaro and Torelli 2004 Paggiaro, A., and N. Torelli. 2004. The effect of classification errors in survival data analysis. Statistical Methods and Applications 13:213–25.[Crossref] , [Google Scholar]). In this context, we developed a parametric model that incorporates sensitivity and specificity into a grouped survival data analysis (a case of interval-censored data in which all subjects are tested at the same predetermined time points). Inferential aspects and properties of the methodology, such as the likelihood function and identifiability, are discussed in this article. Assuming known and non differential misclassification, Monte Carlo simulations showed that the proposed model performed well in the case of mismeasured outcomes; the estimates of the relative bias of the model were lower than those provided by the naive method that assumes perfect sensitivity and specificity. The proposed methodology is illustrated by a study related to mango tree lifetimes.     

Analysis of Gamma and Weibull lifetime data under a general censoring scheme and in the presence of covariates

We consider the problem of estimating the lifetime distributions of survival times subject to a general censoring scheme called “middle censoring”. The lifetimes are assumed to follow a parametric family of distributions, such as the Gamma or Weibull distributions, and is applied to cases when the lifetimes come with covariates affecting them. For any individual in the sample, there is an independent, random, censoring interval. We will observe the actual lifetime if the lifetime falls outside of this censoring interval, otherwise we only observe the interval of censoring. This censoring mechanism, which includes both right- and left-censoring, has been called “middle censoring” (see Jammalamadaka and Mangalam, 2003 Jammalamadaka, S. Rao, Mangalam, V. (2003). Nonparametric estimation for middle censored data. J. Nonparamet. Stat. 15(2):253–265.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]). Maximum-likelihood estimation of the parameters as well as their large-sample properties are studied under this censoring scheme, including the case when covariates are available. We conclude with an application to a dataset from Environmental Economics dealing with ContingentValuation of natural resources.     

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.