Cox regression analysis in presence of collinearity: an application to assessment of health risks associated with occupational radiation exposure

This paper considers the analysis of time to event data in the presence of collinearity between covariates. In linear and logistic regression models, the ridge regression estimator has been applied as an alternative to the maximum likelihood estimator in the presence of collinearity. The advantage of the ridge regression estimator over the usual maximum likelihood estimator is that the former often has a smaller total mean square error and is thus more precise. In this paper, we generalized this approach for addressing collinearity to the Cox proportional hazards model. Simulation studies were conducted to evaluate the performance of the ridge regression estimator. Our approach was motivated by an occupational radiation study conducted at Oak Ridge National Laboratory to evaluate health risks associated with occupational radiation exposure in which the exposure tends to be correlated with possible confounders such as years of exposure and attained age. We applied the proposed methods to this study to evaluate the association of radiation exposure with all-cause mortality.     

Two-stage Huber estimation

In this paper we propose a new robust estimator in the context of two-stage estimation methods directed towards the correction of endogeneity problems in linear models. Our estimator is a combination of Huber estimators for each of the two stages, with scale corrections implemented using preliminary median absolute deviation estimators. In this way we obtain a two-stage estimation procedure that is an interesting compromise between concerns of simplicity of calculation, robustness and efficiency. This method compares well with other possible estimators such as two-stage least-squares (2SLS) and two-stage least-absolute-deviations (2SLAD), asymptotically and in finite samples. It is notably interesting to deal with contamination affecting more heavily the distribution tails than a few outliers and not losing as much efficiency as other popular estimators in that case, e.g. under normality. An additional originality resides in the fact that we deal with random regressors and asymmetric errors, which is not often the case in the literature on robust estimators.     

A test for abrupt change in hazard regression models with Weibull baselines

We develop a likelihood ratio test for an abrupt change point in Weibull hazard functions with covariates, including the two-piece constant hazard as a special case. We first define the log-likelihood ratio test statistic as the supremum of the profile log-likelihood ratio process over the interval which may contain an unknown change point. Using local asymptotic normality (LAN) and empirical measure, we show that the profile log-likelihood ratio process converges weakly to a quadratic form of Gaussian processes. We determine the critical values of the test and discuss how the test can be used for model selection. We also illustrate the method using the Chronic Granulomatous Disease (CGD) data.     

A class of weighted normal distributions and its variants useful for inequality constrained analysis

This article develops a class of the weighted normal distributions for which the probability density function has the form of a product of a normal density and a weight function. The class constitutes marginal distributions obtained from various kinds of doubly truncated bivariate normal distributions. This class of distributions strictly includes the normal, skew–normal and two-piece skew–normal and is useful for selection modelling and inequality constrained normal mean analysis. Some distributional properties and Bayesian perspectives of the class are given. Probabilistic representation of the distributions is also given. The representation is shown to be straightforward to specify distribution and to implement computation, with output readily adapted for required analysis. Necessary theories and illustrative examples are provided.     

Number of Replications Required in Control Chart Monte Carlo Simulation Studies

Monte Carlo simulations have been used extensively in studying the performance of control charts. Researchers have used various numbers of replications in their studies, but almost none of them provided justifications for the number of replications used. Currently, there are no empirically based recommendations regarding the required number of replications to ensure accurate results. This research examined six recently published studies to develop recommendations for the minimum number of replications necessary to reproduce the reported results within a specified degree of accuracy. The results of this study indicated that using 10,000 replications was unnecessarily large and a smaller number of replications could be used to reproduce the target ARLs within the 2% error bands satisfying the modified Mundfrom's criteria. In many cases, only 5,000 replications or fewer were required. In general, the number of replications required to reproduce the target ARL decreased as the shift size increased. In addition, the results of this study provide general recommendations for the required number of replications to use in future SPC simulation studies.     

The estimation of treatment effect for the censored bivariate data under the univariate censoring

This paper establishes a nonparametric estimator for the treatment effect on censored bivariate data under unvariate censoring. This proposed estimator is based on the one from Lin and Ying(1993)'s nonparametric bivariate survival function estimator, which is itself a generalized version of Park and Park(1995)' quantile estimator. A Bahadur type representation of quantile functions were obtained from the marginal survival distribution estimator of Lin and Ying' model. The asymptotic property of this estimator is shown below and the simulation studies are also given     

Multivariate outliers and decompositions of mahalanobis distance

Two decompositions of the Mahalanobis distance are considered. These decompositions help to explain some reasons for the outlyingness of multivari- ate observations. They also provide a graphical tool for identifying outliers including those that have a large influence on the multiple correlation coefficient, Illustrative examples are given.     

A family of estimators of population variance in two-occasion rotation patterns

Abstract

In this article, we have considered the problem of estimation of population variance on current (second) occasion in two occasion successive (rotation) sampling. A class of estimators of population variance has been proposed and its asymptotic properties have been discussed. The proposed class of estimators is compared with the sample variance estimator when there is no matching from the previous occasion and the Singh et al. (2013) estimator. Optimum replacement policy is discussed. It has been shown that the suggested estimator is more efficient than the Singh et al. (2013) estimator and a usual unbiased estimator when there is no matching. An empirical study is carried out in support of the present study.     

A Class of Multivariate Bilateral Selection t Distributions and Its Properties

This article proposes a class of multivariate bilateral selection t distributions useful for analyzing non-normal (skewed and/or bimodal) multivariate data. The class is associated with a bilateral selection mechanism, and it is obtained from a marginal distribution of the centrally truncated multivariate t. It is flexible enough to include the multivariate t and multivariate skew-t distributions and mathematically tractable enough to account for central truncation of a hidden t variable. The class, closed under linear transformation, marginal, and conditional operations, is studied from several aspects such as shape of the probability density function, conditioning of a distribution, scale mixtures of multivariate normal, and a probabilistic representation. The relationships among these aspects are given, and various properties of the class are also discussed. Necessary theories and two applications are provided.     

Hypothesis Testing Problems in an Unbalanced Longitudinal Ophthalmology Study

ABSTRACT

The analysis of clustered data in a longitudinal ophthalmology study is complicated by correlations between repeatedly measured visual outcomes of paired eyes in a participant and missing observations due to the loss of follow-up. In the present article we consider hypothesis testing problems in an ophthalmology study, where eligible eyes are randomized to two treatments (when two eyes of a participant are eligible, the paired eyes are assigned to different treatments), and vision function outcomes are repeatedly measured over time. A large sample-based nonparametric test statistic and a nonparametric Bootstrap test analog are proposed for testing an interaction effect of two factors and testing an effect of a eye-specific factor within a level of the other person-specific factor on visual function outcomes. Both test statistics allow for missing observations, correlations between repeatedly measured outcomes on individual eyes, and correlations between repeatedly measured outcomes on both eyes of each participant. A simulation study shows that these proposed test statistics maintain nominal significance levels approximately and comparable powers to each other, as well as higher powers than the naive test statistic ignoring correlations between repeated bilateral measurements of both eyes in the same person. For illustration, we apply the proposed test statistics to the changes of visual field defect score in the Advanced Glaucoma Intervention Study.