Reducing component estimation for varying coefficient models

The estimation problem for varying coefficient models has been studied by many authors. We consider the problem in the case that the unknown functions admit different degrees of smoothness. In this paper we propose a reducing component local polynomial method to estimate the unknown functions. It is shown that all of our estimators achieve the optimal convergence rates. The asymptotic distributions of our estimators are also derived. The established asymptotic results and the simulation results show that our estimators outperform the the existing two-step estimators when the coefficient functions admit different degrees of smoothness. We also develop methods to speed up the estimation of the model and the selection of the bandwidths.     

Bayes sequential estimation for a Poisson process under a LINEX loss function

In this paper, within the framework of a Bayesian model, we consider the problem of sequentially estimating the intensity parameter of a homogeneous Poisson process with a linear exponential (LINEX) loss function and a fixed cost per unit time. An asymptotically pointwise optimal (APO) rule is proposed. It is shown to be asymptotically optimal for the arbitrary priors and asymptotically non-deficient for the conjugate priors in a similar sense of Bickel and Yahav [Asymptotically pointwise optimal procedures in sequential analysis, in Proceedings of the Fifth Berkeley Symposium on Mathematical Statistics and Probability, Vol. 1, University of California Press, Berkeley, CA, 1967, pp. 401–413; Asymptotically optimal Bayes and minimax procedures in sequential estimation, Ann. Math. Statist. 39 (1968), pp. 442–456] and Woodroofe [A.P.O. rules are asymptotically non-deficient for estimation with squared error loss, Z. Wahrsch. verw. Gebiete 58 (1981), pp. 331–341], respectively. The proposed APO rule is illustrated using a real data set.     

Bounds on expected generalized order statistics

We establish the upper nonpositive and all the lower bounds on the expectations of generalized order statistics based on a given distribution function with the finite mean and central absolute moment of a fixed order. We also describe the distributions for which the bounds are attained. The methods of deriving the lower nonpositive (upper nonnegative) and lower nonnegative (upper nonpositive) bounds are totally different. The first one, the greatest convex minorant method is the combination of the Moriguti and well-known Hölder inequalities and the latter one is based on the maximization of some norm on the properly chosen convex set. The paper completes the results of Cramer et al. [Evaluations of expected generalized order statistics in various scale units. Appl Math. 2002;29:285–295].     

Fitting of self-exciting threshold autoregressive moving average nonlinear time-series model through genetic algorithm and development of out-of-sample forecasts

Self-exciting threshold autoregressive moving average (SETARMA) nonlinear time-series model is considered here. Sufficient conditions for invertibility and stationarity are derived. Parameter estimation algorithm is developed by employing real-coded genetic algorithm stochastic optimization procedure. A significant feature of the work done is that optimal out-of-sample forecasts up to three-step ahead and their forecast error variances are derived analytically. Relevant computer programs are written in statistical analysis system (SAS) and C. As an illustration, annual mackerel catch time-series data are considered. Forecast performance of the fitted model for hold-out data is evaluated by using Naive and Monte Carlo approaches. It is found that optimal out-of-sample forecast values are quite close to actual values and estimated variances are quite close to theoretical values. Superiority of the SETARMA model over the SETAR model for equal predictive ability through Diebold–Mariano test is also established.     

Alternative Factors for Control Chart Construction

Conflicting decision signals resulting from Shewhart Mean and Variability Control Charts produce undesirable consequences. Therefore, we construct “six-sigma” control charts by alternative methods to avoid problems associated with Control Charts yielding different signals. By adjusting methods for control chart construction for dispersion charts, we produce new reference tables. The new tables provide consistent signals and appear similar and familiar to users of historical tables.     

A Multivariate Synthetic Control Chart for Monitoring the Process Mean Vector of Skewed Populations Using Weighted Standard Deviations

This article proposes a multivariate synthetic control chart for skewed populations based on the weighted standard deviation method. The proposed chart incorporates the weighted standard deviation method into the standard multivariate synthetic control chart. The standard multivariate synthetic chart consists of the Hotelling's T ² chart and the conforming run length chart. The weighted standard deviation method adjusts the variance–covariance matrix of the quality characteristics and approximates the probability density function using several multivariate normal distributions. The proposed chart reduces to the standard multivariate synthetic chart when the underlying distribution is symmetric. In general, the simulation results show that the proposed chart performs better than the existing multivariate charts for skewed populations and the standard T ² chart, in terms of false alarm rates as well as moderate and large mean shift detection rates based on the various degrees of skewnesses.     

The computational aspects of multiple attribute sampling procedures

One method of controlling the quality of incoming lots is through attribute sampling. To simultaneously control several (possibly dependent) attributes, properly chosen single attribute sampling plans can be merged into a multiple attribute sampling plan. The general form of such a plan is given and various alternatives are discussed. The multinomial distribution is used to develop formulae necessary for an analysis of a multiple attribute plan. Due to the lengthy nature of the calculations involved, a computer algorithm is outlined.     

Determining Optimal Number of Samples for Constructing Multivariate Control Charts

Normally, an average run length (ARL) is used as a measure for evaluating the detecting performance of a multivariate control chart. This has a direct impact on the false alarm cost in Phase II. In this article, we first conduct a simulation study to calculate both in-control and out-of-control ARLs under various combinations of process shifts and number of samples. Then, a trade-off analysis between sampling inspection and false alarm costs is performed. Both the simulation results and trade-off analysis suggest that the optimal number of samples for constructing a multivariate control chart in Phase I can be determined.     

The Effect of Methods for Handling Missing Values on the Performance of the MEWMA Control Chart

This article is concerned with the effect of the methods for handling missing values in multivariate control charts. We discuss the complete case, mean substitution, regression, stochastic regression, and the expectation–maximization algorithm methods for handling missing values. Estimates of mean vector and variance–covariance matrix from the treated data set are used to build the multivariate exponentially weighted moving average (MEWMA) control chart. Based on a Monte Carlo simulation study, the performance of each of the five methods is investigated in terms of its ability to obtain the nominal in-control and out-of-control average run length (ARL). We consider three sample sizes, five levels of the percentage of missing values, and three types of variable numbers. Our simulation results show that imputation methods produce better performance than case deletion methods. The regression-based imputation methods have the best overall performance among all the competing methods.     

Multivariate Control Charts Based on Hybrid Novelty Scores

We propose a new nonparametric multivariate control chart that integrates a novelty score. The proposed control chart uses as its monitoring statistic a hybrid novelty score, calculated based on the distance to local observations as well as on the distance to the convex hull constructed by its neighbors. The control limits of the proposed control chart were established based on a bootstrap method. A rigorous simulation study was conducted to examine the properties of the proposed control chart under various scenarios and compare it with existing multivariate control charts in terms of average run length (ARL) performance. The simulation results showed that the proposed control chart outperformed both the parametric and nonparametric Hotelling's T ² control charts, especially in nonnormal situations. Moreover, experimental results with real semiconductor data demonstrated the applicability and effectiveness of the proposed control chart. To increase the capability to detect small mean shift, we propose an exponentially weighted hybrid novelty score control chart. Simulation results indicated that exponentially weighted hybrid score charts outperformed the hybrid novelty score based control charts.