A Multivariate Adaptive Exponentially Weighted Moving Average Control Chart

A multivariate extension of the adaptive exponentially weighted moving average (AEWMA) control chart is proposed. The new multivariate scheme can detect small and large shifts in the process mean vector effectively. The proposed scheme can be viewed as a smooth combination of a multivariate exponentially weighted moving average (MEWMA) chart and a Shewhart χ²-chart. The optimal design of the proposed chart is given according to a pre-specified in-control average run length and two shift sizes; a small and large shift each measured in terms of the non centrality parameter. The signal resistance of the newly proposed multivariate chart is also given. Comparisons among the new chart, the MEWMA chart, and the combined Shewhart-MEWMA (S-MEWMA) chart in terms of the standard and worst-case average run length profiles are presented. In addition, the three charts are compared with respect to their worst-case signal resistance values. The proposed chart gives somewhat better worst-case ARL and signal resistance values than the competing charts. It also gives better standard ARL performance especially for moderate and large shifts. The effectiveness of our proposed chart is illustrated through an example with simulated data set.     

A New Goodness-of-Fit Test Based on the Empirical Characteristic Function

In this article, we present a goodness-of-fit test for a distribution based on some comparisons between the empirical characteristic function c_n(t) and the characteristic function of a random variable under the simple null hypothesis, c₀(t). We do this by introducing a suitable distance measure. Empirical critical values for the new test statistic for testing normality are computed. In addition, the new test is compared via simulation to other omnibus tests for normality and it is shown that this new test is more powerful than others.     

Inverse Regression in Binary Response LDV Model

Sliced inverse regression, a link-free and distribution-free method, is applied to binary response limited dependent variable models. An inverse regression property of binary response LDV model is found. Based on this property, if the distributions of X _j (j = 1, 2,…, p) satisfy the linearity condition, then β can be estimated up to a positive multiplicative scalar without any assumptions on the distribution of error ε. Moreover, the estimator can be proved to be asymptotically normal based on which testing hypotheses are considered. Simulations results are reported.     

Generalized Linear Failure Rate Distribution

The exponential and Rayleigh are the two most commonly used distributions for analyzing lifetime data. These distributions have several desirable properties and nice physical interpretations. Unfortunately, the exponential distribution only has constant failure rate and the Rayleigh distribution has increasing failure rate. The linear failure rate distribution generalizes both these distributions which may have non increasing hazard function also. This article introduces a new distribution, which generalizes linear failure rate distribution. This distribution generalizes the well-known (1) exponential distribution, (2) linear failure rate distribution, (3) generalized exponential distribution, and (4) generalized Rayleigh distribution. The properties of this distribution are discussed in this article. The maximum likelihood estimates of the unknown parameters are obtained. A real data set is analyzed and it is observed that the present distribution can provide a better fit than some other very well-known distributions.     

Estimation of Survivorship Function Based on Doubly Censored Data with Application to HIV/AIDS

Left censoring concept has been defined in different ways in statistical applications. Turnbull (1974 Turnbull , B. W. ( 1974 ). Nonparametric estimation of a survivorship function with doubly censored data . J. Amer. Statist. Assoc. 69 : 169 – 173 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) defines it in a particular way. Whereas in recent literature, especially in epidemiological studies, it has been defined in another way. This difference between the two approaches is the main reason that despite simplicity, Turnbull method cannot be applicable in all cases of doubly censored data. In this article we present a modified Turnbull method for analysis of doubly censored data adequate with recent definition. Comparison has been done with other statistical methods, including imputation estimator, full likelihood-based and conditional likelihood-based approach using Iranian HIV data.     

Inference for a Simple Step-Stress Model with Type-I Censoring and Lognormally Distributed Lifetimes

Accelerated life-testing (ALT) is a very useful technique for examining the reliability of highly reliable products. It allows the experimenter to obtain failure data more quickly at increased stress levels than under normal operating conditions. A step-stress model is one special class of ALT, and in this article we consider a simple step-stress model under the cumulative exposure model with lognormally distributed lifetimes in the presence of Type-I censoring. We then discuss inferential methods for the unknown parameters of the model by the maximum likelihood estimation method. Some numerical methods, such as the Newton–Raphson and quasi-Newton methods, are discussed for solving the corresponding non-linear likelihood equations. Next, we discuss the construction of confidence intervals for the unknown parameters based on (i) the asymptotic normality of the maximum likelihood estimators (MLEs), and (ii) parametric bootstrap resampling technique. A Monte Carlo simulation study is carried out to examine the performance of these methods of inference. Finally, a numerical example is presented in order to illustrate all the methods of inference developed here.     

Renyi Entropy Properties of Order Statistics

In this article, we discuss some properties of Renyi entropy and Renyi information of order statistics. Some bounds for Renyi entropy of order statistics are obtained. Also, we relate Renyi entropy ordering of order statistics to Renyi entropy ordering and other well known orderings of parent random variables. Then it is proved that the Renyi information between order statistics and parent random variable is distribution free, and it is shown, as expected, the distance is minimum for the median.     

Free Knot Splines with RJMCMC in Survival Data Analysis

The B-spline representation is a common tool to improve the fitting of smooth nonlinear functions, it offers a fitting as a piecewise polynomial. The regions that define the pieces are separated by a sequence of knots. The main difficulty in this type of modeling is the choice of the number and the locations of these knots. The Reversible Jump Markov Chain Monte Carlo (RJMCMC) algorithm provides a solution to simultaneously select these two parameters by considering the knots as free parameters. This algorithm belongs to the MCMC techniques that allow simulations from target distributions on spaces of varying dimension. The aim of the present investigation is to use this algorithm in the framework of the analysis of survival time, for the Cox model in particular. In fact, the relation between the hazard ratio function and the covariates being assumed to be log-linear, this assumption is too restrictive. Thus, we propose to use the RJMCMC algorithm to model the log hazard ratio function by a B-spline representation with an unknown number of knots at unknown locations. This method is illustrated with two real data sets: the Stanford heart transplant data and lung cancer survival data. Another application of the RJMCMC is selecting the significant covariates, and a simulation study is performed.     

Minimum and Maximum Information Censoring Plans in Progressive Censoring

Expressions for the entropy, the Kullback-Leibler information, and the I _α-information are established for distributions of progressively Type-II censored order statistics. These results are used to identify minimum and maximum information censoring plans. In particular, we find minimum and maximum entropy plans for DFR, exponential, Pareto, reflected power, and Weibull distributions. The results for Kullback-Leibler and I _α-information hold for any continuous distribution.     

A Weighting Approach for GEE Analysis with Missing Data

We propose a new weighting (WT) method to handle missing categorical outcomes in longitudinal data analysis using generalized estimating equations (GEE). The proposed WT provides a valid GEE estimator when the data are missing at random (MAR), and has more stable weights and shows advantage in efficiency compared to the inverse probability weighing method in the presence of small observation probabilities. The WT estimator is similar to the stabilized weighting (SWT) estimator under mild conditions, but it is more stable and efficient than SWT when the associations of the outcome with the observation probabilities and the covariate are strong.