Elliptical safety region plots for C pk

The process capability index C _pk is widely used when measuring the capability of a manufacturing process. A process is defined to be capable if the capability index exceeds a stated threshold value, e.g. C _pk >4/3. This inequality can be expressed graphically using a process capability plot, which is a plot in the plane defined by the process mean and the process standard deviation, showing the region for a capable process. In the process capability plot, a safety region can be plotted to obtain a simple graphical decision rule to assess process capability at a given significance level. We consider safety regions to be used for the index C _pk . Under the assumption of normality, we derive elliptical safety regions so that, using a random sample, conclusions about the process capability can be drawn at a given significance level. This simple graphical tool is helpful when trying to understand whether it is the variability, the deviation from target, or both that need to be reduced to improve the capability. Furthermore, using safety regions, several characteristics with different specification limits and different sample sizes can be monitored in the same plot. The proposed graphical decision rule is also investigated with respect to power.     

De l'unicité des estimateurs robustes en régression lorsque le paramètre d'échelle et le paramètre de la régression sont estimés simultanément

Several methods have been suggested to calculate robust M- and G-M -estimators of the regression parameter β and of the error scale parameter σ in a linear model. This paper shows that, for some data sets well known in robust statistics, the nonlinear systems of equations for the simultaneous estimation of β, with an M-estimate with a redescending ψ-function, and σ, with the residual median absolute deviation (MAD), have many solutions. This multiplicity is not caused by the possible lack of uniqueness, for redescending ψ-functions, of the solutions of the system defining β with known σ; rather, the simultaneous estimation of β and σ together creates the problem. A way to avoid these multiple solutions is to proceed in two steps. First take σ as the median absolute deviation of the residuals for a uniquely defined robust M-estimate such as Huber's Proposal 2 or the L₁-estimate. Then solve the nonlinear system for the M-estimate with σ equal to the value obtained at the first step to get the estimate of β. Analytical conditions for the uniqueness of M and G-M-estimates are also given.     

Estimation for semi-Markov models with partial observations via self-consistency equations

Nonparametric estimators of the survival function S(t) = P(T ≥ t) for a partially observed time variable T have been defined by several methods, in particular, by integral self-consistency equations. The author establishes explicit expressions of the estimators in an additive form and extend this approach to several cases: a left-truncated and right-censored variable and the left-censored or left-truncated sojourn times of a right-censored semi-Markov process. These estimators are always identical to the product-limit estimators if hazard functions may be defined.     

Conditional design of the EWMA median chart with estimated parameters

The exponentially weighted moving average (EWMA) chart is often designed assuming the process parameters are known. In practice, the parameters are rarely known and need to be estimated from Phase I samples. Different Phase I samples are used when practitioners construct their own control chart's limits, which leads to the “Phase I between-practitioners” variability in the in-control average run length (ARL) of control charts. The standard deviation of the ARL (SDARL) is a good alternative to quantify this variability in control charts. Based on the SDARL metric, the performance of the EWMA median chart with estimated parameters is investigated in this paper. Some recommendations are given based on the SDARL metric. The results show that the EWMA median chart requires a much larger amount of Phase I data in order to reduce the variation in the in-control ARL up to a reasonable level. Due to the limitation of the amount of the Phase I data, the suggested EWMA median chart is designed with the bootstrap method which provides a good balance between the in-control and out-of-control ARL values.     

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.     

Estimation and comparison of lognormal parameters in the presence of censored data

It is assumed that k(k?>?2) independent samples of sizes n _i (i?=?1, …, k) are available from k lognormal distributions. Four hypothesis cases (H ₁ – H ₄) are defined. Under H ₁, all k median parameters as well as all k skewness parameters are equal; under H ₂, all k skewness parameters are equal but not all k median parameters are equal; under H ₃, all k median parameters are equal but not all k skewness parameters are equal; under H ₄, neither the k median parameters nor the k skewness parameters are equal. The Expectation Maximization (EM) algorithm is used to obtain the maximum likelihood (ML) estimates of the lognormal parameters in each of these four hypothesis cases. A (2k???1) degree polynomial is solved at each step of the EM algorithm for the H ₃ case. A two-stage procedure for testing the equality of the medians either under skewness homogeneity or under skewness heterogeneity is also proposed and discussed. A simulation study was performed for the case k?=?3.     

A new generation of process capability indices based on fuzzy measurements

Process capability indices (PCIs) provide numerical measures on whether a process conforms to the defined manufacturing capability prerequisite. These have been successfully applied by companies to compete with and to lead high-profit markets by evaluating the quality and productivity performance. The PCI C_p compares the output of a process to the specification limits (SLs) by forming the ratio of the width between the process SLs with the width of the natural tolerance limits which is measured by six process standard deviation units. As another common PCI, C_pm incorporates two variation components which are variation to the process mean and deviation of the process mean from the target. A meaningful generalized version of above PCIs is introduced in this paper which is able to handle in a fuzzy environment. These generalized PCIs are able to measure the capability of a fuzzy-valued process in producing products on the basis of a fuzzy quality. Fast computing formulas for the generalized PCIs are computed for normal and symmetric triangular fuzzy observations, where the fuzzy quality is defined by linear and exponential fuzzy SLs. A practical example is presented to show the performance of proposed indices.     

A Monte Carlo method for filtering a marked doubly stochastic Poisson process

The author provides an approximated solution for the filtering of a state-space model, where the hidden state process is a continuous-time pure jump Markov process and the observations come from marked point processes. Each state k corresponds to a different marked point process, defined by its conditional intensity function λ _k (t). When a state is visited by the hidden process, the corresponding marked point process is observed. The filtering equations are obtained by applying the innovation method and the integral representation theorem of a point process martingale. Since the filtering equations belong to the family of Kushner–Stratonovich equations, an iterative solution is calculated. The theoretical solution is approximated and a Monte Carlo integration technique employed to implement it. The sequential method has been tested on a simulated data set based on marked point processes widely used in the statistical analysis of seismic sequences: the Poisson model, the stress release model and the Etas model.     

Convergence of non-linear functionals of smoothed empirical processes and kernel density estimates

We consider regularizations by convolution of the empirical process and study the asymptotic behaviour of non-linear functionals of this process. Using a result for the same type of non-linear functionals of the Brownian bridge, shown in a previous paper [4], and a strong approximation theorem, we prove several results for the p-deviation in estimation of the derivatives of the density. We also study the asymptotic behaviour of the number of crossings of the smoothed empirical process defined by Yukich [17] and of a modified version of the Kullback deviation.     

Optimal Statistical Design of a Multivariate EWMA Chart Based on ARL and MRL

Statistical design is applied to a multivariate exponentially weighted moving average (MEWMA) control chart. The chart parameters are control limit H and smoothing constant r. The choices of the parameters depend on the number of variables p and the size of the process mean shift δ. The MEWMA statistic is modeled as a Markov chain and the Markov chain approach is used to determine the properties of the chart. Although average run length has become a traditional measure of the performance of control schemes, some authors have suggested other measures, such as median and other percentiles of the run length distribution to explain run length properties of a control scheme. This will allow a thorough study of the performance of the control scheme. Consequently, conclusions based on these measures would provide a better and comprehensive understanding of a scheme. In this article, we present the performance of the MEWMA control chart as measured by the average run length and median run length. Graphs are given so that the chart parameters of an optimal MEWMA chart can be determined easily.     

Modified repeated median filters

We discuss moving window techniques for fast extraction of a signal composed of monotonic trends and abrupt shifts from a noisy time series with irrelevant spikes. Running medians remove spikes and preserve shifts, but they deteriorate in trend periods. Modified trimmed mean filters use a robust scale estimate such as the median absolute deviation about the median (MAD) to select an adaptive amount of trimming. Application of robust regression, particularly of the repeated median, has been suggested for improving upon the median in trend periods. We combine these ideas and construct modified filters based on the repeated median offering better shift preservation. All these filters are compared w.r.t. fundamental analytical properties and in basic data situations. An algorithm for the update of the MAD running in time O(log n) for window width n is presented as well.     

La méthode de vraisemblance pour les processus linéaires spatiaux définis sur un treillis triangulaire

Spatial linear processes {X_s, s ? T} where T is a triangular lattice in R² are considered. Special attention is given to the class of spatial moving-average processes. Precisely, for each site s T, the variable X_s is defined as a linear combination of real-valued random shocks located at the vertices of regular concentric hexagons centered at s. For Gaussian random shocks, the process is also Gaussian, and estimates of its parameters are obtained by maximizing the exact likelihood. For non-Gaussian random shocks, the exact likelihood is difficult to obtain; however, the Gaussian likelihood is still used giving the pseudo-Gaussian likelihood estimates. The behaviour of these estimates is analyzed through the study of asymptotic properties and some simulation experiments based on an isotropic model defined with one coefficient.     

A Phase I nonparametric Shewhart-type control chart based on the median

A nonparametric Shewhart-type control chart is proposed for monitoring the location of a continuous variable in a Phase I process control setting. The chart is based on the pooled median of the available Phase I samples and the charting statistics are the counts (number of observations) in each sample that are less than the pooled median. An exact expression for the false alarm probability (FAP) is given in terms of the multivariate hypergeometric distribution and this is used to provide tables for the control limits for a specified nominal FAP value (of 0.01, 0.05 and 0.10, respectively) and for some values of the sample size (n) and the number of Phase I samples (m). Some approximations are discussed in terms of the univariate hypergeometric and the normal distributions. A simulation study shows that the proposed chart performs as well as, and in some cases better than, an existing Shewhart-type chart based on the normal distribution. Numerical examples are given to demonstrate the implementation of the new chart.     

A mean deviation-based approach to monitor process variability

The study proposes a Shewhart-type control chart, namely an MD chart, based on average absolute deviations taken from the median, for monitoring changes (especially moderate and large changes – a major concern of Shewhart control charts) in process dispersion assuming normality of the quality characteristic to be monitored. The design structure of the proposed MD chart is developed and its comparison is made with those of two well-known dispersion control charts, namely the R and S charts. Using power curves as a performance measure, it has been observed that the design structure of the proposed MD chart is more powerful than that of the R chart and is very close competitor to that of the S chart, in terms of discriminatory power for detecting shifts in the process dispersion. The non-normality effect is also examined on design structures of the three charts, and it has been observed that the design structure of the proposed MD chart is least affected by departure from normality.     

Globally applicable control chart for online monitoring of stability of process mean

The shape features of run chart patterns of the most recent m observations arising from stable and unstable processes are different. Using this fact, a new monitoring statistic is defined whose value for given m depends on the pattern parameters but not on the process parameters. A control chart for this statistic for given m, therefore, will be globally applicable to normal processes. The simulation study reveals that the proposed statistic approximately follows normal distribution. The performances of the globally applicable control chart in terms of average run lengths (ARLs) are evaluated and compared with the X chart. Both in-control ARL and out-of-control ARLs with respect to different abnormal process conditions are found to be larger than the X chart. However, the proposed concept is promising because it can eliminate the burden of designing separate control charts for different quality characteristics or processes in a manufacturing set-up.     

La médiane simpliciale d'Oja: Existence,unicité et stabilité

Oja (1983) examined various ways of measuring location, scatter, skewness, and kurtosis for multivariate distributions. Among other measures of location, he introduced a generalised median known in this paper under the name of the Oja median. In our study of the existence of that median, we show that Oja's definition can only be applied to distributions having a mean. In dimension d θ 2, we establish that the usual method of extension breaks down, which raises the question of the validity of the concept as a notion of median. Two fundamental theoretical properties of that median are also considered: uniqueness and consistency.     

A J‐function for Inhomogeneous Spatio‐temporal Point Processes

We propose a new summary statistic for inhomogeneous intensity‐reweighted moment stationarity spatio‐temporal point processes. The statistic is defined in terms of the n‐point correlation functions of the point process, and it generalizes the J‐function when stationarity is assumed. We show that our statistic can be represented in terms of the generating functional and that it is related to the spatio‐temporal K‐function. We further discuss its explicit form under some specific model assumptions and derive ratio‐unbiased estimators. We finally illustrate the use of our statistic in practice. © 2014 Board of the Foundation of the Scandinavian Journal of Statistics     

On Semiparametric Mode Regression Estimation

It has been found that, for a variety of probability distributions, there is a surprising linear relation between mode, mean, and median. In this article, the relation between mode, mean, and median regression functions is assumed to follow a simple parametric model. We propose a semiparametric conditional mode (mode regression) estimation for an unknown (unimodal) conditional distribution function in the context of regression model, so that any m-step-ahead mean and median forecasts can then be substituted into the resultant model to deliver m-step-ahead mode prediction. In the semiparametric model, Least Squared Estimator (LSEs) for the model parameters and the simultaneous estimation of the unknown mean and median regression functions by the local linear kernel method are combined to infer about the parametric and nonparametric components of the proposed model. The asymptotic normality of these estimators is derived, and the asymptotic distribution of the parameter estimates is also given and is shown to follow usual parametric rates in spite of the presence of the nonparametric component in the model. These results are applied to obtain a data-based test for the dependence of mode regression over mean and median regression under a regression model.     

Locking the Wiener Process to its Level-Crossing Time

We consider the specific transformation of a Wiener process {X(t), t ≥ 0} in the presence of an absorbing barrier a that results when this process is “time-locked” with respect to its first passage time T _a through a criterion level a, and the evolution of X(t) is considered backwards (retrospectively) from T _a . Formally, we study the random variables defined by Y(t) ≡ X(T _a  ? t) and derive explicit results for their density and mean, and also for their asymptotic forms. We discuss how our results can aid interpretations of time series “response-locked” to their times of crossing a criterion level.     

Semiparametric median residual life model and inference

For randomly censored data, the authors propose a general class of semiparametric median residual life models. They incorporate covariates in a generalized linear form while leaving the baseline median residual life function completely unspecified. Despite the non‐identifiability of the survival function for a given median residual life function, a simple and natural procedure is proposed to estimate the regression parameters and the baseline median residual life function. The authors derive the asymptotic properties for the estimators, and demonstrate the numerical performance of the proposed method through simulation studies. The median residual life model can be easily generalized to model other quantiles, and the estimation method can also be applied to the mean residual life model. The Canadian Journal of Statistics 38: 665–679; 2010 © 2010 Statistical Society of Canada