On the joint distribution of placement statistics under progressive censoring and applications to precedence test

In reliability and survival analysis, comparison of two or more populations is an important problem. For example, while comparing a treatment group with a control group, one may be interested in determining whether the observations in the treatment group have a longer lifetime than those from the control group, that is, whether the treatment is effective or not. In such a study, it would be extremely valuable to make a decision based on early failures. In this paper, we consider independent progressively Type-II censored random samples from two populations with cumulative distribution function's (cdf) F(·)$F(·)$ and G(·)$G(·)$ respectively, and discuss a precedence test for testing the equality of the two distributions based on placements. For this purpose, we derive the joint distribution of the first l   placement statistics from the progressively censored sample from F(·)$F(·)$. We then derive the exact null distribution of the precedence test statistic which is simply the sum of the placements. We provide the rejection regions for fixed levels of significance and various sample sizes and different progressive censoring schemes.     

Robust signed-rank estimation and variable selection for semi-parametric additive partial linear models

A fully nonparametric model may not perform well or when the researcher wants to use a parametric model but the functional form with respect to a subset of the regressors or the density of the errors is not known. This becomes even more challenging when the data contain gross outliers or unusual observations. However, in practice the true covariates are not known in advance, nor is the smoothness of the functional form. A robust model selection approach through which we can choose the relevant covariates components and estimate the smoothing function may represent an appealing tool to the solution. A weighted signed-rank estimation and variable selection under the adaptive lasso for semi-parametric partial additive models is considered in this paper. B-spline is used to estimate the unknown additive nonparametric function. It is shown that despite using B-spline to estimate the unknown additive nonparametric function, the proposed estimator has an oracle property. The robustness of the weighted signed-rank approach for data with heavy-tail, contaminated errors, and data containing high-leverage points are validated via finite sample simulations. A practical application to an economic study is provided using an updated Canadian household gasoline consumption data.     

The Challenges in Developing an Online Applied Statistics Program: Lessons Learned at Penn State University

Numerous professional fields have an increasing need for individuals trained in statistics and other quantitative analysis techniques. Today there exists great potential to fulfill this need by providing opportunities through online learning. However, to provide a high-quality education for returning adult professionals seeking advanced degrees in applied statistics online, many challenges need to be overcome. Based on our experience developing Penn State University’s online program in applied statistics, we discuss the evolution of the program’s curriculum, recruitment and development of online faculty, and meeting the requirements of students as important areas that require consideration in the development of an online program. We also highlight program evaluation strategies employed to ensure innovation and improvement in online education as cornerstones to a program’s success.     

Mean Time to Failure of Weighted k-out-of-n: G Systems

The purpose of this article is to develop a Monte-Carlo simulation algorithm for computing mean time to failure (MTTF) of weighted-k-out-of-n:G and linear consecutive-weighted-k-out-of-n:G systems. Our algorithm is based on the use of appropriately defined stochastic process which represents the total weight of the system at time t. These stochastic processes are explicitly defined and used along with the ordered component lifetimes to simulate MTTF of the systems with weighted components.     

Power analysis of independence testing for three-way contingency tables of small sizes

The first aim of this paper is to introduce a modular test for the three-way contingency table (TT). The second aim is to describe the procedure of generating TT using the bar method. The third aim is on the one hand to suggest the measure of untruthfulness of H₀ and on the other hand to compare the quality of independence tests by using their power. Critical values for analyzed statistics were determined by simulating the Monte Carlo method.     

A Bayesian analysis for the Wilcoxon signed-rank statistic

A Bayesian analysis is provided for the Wilcoxon signed-rank statistic (T⁺). The Bayesian analysis is based on a sign-bias parameter φ on the (0, 1) interval. For the case of a uniform prior probability distribution for φ and for small sample sizes (i.e., 6 ? n ? 25), values for the statistic T⁺ are computed that enable probabilistic statements about φ. For larger sample sizes, approximations are provided for the asymptotic likelihood function P(T⁺|φ) as well as for the posterior distribution P(φ|T⁺). Power analyses are examined both for properly specified Gaussian sampling and for misspecified non Gaussian models. The new Bayesian metric has high power efficiency in the range of 0.9–1 relative to a standard t test when there is Gaussian sampling. But if the sampling is from an unknown and misspecified distribution, then the new statistic still has high power; in some cases, the power can be higher than the t test (especially for probability mixtures and heavy-tailed distributions). The new Bayesian analysis is thus a useful and robust method for applications where the usual parametric assumptions are questionable. These properties further enable a way to do a generic Bayesian analysis for many non Gaussian distributions that currently lack a formal Bayesian model.     

Semiparametric inference for estimating equations with nonignorably missing covariates

We consider statistical inference of unknown parameters in estimating equations (EEs) when some covariates have nonignorably missing values, which is quite common in practice but has rarely been discussed in the literature. When an instrument, a fully observed covariate vector that helps identifying parameters under nonignorable missingness, is available, the conditional distribution of the missing covariates given other covariates can be estimated by the pseudolikelihood method of Zhao and Shao [(2015), ‘Semiparametric pseudo likelihoods in generalised linear models with nonignorable missing data’, Journal of the American Statistical Association, 110, 1577–1590)] and be used to construct unbiased EEs. These modified EEs then constitute a basis for valid inference by empirical likelihood. Our method is applicable to a wide range of EEs used in practice. It is semiparametric since no parametric model for the propensity of missing covariate data is assumed. Asymptotic properties of the proposed estimator and the empirical likelihood ratio test statistic are derived. Some simulation results and a real data analysis are presented for illustration.