Next Article in Journal
On a New Approach for Stability and Controllability Analysis of Functional Equations
Next Article in Special Issue
New Machine-Learning Control Charts for Simultaneous Monitoring of Multivariate Normal Process Parameters with Detection and Identification
Previous Article in Journal
An Enhanced Multi-Constraint Optimization Algorithm for Efficient Network Topology Generation
Previous Article in Special Issue
A Three-Stage Nonparametric Kernel-Based Time Series Model Based on Fuzzy Data
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Process Capability and Performance Indices for Discrete Data

by
Vasileios Alevizakos
Department of Mathematics, National Technical University of Athens, Zografou, 15773 Athens, Greece
Mathematics 2023, 11(16), 3457; https://doi.org/10.3390/math11163457
Submission received: 17 July 2023 / Revised: 29 July 2023 / Accepted: 1 August 2023 / Published: 9 August 2023

Abstract

:
Process capability and performance indices (PCIs and PPIs) are used in industry to provide numerical measures for the capability and performance of several processes. The majority of the literature refers to PCIs and PPIs for continuous data. The aim of this paper is to compute the classical indices for discrete data following Poisson, binomial or negative binomial distribution using various transformation techniques. A simulation study under different situations of a process and comparisons with other existing PCIs for discrete data are also presented. The methodology of computing the indices is easy to use, and as a result, one can have an assessment of the process capability and performance without difficulty. Three examples are further provided to illustrate the application of the transformation techniques.

1. Introduction

Process capability indices (PCIs) are widely used in industry as they are considered practical tools for continuous quality improvement. They provide numerical measures of the process ability to produce output within specification limits on a quality characteristic of interest. The four most broadly used PCIs are C p (Juran [1]), C p k (Kane [2]), C p m (Chan, Cheng, and Spiring [3], Hsiang and Taguchi [4]) and C p m k (Choi and Owen [5], Pearn, Kotz, and Johnson [6]) and are defined as follows:
C p = U L 6 σ ,
C p k = m i n { U μ , μ L } 3 σ ,
C p m = U L 6 σ 2 + ( μ T ) 2 ,
C p m k = m i n { U μ , μ L } 3 σ 2 + ( μ T ) 2 ,
where U and L are the upper and the lower specification limits, respectively; T is the target value of the process; μ is the process mean; and σ is the process standard deviation. If μ and σ are unknown, one can replace them with the sample mean X ¯ and the sample standard deviation S, respectively. The properties, the estimation methods and the confidence limits of PCIs have been examined by many authors, such as Rodriguez [7], Kotz and Johnson [8,9], Kotz and Lovelace [10], and Pearn and Kotz [11].
In 1991, the Automotive Industry Action Group (AIAG) recommended using the PCIs C p and C p k when the process is in control with the process standard deviation estimated using the formula σ ^ = R ¯ / d 2 (s-short-term). On the other hand, when the process is not in control, AIAG recommended using the process performance indices (PPIs) P p and P p k , given by
P p = U L 6 s ,
P p k = m i n { U μ , μ L } 3 s ,
where s is computed using s = i = 1 n ( X i X ¯ ) 2 n 1 (s-long-term). The main difference between PCIs and PPIs is that PCIs are used to predict the capability of a process to produce parts conforming to specifications, while PPIs are used to evaluate the behaviour of a process. Note that for a stable process, the differences between C p ( C p k ) and P p ( P p k ) are negligible (Montgomery [12]).
The use of the four well-known PCIs presupposes that the process is in control as well as that the quality characteristic of interest is normally distributed. However, there are many processes in which the distribution of the quality characteristic is continuous but non-normal. Many authors have studied PCIs on continuous non-normal distributions. Clements [13] used a percentile method (known as “Clements’ method”) in order to calculate the C p and C p k indices for non-normal Pearsonian populations. Pearn and Kotz [14] applied Clements’ method for computing the C p m and C p m k indices. Rivera, Hubele and Lawrence [15] studied the effect of several transformation methods on the estimate of C p k index for non-normal data. Pearn and Chen [16] proposed a modification of Clements’ method in order to calculate the four well-known PCIs for non-normal Pearsonian populations. Castagliola [17] proposed a method based on Burr’s distributions to estimate the C p and C p k indices independently of the real process distribution. Hosseinifard et al. [18] proposed a root transformation technique to estimate PCIs of right skewed processes. Hosseinifard, Abbasi and Niaki [19] studied and compared various methods for estimating PCIs of non-normal processes.
In practice, there are many processes that the quality characteristic of interest is discrete and, unfortunately, the above studies cannot be used. Few studies have been focused on PCIs for discrete quality characteristics. Bothe [20] computed the C p k index based on the percentage of nonconforming parts of a process with attribute data. Yeh and Bhattcharya [21] proposed a PCI that is related to the probability of non-conformance of a process and can be applied to continuous as well as discrete distributions. Borges and Ho [22] proposed a PCI based on the process fraction defective. Perakis and Xekalaki [23,24] proposed a PCI based on the proportion of conformance of the process, which can be used for continuous or discrete processes. Maiti, Saha and Nanda [25] defined a generalized index, named C p y , based on the ratio of the proportion of specification conformance to the proportion of desired conformance, which can be used for continuous as well as discrete processes. Moreover, they introduced two modifications of the C p y index, the C p y k and C p T k indices, for studying off-centered and off-target processes, respectively. Maravelakis [26] proposed the Q transformation technique to estimate the four well-known PCIs for Poisson and binomial data. Dey and Saha [27] considered three bootstrap confidence intervals of the C p y k index using different methods of estimation for the logistic exponential distribution. Pal and Gauri [28] studied several approaches for estimating the PCIs C p u and C p l of a binomial process. In addition, Gauri and Pal [29] assessed the relative goodness of some generalized PCIs in quantifying the capability of a Poisson or binomial process.
Motivated by the work of Maravelakis [26] and of Pal and Gauri [28], in this paper, we expand the transformation techniques for Poisson and binomial data and we also apply them to negative binomial data in order to compute the four classical PCIs and the two PPIs. Moreover, we compare these indices with other existing PCIs for discrete data to make conclusions about the proposed methodology. Finally, we illustrate the application of the transformation techniques to real data from Poisson, binomial and negative binomial processes.
The paper is organized as follows. In Section 2, we present the transformation techniques for Poisson, binomial and negative binomial data, while a brief overview of existing PCIs for discrete data is given in Section 3. Section 4 presents the results of a simulation study that was performed to compare the indices using different transformation methods. In Section 5, we present three illustrative examples, one for each discrete distribution. Finally, we conclude the paper in Section 6.

2. Transformation Techniques to Compute PCIs

In this section, we present several transformation techniques for Poisson, binomial and negative binomial data, in order to transform them into being approximately independent and normally distributed.

2.1. Transformation Techniques for Poisson Data

Let X be the number of defects or nonconformities in an inspection unit with a standard size. If the process is stable, i.e., the defects rate c is constant, and each inspection is performed independently, then the variable X follows the Poisson distribution with probability mass function (pmf) given by
f ( x ; c ) = e c c x x ! , x = 0 , 1 , 2 , ,
where c > 0 is the parameter of the Poisson distribution. The mean and the variance of X are equal to c.

2.1.1. Anscombe’s Transformation

Let x i , i = 1 , 2 , be observations of the Poisson variable X. Bartlett [30] used two nonlinear transformations based on the square root in order to transform Poisson data into normal data. Anscombe [31] proposed a variation of Balett’s transformation afterwards, and he demonstrated that the statistics y i = x i + 3 / 8 , i = 1 , 2 , are approximately normally distributed with mean and variance given by c + 3 / 8 1 / ( 8 c ) and 0.25 , respectively. In this paper, we use a modification of this transformation in order to obtain a variance close to 1, that is
y i = 2 x i + 3 8 , for i = 1 , 2 ,

2.1.2. Freeman–Tukey’s Transformation

Freeman and Tukey [32] proposed another square root transformation for approximately normalizing Poisson data. They displayed that the statistics
y i = x i + x i + 1 , i = 1 , 2 ,
are approximately normal with mean and variance given by 2 c and 1, respectively. Ryan and Schwertman [33] applied this transformation to the classical Shewhart c-chart, so as to study its performance as well as to investigate the optimal control limits.

2.1.3. The Q Transformation

Quesenberry [34,35] defined the Q statistic and introduced the Poisson Q-chart for data with known and unknown parameter c. For c = c 0 known a priori, define
u i = e c 0 k = 0 x i c 0 k ( k ! ) , Q i = Φ 1 ( u i ) , for i = 1 , 2 ,
where Φ 1 ( · ) is the inverse function of the standard normal distribution function. Then, the values Q 1 , Q 2 , are approximately standard normal.
If the defects rate c is unknown, then we consider that the data are in a sequence of values ( n i , x i ) , for i = 1 , 2 , Let t i = j = 1 i x j , N i = j = 1 i n j and define
u i = k = 0 x i t i k n i N i k 1 n i N i n i k Q i = Φ 1 ( u i ) , for i = 2 , 3 ,
Then, the values Q 2 , Q 3 , are approximately standard normal. The Poisson Q-chart is obtained by plotting the statistics Q with control limits UCL = 3, CL = 0 and LCL = −3. We notice that for the case where c is unknown, the first plotting statistic is the Q 2 . Moreover, according to Quesenberry [36], this approximation performs satisfactorily after the first observations, as in an effort to approximate c, past observations are used to update the Q statistic.

2.2. Transformation Techniques for Binomial Data

Let X be the number of noncomforming products found when n products are inspected. If the process is stable, i.e., the probability p that any product does not conform to specifications is constant, and each inspection is performed independently, then the variable X follows the binomial distribution with parameters n and p and the pmf given by
f ( x ; n , p ) = n x p x ( 1 p ) n x , x = 0 , 1 , , n
The mean and the variance of X are equal to n p and n p ( 1 p ) , respectively.

2.2.1. Freeman–Tukey’s Transformation

Let x i be the number of nonconforming items produced in a sample of size n with known and constant probability p = p 0 of observing a nonconforming product. Freeman and Tukey [32] used the average angular transformation to transform binomial data to normal. The transformation is defined as
y i = sin 1 x i n + 1 + sin 1 x i + 1 n + 1 , for i = 1 , 2 ,
The statistics y i are normally distributed with mean and variance equal to s i n 1 p and 1 / ( n + 0.5 ) , respectively. Ryan [37] used this transformation to construct a control chart for monitoring p.

2.2.2. Chen’s Transformation

Johnson and Kotz [38] showed that the statistics
y i = sin 1 x i + 3 8 n + 3 4 for i = 1 , 2 ,
are approximately normally distributed with mean and variance equal to s i n 1 p and 1 / ( 4 n ) , respectively. Ryan [37] used this transformation to construct a control chart for monitoring p. Chen [39] used a variation of Johnson and Kotz’s transformation for monitoring p, given by
y i = 2 n i sin 1 x i + 3 8 n + 3 4 sin 1 ( p ) , for i = 1 , 2 , ,
for the case where the parameter p is known. If the parameter p is unknown, then we use the transformation
y i = 2 n i sin 1 x i + 3 8 n + 3 4 sin 1 ( p ^ i 1 ) , for i = 2 , 3 , ,
where p ^ i = ( x 1 + x 2 + + x i ) / ( n 1 + n 2 + + n i ) . The arcsine-p chart, proposed by Chen [39], is obtained by plotting the statistics y i ’s with control limits at UCL = 3, CL = 0 and LCL = −3. We notice that for the case where p is unknown, the first plotted statistic is y 2 .

2.2.3. The Q Transformation

Quesenberry [40] introduced the binomial Q-chart for data with known and unknown p. For p = p 0 known a priori, we define
u i = k = 0 x i n i k p 0 k ( 1 p 0 ) ( n i k ) , Q i = Φ 1 ( u i ) , for i = 1 , 2 ,
According to Quesenberry [40], the values Q 1 , Q 2 , , are approximately standard normal.
For p unknown, let t i = j = 1 i x j and N i = j = 1 i n j and define
u i = k = h i x i n i k t i k N i 1 n i + N i 1 t i Q i = Φ 1 ( u i ) , for i = 1 , 2 ,
where h i = max ( 0 , t i N i 1 ) . Then, the values Q 2 , Q 3 , are approximately standard normal. The binomial Q-chart is obtained by plotting the statistics Q on a chart with control limits UCL = 3, CL = 0 and LCL = −3. We notice that for the case where p is unknown, the first plotting statistic is Q 2 . Furthermore, the approximation has adequate performance after the first observations, as in the case of Poisson data.

2.3. Transformation Techniques for Negative Binomial Data

Let X be the total number of products inspected until a nonconforming r ( r 1 ) is counted and p be the probability of observing a nonconforming one. Then, the variable X follows the negative binomial distribution with parameters r and p with a pmf given by
f ( x ; r , p ) = x 1 r 1 p r ( 1 p ) x r , x = r , r + 1 ,
The mean and variance of X are equal to r / p and r ( 1 p ) / p 2 , respectively. Note that the geometric distribution is a special case of the negative binomial distribution for r = 1 . The CCC-r chart, proposed by Xie et al. [41], is used to monitor observations of the variable X and it is suggested when the parameter p is low. For more details, the reader can refer to the work of Xie, Goh and Kuralmani [42] and that of Joekes, Smrekar and Righetti [43].

2.3.1. Anscombe’s Transformation

Let x i be the total number of products inspected until a nonconforming r ( r 1 ) is counted. Anscombe [31] showed that the transformation
y i = ln ( x i + 1 2 r ) , for i = 1 , 2 ,
for a large mean m = r / p and r 1 is approximately normal with variance approximated by 1 / ( r 0.5 ) for a large r.

2.3.2. The Box–Cox Transformation

Box and Cox [44] proposed a family of power transformations on a positive variable X given by
Y ( L ) = X L 1 L , for L 0 , ln ( X ) , for L = 0 ,
where L is a single parameter that can be estimated using the maximum likelihood estimation (MLE), so that Y ( L ) is approximately normally distributed. Applying this transformation to observations x i , we have the transformed data y i , which follow approximately the normal distribution.

3. A Brief Overview of PCIs for Discrete Distributions

Yeh and Bhattcharga [21] proposed an index based on the proportion of nonconforming products. This index is defined as follows:
C f = m i n p 0 L p L , p 0 U p U ,
where p L = P ( X < L ) , p U = P ( X > U ) and p 0 L , p 0 U are the expected proportions of nonconforming products that the manufacturer can tolerate on L and U, respectively.
Perakis and Xekalaki [23] proposed an index that can be used for continuous as well as for discrete distributions and is defined as follows:
C p c = 1 p 0 1 p ,
where p = P ( L < X < U ) and p 0 is the minimum allowable proportion of conformance, usually equal to 0.9973. Perakis and Xekalaki [24] investigated the above index for discrete data where only one specification limit is set. More specifically, if only the U is set, then the index is denoted by C p c u while if only the L is set, the index is denoted by C p c l . The two indices are defined as follows:
C p c u = 0.0027 1 P ( X < U ) , C p c l = 0.0027 1 P ( X > L ) .
Alternative forms of the two indices for a Poisson process can be obtained by using the chi-square distribution. Thus, they can be written as follows:
C p c u = 0.0027 P ( x 2 U 2 < 2 c ) , C p c l = 0.0027 P ( x 2 ( L + 1 ) 2 > 2 c ) .
Maiti, Saha and Nanda [25] proposed an index based on the process yield, which is defined as follows:
C p y = p p 0 ,
where p = P ( L X U ) and p 0 is the desired yield equal to 1 p 0 L p 0 U .
For an off-centered process, i.e., F ( U ) + F ( L ) 1 , where F ( · ) is the cumulative density function (cdf), the index is defined as follows:
C p y k = m i n F ( U ) F ( M ) 1 2 p 0 U , F ( M ) F ( L ) 1 2 p 0 L = m i n F ( U ) 1 2 1 2 p 0 U , 1 2 F ( L ) 1 2 p 0 L ,
where M is the median of the distribution and F ( M ) = F ( U ) + F ( L ) 2 = 1 2 for a centered process. We notice that for a centered process, C p y = C p y k , while if C p y k < C p y , the process median is off-centered.
For generalized asymmetric tolerance, i.e., T M , the index is defined as follows:
C p T k = m i n F ( U ) F ( T ) 1 2 p 0 U , F ( T ) F ( L ) 1 2 p 0 L .
We point out that C p T k = C p y k when M = T .
The above PCIs have been studied for various non-normal as well as for Poisson distributions. However, these indices can be easily extended to other discrete data following a binomial or negative binomial distribution using the appropriate cdf.
When the value of the parameter of distribution is unknown, then it has to be estimated in order to compute the indices. In most cases, the maximum likelihood estimator (MLE) method is used. If X 1 , X 2 , , X m are independent observations from a Poison process, the MLE of c is given by c ^ = 1 m i = 1 m X i , while if the observations are from a binomial process, the MLE of p is given by p ^ = 1 m n i = 1 m X i . Finally, for a negative binomial process, the MLE of p is given by p ^ = m r i = 1 m X i .

4. A Simulation Study

In this section, a simulation study is performed to evaluate the PCIs and PPIs for different Poisson, binomial and negative binomial processes. Six different processes are considered for each distribution: two centered, two off-centered and two off-target processes. For the centered processes, the process target is approximately equal to the median, and the U and L are chosen so that the condition F ( U ) + F ( L ) = 1 is satisfied. For the off-centered processes, the values of specification limits are chosen so that F ( U ) + F ( L ) 1 . Finally, under off-target processes, the process target is such that F ( T ) F ( U ) + F ( L ) / 2 . These situations of a process are also considered in Maiti, Saha and Nanda [25].
In our simulation study, random variables of each discrete distribution are generated for specified values of parameters and for sample of sizes m = 50 and 100 to compute the mean and the standard deviation of each index. In order to estimate them, first, the appropriate transformation method is applied to the data. Then, the specification limits and the target value are transformed via the transformation applied on the data. Next, the obtained specification limits, the target value, as well as the mean and the standard deviation of the transformed data are used to estimate each index. As C p has an indirect association with C p c and C p y , C p k with C p y k and C p m k with C p T k , comparisons among them and the classical PCIs are also presented. The estimators of the C p c , C p y , C p y k and C p T k are obtained using the MLE method. Moreover, in order to make the latter indices comparable with the classical PCIs, we use p 0 = 0.9973 and p 0 L = p 0 U = 0.00135 . All the computations are carried out using the R version 4.3.0. software with 10,000 iterations.
The results of the simulation study are presented in Table 1, Table 2, Table 3, Table 4, Table 5 and Table 6. The first two processes in each table are centered, the other two are off-centered and the last two are off-target. The first columns in these tables show the values of the parameters of each distribution, the specification limits and the target value of the processes. In each process, the first row shows the mean of the index estimators and the second row shows their standard deviations.
As can be observed in the aforementioned tables, we observe that for each process, the differences in the mean values of PCI and PPI estimators applying different transformation methods are minor. More specifically, the differences are very small when applying Anscombe’s and Freeman–Tukey’s transformations for Poisson processes and Freeman–Tukey’s and Chen’s transformations for binomial processes. Moreover, for all processes, the values of C ^ p are very close to the corresponding values of C ^ p y , while at almost all of them, the values of C ^ p c differ from those of C ^ p and C ^ p y . This can be explained as follows: The C p c index is based on the ratio of the allowable proportion of nonconforming products to the observed proportion of nonconforming products, while the C p y index is based on the ratio of the observed proportion of conformance to the allowable proportion of conformance. In other words, the C p y index is the reciprocal of the complementary proportions of the C p c index. For off-centered processes, the values of C ^ p k indices are slightly smaller than those of C ^ p y k . Furthermore, for binomial and negative binomial off-target processes, the values of C ^ p m k indices are very close to the corresponding C ^ p T k indices, while for Poisson processes, the differences are larger but both of them indicate incapability processes. Finally, it is observed that the values of P ^ p and P ^ p k are approximately similar to those of C ^ p and C ^ p k , respectively. Thus, the processes can be considered stable.
In addition, we study the impact of sample size and parameters of distribution on the results. For a certain process, the differences in the mean values of each index estimator are negligible for different values of sample sizes. This entails the fact that one should use a small value of sample size (say m = 50 ) to obtain safe conclusions about the process capability/performance. On the other hand, the value of the standard deviation of each index estimator for m = 50 is always higher than the corresponding value for m = 100 . It should be noticed that the standard deviations of C ^ p y , C ^ p y k and C ^ p T k are smaller than that of the classical PCIs estimators. Moreover, for certain values of U, L and T, the capability of the process depends on the parameters of distribution. For example, consider the second and the fifth processes in Table 1, where ( U , L , T ) = ( 15 , 0 , 6 ) , while c = 6.5 for the second process and c = 4.5 for the fifth one. It is recalled that the second process is centered and the fifth process is off-target. As can be seen in Table 1, we observe that the differences in the mean and standard deviation values of the C ^ p and P ^ p indices for the two processes are negligible, while the corresponding differences are significant for the other three classical PCIs and P ^ p k . This was expected as the C p and P p indices do not depend on the mean value of the process.

5. Applications

To illustrate the application of the above transformations, three examples are presented. The values of the PCIs and PPIs are computed with particular emphasis on describing Anscombe’s transformation for Poisson and negative binomial data and the Freeman–Tukey’s transformation for binomial data. Similarly, the indices can be calculated using the other transformations techniques. The C p y , C p y k and C p T k indices are also computed for comparison reasons. In the following examples, we present the tools and we guide practitioners on how to apply the proposed methodology. The control charts are provided in order to confirm that the process is in-control. In addition, we apply the transformation techniques to compute the PCIs and PPIs, and we discuss the results afterwards.
Example 1. 
We use the data on the number of nonconformities in samples of 100 printed circuit boards, which is described in more detail in Chapter 7 by Montgomery [12]. The dataset consists of 46 samples. The first 26 samples of Table 7 are from Table 7.7 of Montgomery [12] and the last 20 samples are from the Table 7.8. of Montgomery [12] The data are presented in Table 7. Moreover, it is assumed that the target value is T = 18 nonconformities in each sample, and the upper and lower specification limits are U = 35 and L = 5 nonconformities in each sample, respectively.
After a retrospective data analysis, two values (in bold) in Table 7 were eliminated due to assignable causes. Fitting in the remaining 44 samples the Poisson distribution, we find that c ^ = 19.0455 is the MLE (p-value = 0.4669). In Figure 1, the c-chart is presented. Note that the control limits of the c-chart are computed with
U C L = c ^ + 3 c ^ = 19.0455 + 3 19.0455 = 32.1378 ,
C L = c ^ = 19.0455 ,
L C L = c ^ 3 c ^ = 19.0455 3 19.0455 = 5.9532 .
As no point and no systematic behaviour can be seen from the c-chart, we can conclude that the process is in control. In addition, we compute all the y i ’s values, for i = 1 , 2 , , 44 using Anscombe’s transformation given by (2). The mean of the transformed data is μ ^ = 8.7371 . In order to estimate the standard deviation of the transformed data, we compute their moving range. Hence, the standard deviation σ ^ is given by
σ ^ = M R ¯ d 2 = 1.2517 1.1284 = 1.1093
An I M R chart is presented in Figure 2.
The control limits for the I-chart are computed with
U C L = μ ^ + 3 σ ^ = 8.7371 + 3 · 1.1093 = 12.065 ,
C L = μ ^ = 8.7371 ,
L C L = μ ^ 3 σ ^ = 8.7371 3 · 1.1093 = 5.4092 ,
and those for the M R -chart are computed with
U C L = D 4 · M R ¯ = 3.267 · 1.2517 = 4.0893 ,
C L = M R ¯ = 1.2517 ,
L C L = D 3 · M R ¯ = 0 · 1.2517 = 0 .
From the I M R chart in Figure 2, we confirm that the process is in control. The last step before computing the four PCIs is to transform the values of T, U and L, using transformation (2), and as a result, we have T = 8.5732 , U = 12.2270 and L = 4.1833 . The values of the four PCIs are computed with
C p = U L 6 σ ^ = 12.2270 4.1833 6 · 1.1093 = 1.2085 ,
C p k = m i n { U μ ^ , μ ^ L } 3 σ ^ = m i n { 12.2270 8.7371 , 8.7371 4.1833 } 3 · 1.1093 = 1.0487 ,
C p m = U L 6 σ ^ 2 + ( μ ^ T ) 2 = 12.2270 4.1833 6 1 . 1093 2 + ( 8.7371 8.5732 ) 2 = 1.1955 ,
C p m k = m i n { U μ ^ , μ ^ L } 3 σ ^ 2 + ( μ ^ T ) 2 = m i n { 12.2270 8.7371 , 8.7371 4.1833 } 3 1 . 1093 2 + ( 8.7371 8.5732 ) 2 = 1.0374 .
The long-term standard deviation of the transformed data is computed with
σ s ^ = i = 1 44 ( y i μ ^ ) 2 n 1 ,
where y i are the transformed data. Thus, σ s ^ = 1.1732 . The values of the two PPIs are equal to
P p = U L 6 σ s ^ = 12.2270 4.1833 6 · 1.1732 = 1.1427 ,
P p k = m i n { U μ ^ , μ ^ L } 3 σ s = m i n { 12.2270 8.7371 , 8.7371 4.1833 } 3 · 1.1732 = 0.9916 .
Similarly, we compute the values of the indices with the other transformation techniques as well as the values of the C p y , C p y k and C p T k indices. The results are presented in Table 8. From this table, we conclude that the differences in the values of the indices applying different transformation techniques are very small. Moreover, the values of C p , C p k and C p m k are very close to the corresponding values of C p y , C p y k and C p T k , while the values of P p and P p k are approximately similar to those of C p and C p k . Thus, one should consider that the process is stable. As the PCIs are higher than one, except for C p T k , one should conclude that the process is capable. Consequently, the transformation techniques give similar results and one should apply them in order to make conclusions about the capability of the process.
Example 2. 
We use the data on a number of nonconforming cans with frozen orange juice in samples of size n = 50 cans from an industrial process, which is also described in detail, as the previous example, in Chapter 7 by Montgomery [12]. The dataset consists of 40 samples and is presented in Table 9. It is also assumed that the target value is T = 5 defective bellows in each sample and that the upper and lower specification limits are U = 13 and L = 1 defective bellows per sample, respectively.
Applying the binomial distribution to these data, the MLE of the parameter p is 0.109 (p-value = 0.9999).
In Figure 3, the p-chart is presented. Note that the control limits of the p-chart are computed with
U C L = p ^ + 3 p ^ ( 1 p ^ ) n = 0.109 + 3 0.109 ( 1 0.109 ) 50 = 0.2412 ,
C L = p ^ = 0.109 ,
L C L = p ^ 3 p ^ ( 1 p ^ ) n = 0.109 3 0.109 ( 1 0.109 ) 50 = 0.0232 L C L = 0 .
From the p-chart in Figure 3, we conclude that the process is in control as no point and no systematic behaviour can be seen. In addition, we compute all the y i ’s values, for i = 1 , 2 , , 40 using Freeman–Tukey’s transformation, given in (7). The mean and standard deviation of the transformed data are μ ^ = 0.6846 and σ ^ = 0.1176 , respectively. An I M R chart for the transformed data is presented in Figure 4.
The control limits for the I-chart are calculated as
U C L = μ ^ + 3 σ ^ = 0.6846 + 3 · 0.1176 = 1.0374 ,
C L = μ ^ = 0.6846 ,
L C L = μ ^ 3 σ ^ = 0.6846 3 · 0.1176 = 0.3318 ,
and those for the M R -chart are calculated with
U C L = D 4 · M R ¯ = 3.267 · 0.1327 = 0.4335 ,
C L = M R ¯ = 0.1327 ,
L C L = D 3 · M R ¯ = 0 · 0.1327 = 0 .
From the I M R chart in Figure 4, we confirm that the process is in control. The transformed values of T, U and L are T = 0.6686 , U = 1.0807 and L = 0.3398 , respectively, and the values of the four PCIs are computed with
C p = U L 6 σ ^ = 1.0807 0.3398 6 · 0.1176 = 1.0500 ,
C p k = m i n { U μ ^ , μ ^ L } 3 σ ^ = m i n { 1.0807 0.6846 , 0.6846 0.3398 } 3 · 0.1176 = 0.9773 ,
C p m = U L 6 σ ^ 2 + ( μ ^ T ) 2 = 1.0807 0.3398 6 0 . 1176 2 + ( 0.6846 0.6686 ) 2 = 1.0404 ,
C p m k = m i n { U μ ^ , μ ^ L } 3 σ ^ 2 + ( μ ^ T ) 2 = m i n { 1.0807 0.6846 , 0.6846 0.3398 } 3 0 . 1176 2 + ( 0.6846 0.6686 ) 2 = 0.9684 .
The long-term standard deviation of the transformed data is equal to σ s ^ = 0.1415 . Thus, the values of the PPIs are
P p = U L 6 σ s ^ = 1.0807 0.3398 6 · 0.1415 = 0.8725 ,
P p k = m i n { U μ ^ , μ ^ L } 3 σ s = m i n { 1.0807 0.6846 , 0.6846 0.3398 } 3 · 0.1415 = 0.8121 .
The values of the indices with the other transformations techniques as well as the values of the C p y , C p y k and C p T k indices are presented in Table 10. We conclude that the differences in the values of the PCIs applying different transformation techniques are very small. The values of C p , C p k and C p m k are close to the values of C p y , C p y k and C p T k , respectively. As C p k and C p m k and of course C p y k and C p T k are lower than one, one should conclude that the process is incapable. Moreover, the differences between C p and P p as well as C p k and P p k are significant. Thus, the process is unstable. Similarly to the previous example, the transformation techniques give similar results and one should apply them in order to make conclusions about the capability/performance of the process. However, one should be careful, as C p and C p m are larger than one while C p k , C p m k , P p and P p k are smaller than one; thus, one should evaluate all the indices to make safe conclusions about the capability/performance of the process.
Example 3. 
We assume a process in which we count the total number of inspected products until r = 5 nonconforming items are found. It is also known that the probability p of observing a nonconforming item is equal to 0.1; the target value is T = 45 inspected products; and the upper and lower specification limits are U = 140 and L = 5 inspected products, respectively. The dataset can be found in Chapter 3 of Xie, Goh, and Kuralmani [42], and it is presented row by row in Table 11.
The control limits of a CCC-r chart are given as the solution of the above equations
F ( U C L r , r , p ) = i = r U C L r i 1 r 1 p r ( 1 p ) i r = 1 a 2 ,
F ( C L r , r , p ) = i = r C L r i 1 r 1 p r ( 1 p ) i r = 0.5 ,
F ( L C L r , r , p ) = i = r L C L r i 1 r 1 p r ( 1 p ) i r = a 2 ,
where α is the acceptable risk of a false alarm. In this example, the control limits of the CCC-5 chart under a standard false alarm probability level α = 0.0027 are U C L 5 = 134 , C L 5 = 41 and L C L 5 = 10 .
In Figure 5, the CCC-5 chart is presented from which we conclude that the process is in control as no point and no systematic behaviour can be seen. In addition, we compute the y i values, for i = 1 , 2 , , 100 using Anscombe’s transformation, given in (14). The mean and standard deviation of the transformed data are μ ^ = 3.7978 and σ ^ = 0.4885 , respectively. An I M R chart for the transformed data is presented in Figure 6.
The control limits for the I-chart are computed with
U C L = μ ^ + 3 σ ^ = 3.7978 + 3 · 0.4885 = 5.2633 ,
C L = μ ^ = 3.7978 ,
L C L = μ ^ 3 σ ^ = 3.7978 3 · 0.4885 = 2.3323 ,
and those for the M R -chart are computed with
U C L = D 4 · M R ¯ = 3.267 · 0.5512 = 1.8008 ,
C L = M R ¯ = 0.5512 ,
L C L = D 3 · M R ¯ = 0 · 0.5512 = 0 .
From the I M R chart in Figure 6, we note that one point, that of sample 39, is plotted over the U C L in the M R -chart. However, no assignable cause is determined, and therefore, we decide that the process is in-control. The transformed values of T, U and L are T = 3.8607 , U = 4.9593 and L = 2.0149 , respectively, and the values of the four PCIs are computed with
C p = U L 6 σ ^ = 4.9593 2.0149 6 · 0.4885 = 1.0046 ,
C p k = m i n { U μ ^ , μ ^ L } 3 σ ^ = m i n { 4.9593 3.7978 , 3.7978 2.0149 } 3 · 0.4885 = 0.7926 ,
C p m = U L 6 σ ^ 2 + ( μ ^ T ) 2 = 4.9593 2.0149 6 0 . 4885 2 + ( 3.7978 3.8607 ) 2 = 0.9963 ,
C p m k = m i n { U μ ^ , μ ^ L } 3 σ ^ 2 + ( μ ^ T ) 2 = m i n { 4.9593 3.7978 , 3.7978 2.0149 } 3 0 . 4885 2 + ( 3.7978 3.8607 ) 2 = 0.7861 .
The long-term standard deviation of the transformed data is equal to σ s ^ = 0.4561 . Thus, the two PPIs are computed with
P p = U L 6 σ s ^ = 4.9593 2.0149 6 · 0.4561 = 1.0760 ,
P p k = m i n { U μ ^ , μ ^ L } 3 σ s = m i n { 4.9593 3.8607 , 3.8607 2.0149 } 3 · 0.4561 = 0.8489 .
The values of the indices applying the two transformation methods are presented in Table 12. Moreover, the values of C p y , C p y k and C p T k are also presented in the same table. We note that the value of parameter L using the Box–Cox transformation is equal to −0.5415. Minor differences are observed in the C p and C p m values applying the two transformation methods. However, the values of the indices indicate that the process is incapable and unstable. In this example, there are small differences in the computation of the four PCIs and two PPIs. However, similar to the previous example, one should evaluate all of the indices as only C p and P p are larger than the one for Anscombe’s transformation, while C p , C p m and P p are larger than the one for the Q transformation.
As a general conclusion from the above three applications, one practitioner should compute all four PCIs and two PPIs in order to make safe conclusions about the capability and the performance of the process. It should also be pointed out that if a practitioner evaluates only C p or C p m , then they may wrongly conclude that the process is capable. Similarly, one should evaluate not only the P p but also the P p k index to conclude about the behaviour of the process. The results from different transformation methods are similar. Therefore, one should compute the PCIs and PPIs by applying any method. We also note that the conclusions about the capability and the performance of the processes obtained by applying the transformation techniques are the same as those obtained by computing the indices for non-normal distributions; thus, practitioners do not have to know the indices for non-normal data.

6. Conclusions

Process capability and performance indices, named PCIs and PPIs, respectively, are practical and useful tools for continuous quality improvement. There are various methods in the literature to estimate PCIs for non-normal and continuous processes. In the case of discrete processes, new PCIs have been defined by researchers, with their use causing difficulties to many practitioners. In the present paper, we proposed various transformation techniques for data following Poisson, binomial or negative binomial distributions. The transformed data are then used to compute the well-known PCIs and PPIs. Through a simulation study, we demonstrated that the differences in the indices applying different transformation methods are very small. Comparisons with other PCIs for discrete data showed that the proposed methodology can safely be applied without any misinterpretation of the process capability. Thus, we suggest that practitioners who are not familiar with indices for non-normal data use the transformation methodology in order to obtain safe conclusions about the process capability and performance. Finally, we applied the transformation techniques in three illustrative examples for data, following different distribution.
For future research, we recommend the extension of these transformation techniques to multivariate discrete processes. Moreover, one should propose indices for data following the Conway–Maxwell Poisson distribution, which encompasses the geometric, Poisson and Bernoulli distributions as special cases, or zero-inflated distributions.

Funding

This research received no external funding.

Data Availability Statement

The data that support the findings of this study are available from the respective references as mentioned in the main text.

Acknowledgments

The author thank the editors and the referees for their useful comments, which resulted in improving the quality of this article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Juran, J.M. Juran’s Quality Control Handbook, 3rd ed.; McGraw-Hill: New York, NY, USA, 1974. [Google Scholar]
  2. Kane, V.E. Process capability indices. J. Qual. Technol. 1986, 18, 41–52. [Google Scholar] [CrossRef]
  3. Chan, L.K.; Cheng, S.W.; Spiring, F.A. A new measure of process capability: Cpm. J. Qual. Technol. 1988, 20, 162–175. [Google Scholar] [CrossRef]
  4. Hsiang, T.C.; Taguchi, G. A tutorial on quality control and assurance—The Taguchi methods. In Proceedings of the ASA Annual Meeting, Las Vegas, NV, USA, 5–8 August 1985; p. 188. [Google Scholar]
  5. Choi, B.C.; Owen, D.B. A study of a new process capability index. Commun. Stat.-Theory Methods 1990, 19, 1232–1245. [Google Scholar] [CrossRef]
  6. Pearn, W.L.; Kotz, S.; Johnson, N.L. Distributional and inferential properties of process capability indices. J. Qual. Technol. 1992, 24, 216–231. [Google Scholar] [CrossRef]
  7. Rodriguez, R.N. Recent developments in process capability analysis. J. Qual. Technol. 1992, 24, 176–187. [Google Scholar] [CrossRef]
  8. Kotz, S.; Johnson, N.L. Process Capability Indices; Chapman and Hall: London, UK, 1993. [Google Scholar]
  9. Kotz, S.; Johnson, N.L. Process capability indices—A review, 1992–2000. J. Qual. Technol. 2002, 34, 2–19. [Google Scholar] [CrossRef]
  10. Kotz, S.; Lovelace, C.R. Process Capability Indices in Theory and Practice; Arnold: London, UK, 1998. [Google Scholar]
  11. Pearn, W.L.; Kotz, S. Encyclopedia and Handbook of Process Capability Indices—A Comprehensive Exposition of Quality Control Measures; Series on Quality, Reliability and Engineering Statistics; World Scientific Publishing Co.: Singapore, 2006. [Google Scholar]
  12. Montgomery, D.C. Introduction to Statistical Quality Control, 7th ed.; John Wiley & Sons: New York, NY, USA, 2012. [Google Scholar]
  13. Clements, J.A. Process capability calculations for non-normal distributions. Qual. Prog. 1989, 22, 95–100. [Google Scholar]
  14. Pearn, W.L.; Kotz, S. Application of Clements’ method for calculating second and third-generation process capability indices for non-normal Pearsonian populations. Qual. Eng. 1994, 7, 139–145. [Google Scholar] [CrossRef]
  15. Rivera, L.A.R.; Hubele, N.F.; Lawrence, F.P. Cpk index estamation using data transformation. Comput. Ind. Eng. 1995, 29, 55–58. [Google Scholar] [CrossRef]
  16. Pearn, W.L.; Chen, K.S. Estimating process capability indices for non-normal Pearsonian populations. Qual. Reliab. Eng. Int. 1995, 11, 386–388. [Google Scholar] [CrossRef]
  17. Castagliola, P. Evaluation of non-normal process capability indices using Burr’s distributions. Qual. Eng. 1996, 8, 587–593. [Google Scholar] [CrossRef]
  18. Hosseinifard, S.Z.; Abbasi, B.; Ahmad, S.; Abdollahian, M. A transformation technique to estimate the process capability index for non-normal processes. Int. J. Adv. Manuf. Technol. 2009, 40, 512–517. [Google Scholar] [CrossRef]
  19. Hosseinifard, S.Z.; Abbasi, B.; Niaki, S.T.A. Process capability estimation for leukocyte filtering process in blood service: A comparison study. IIE Trans. Healthc. Syst. Eng. 2014, 4, 167–177. [Google Scholar] [CrossRef]
  20. Bothe, D.R. Measuring Process Capability: Techniques and Calculations for Quality and Manufacturing Engineers; McGraw-Hill: New York, NY, USA, 1997. [Google Scholar]
  21. Yeh, A.B.; Bhattacharya, S. A robust process capability index. Commun. Stat. Simul. Comput. 1998, 27, 565–589. [Google Scholar] [CrossRef]
  22. Borges, W.; Ho, L.L. A fraction defective based capability index. Qual. Reliab. Eng. Int. 2001, 17, 447–458. [Google Scholar] [CrossRef]
  23. Perakis, M.; Xekalaki, E. A process capability index that is based on the proportion of conformance. J. Stat. Comput. Simul. 2002, 72, 707–718. [Google Scholar] [CrossRef]
  24. Perakis, M.; Xekalaki, E. A process capability index for discrete processes. J. Stat. Comput. Simul. 2005, 75, 175–187. [Google Scholar] [CrossRef]
  25. Maiti, S.S.; Saha, M.; Nanda, A.K. On generalizing process capability indices. Qual. Technol. Quant. Manag. 2010, 7, 279–300. [Google Scholar] [CrossRef]
  26. Maravelakis, P.E. Process capability indices for data following the Poisson or binomial distribution. Qual. Technol. Quant. Manag. 2016, 13, 197–206. [Google Scholar] [CrossRef]
  27. Dey, S.; Saha, M. Bootstrap confidence intervals of generalized process capability index Cpyk using different methods of estimation. J. Appl. Stat. 2019, 46, 1843–1869. [Google Scholar] [CrossRef]
  28. Pal, S.; Gauri, S.K. Measuring capability of a binomial process. Int. J. Eng. Sci. Technol. 2020, 12, 25–37. [Google Scholar] [CrossRef]
  29. Gauri, S.K.; Pal, S. A Note on the Generalized Indices of Process Capability. Int. J. Res. Ind. Eng. 2020, 9, 286–303. [Google Scholar]
  30. Bartlett, M.S. The square root transformation in analysis of variance. Suppl. J. R. Stat. Soc. 1936, 3, 68–78. [Google Scholar] [CrossRef]
  31. Anscombe, F.J. The transformation of Poisson, binomial and negative-binomial data. Biometrika 1948, 35, 246–254. [Google Scholar] [CrossRef]
  32. Freeman, M.F.; Tukey, J.W. Transformations related to the angular and the square root. Ann. Math. Stat. 1950, 21, 607–611. [Google Scholar] [CrossRef]
  33. Ryan, T.P.; Schwertman, N.C. Optimal limits for attributes control charts. J. Qual. Technol. 1938, 29, 86–98. [Google Scholar] [CrossRef]
  34. Quesenberry, C.P. SPC Q charts for start-up processes and short and long runs. J. Qual. Technol. 1991, 23, 213–224. [Google Scholar] [CrossRef]
  35. Quesenberry, C.P. SPC Q charts for a Poisson parameter λ: Short or long runs. J. Qual. Technol. 1991, 23, 296–303. [Google Scholar] [CrossRef]
  36. Quesenberry, C.P. On properties of Poisson Q charts for attributes. J. Qual. 1995, 27, 293–303. [Google Scholar]
  37. Ryan, T.P. Statistical Methods for Quality Improvement; John Wiley & Sons: New York, NY, USA, 1989. [Google Scholar]
  38. Johnson, N.L.; Kotz, S. Discrete Distributions (Vol. 1 of the 4-Volume Set: Distributions in Statistics); Wiley: New York, NY, USA, 1969. [Google Scholar]
  39. Chen, G. An improved p chart through simple adjustment. J. Qual. Technol. 1998, 30, 142–151. [Google Scholar] [CrossRef]
  40. Quesenberry, C.P. SPC Q charts for a binomial parameter p: Short or long runs. J. Qual. Technol. 1991, 23, 239–246. [Google Scholar] [CrossRef]
  41. Xie, M.; Lu, X.S.; Goh, T.N.; Chan, L.Y. A quality monitoring and decision-making scheme for automated production processes. Int. J. Qual. Reliab. Manag. 1999, 16, 148–157. [Google Scholar] [CrossRef]
  42. Xie, M.; Goh, T.N.; Kuralmani, V. Statistical Models and Control Charts for High Quality Processes; Kluwer Academic Publishers: Boston, MA, USA, 2002. [Google Scholar]
  43. Joekes, S.; Smrekar, M.; Righetti, A.F. A comparative study of two proposed CCC-r charts for high quality processes and their application to an injection molding process. Qual. Eng. 2016, 28, 467–475. [Google Scholar] [CrossRef]
  44. Box, G.E.P.; Cox, D.R. An analysis of transformations. J. R. Stat. Society. Ser. B Methodol. 1964, 26, 211–252. [Google Scholar] [CrossRef]
Figure 1. The c-chart for the Example 1 data.
Figure 1. The c-chart for the Example 1 data.
Mathematics 11 03457 g001
Figure 2. The I M R chart for the Example 1 data.
Figure 2. The I M R chart for the Example 1 data.
Mathematics 11 03457 g002
Figure 3. The p-chart for the Example 2 data.
Figure 3. The p-chart for the Example 2 data.
Mathematics 11 03457 g003
Figure 4. The I M R chart for the Example 2 data.
Figure 4. The I M R chart for the Example 2 data.
Mathematics 11 03457 g004
Figure 5. The CCC-5 chart for the Example 3 data.
Figure 5. The CCC-5 chart for the Example 3 data.
Mathematics 11 03457 g005
Figure 6. The I M R chart for the Example 3 data.
Figure 6. The I M R chart for the Example 3 data.
Mathematics 11 03457 g006
Table 1. The results of the simulation study for Poisson processes (sample size m = 50 ).
Table 1. The results of the simulation study for Poisson processes (sample size m = 50 ).
C ^ p C ^ pc C ^ py C ^ pk C ^ pyk C ^ pm
AnscF-TQAnscF-TQAnscF-TQ
(U, L) = (10, 0)0.89940.94770.84250.14390.99540.73640.73450.78080.97730.88230.92990.8220
(c, T) = (4.5, 4)0.12380.13280.11320.01010.00350.10680.10720.11330.00610.11680.12520.1081
(U, L) = (15, 0)1.13501.18681.04930.93651.00140.92290.92410.98710.99871.11641.16731.0284
(c, T) = (6.5, 6)0.15400.16150.14000.09840.00080.12930.12980.13800.00120.14660.15370.1343
(U, L) = (20, 0)1.34521.39321.38680.25601.00271.06241.16090.89490.97951.31961.36701.3531
(c, T) = (4.5, 4)0.18520.19530.18640.0798 4 · 10 8 0.16030.17720.12890.00710.17470.18400.1779
(U, L) = (12, 0)0.99671.04810.89180.16010.98560.64620.64670.67830.96840.98031.03080.8741
(c, T) = (6.5, 6)0.13520.14260.11900.05070.00690.09560.09590.10330.01380.12870.13580.1141
(U, L) = (15, 0)1.14081.18881.13090.25451.00271.06021.14370.89490.97950.90010.93880.8952
(c, T) = (4.5, 6)0.15710.16660.15200.0779 2 · 10 5 0.15820.16730.12890.00710.10720.11420.0981
(U, L) = (13, 0)1.04451.09600.93680.31490.99500.74180.74260.77860.98710.73590.77230.6619
(c, T) = (6.5, 9)0.14170.14910.08630.08470.00350.10700.10730.07970.00700.07790.08190.0454
C ^ pmk C ^ pTk P ^ p P ^ pk
AnscF-TQAnscF-TQAnscF-TQ
(U, L) = (10, 0)0.72330.72160.76500.88570.88550.92980.83350.72500.72060.7748
(c, T) = (4.5, 4)0.10750.10760.11530.07450.09640.10690.08550.08550.08730.0899
(U, L) = (15, 0)0.90850.90970.96830.90071.12211.17271.04080.91240.91310.9792
(c, T) = (6.5, 6)0.12880.12920.13910.07460.12140.12860.10760.10370.10480.1096
(U, L) = (20, 0)1.04131.13810.87200.89341.32431.36691.37211.04591.13900.8854
(c, T) = (4.5, 4)0.14500.16090.11440.07540.14420.15710.14070.13070.14820.1014
(U, L) = (12, 0)0.63640.63680.66570.88270.98541.03560.88470.63880.63900.6729
(c, T) = (6.5, 6)0.09650.09670.10540.07260.10660.11360.09150.07870.07930.0853
(U, L) = (15, 0)0.83870.90420.71110.34381.12311.16641.11881.04381.12220.8854
(c, T) = (4.5, 6)0.12620.12270.10670.07720.12230.13410.11480.12850.13770.1014
(U, L) = (13, 0)0.52040.52100.54890.23551.03271.08300.93810.73340.73370.7796
(c, T) = (6.5, 9)0.03590.03600.02570.05570.11180.11880.09700.08700.08780.0945
Table 2. The results of the simulation study for Poisson processes (sample size m = 100 ).
Table 2. The results of the simulation study for Poisson processes (sample size m = 100 ).
C ^ p C ^ pc C ^ py C ^ pk C ^ pyk C ^ pm
AnscF-TQAnscF-TQAnscF-TQ
(U, L) = (10, 0)0.89030.93780.83440.14770.99570.72960.72750.78070.97910.87880.92590.8186
(c, T) = (4.5, 4)0.08540.09170.07800.00580.00230.07380.07400.07980.00420.08170.08760.0754
(U, L) = (15, 0)1.12391.17511.03940.97131.00150.91430.91550.98240.99921.11201.16261.0244
(c, T) = (6.5, 6)0.10550.11080.09580.06100.00050.08880.08910.09610.00070.10200.10700.0932
(U, L) = (20, 0)1.33161.37871.37350.24961.00271.05111.14810.88620.98001.31431.36121.3475
(c, T) = (4.5, 4)0.12770.13480.12840.0540 2 · 10 8 0.11060.12240.08880.00490.12220.12870.1241
(U, L) = (12, 0)0.98691.03770.88350.15710.98600.64030.64080.67200.96940.97651.02670.8707
(c, T) = (6.5, 6)0.09270.09780.08140.03650.00470.06570.06590.07100.00950.08950.09450.0792
(U, L) = (15, 0)1.12931.17641.12000.24891.00271.05081.14030.88620.98000.89570.93380.8920
(c, T) = (4.5, 6)0.10830.11500.10470.0532 10 5 0.11040.11780.08880.00490.07440.07950.0680
(U, L) = (13, 0)1.03431.08520.93680.31470.99530.73500.73570.77860.98780.73310.76930.6619
(c, T) = (6.5, 9)0.09710.10230.08630.06420.00240.07350.07370.07970.00480.05400.05690.0454
C ^ pmk C ^ pTk P ^ p P ^ pk
AnscF-TQAnscF-TQAnscF-TQ
(U, L) = (10, 0)0.72050.71870.76640.90640.87770.92070.82700.71920.71420.7737
(c, T) = (4.5, 4)0.07510.07510.08230.05700.06650.07410.05880.05920.06050.0639
(U, L) = (15, 0)0.90500.90610.96870.92321.11241.16211.03260.90490.90540.9760
(c, T) = (6.5, 6)0.08960.08990.09820.05660.08370.08900.07370.07160.07250.0768
(U, L) = (20, 0)1.03701.13310.86890.91561.31271.35361.36131.03611.12720.8783
(c, T) = (4.5, 4)0.10180.11310.08010.05750.09950.10900.09690.09040.10300.0699
(U, L) = (12, 0)0.63390.63430.66270.90180.97681.02630.87770.63370.63370.6677
(c, T) = (6.5, 6)0.06720.06740.07350.05570.07350.07860.06270.05440.05490.0589
(U, L) = (15, 0)0.83460.90580.70720.34171.11321.15501.11001.03591.11960.8783
(c, T) = (4.5, 6)0.08970.08940.07390.05450.08440.09300.07900.09010.09780.0699
(U, L) = (13, 0)0.51990.52040.54890.23381.02371.07320.93070.72750.72760.7736
(c, T) = (6.5, 9)0.02500.02510.02570.03950.07700.08220.06650.06010.06080.0653
Table 3. The results of the simulation study for binomial processes (sample size m = 50 ).
Table 3. The results of the simulation study for binomial processes (sample size m = 50 ).
C ^ p C ^ pc C ^ py C ^ pk C ^ pyk C ^ pm
F-TChenQF-TChenQF-TChenQ
(U, L, T) = (13, 0, 5)1.12391.07240.99830.54751.00090.89640.89560.94790.99571.09191.04200.9641
(n, p) = (113, 0.05)0.15700.14900.13590.04970.00100.12880.12830.13430.00200.14470.13740.1268
(U, L, T) = (22, 1, 10)1.27681.27791.21155.53061.00241.08781.08631.14171.00201.25701.25811.1902
(n, p) = (105, 0.1)0.17310.17330.16270.85690.00020.15150.15120.15890.00030.16600.16620.1566
(U, L, T) = (10, 0, 5)0.96520.91410.82040.09490.97510.57900.57890.60300.94760.93770.88810.7922
(n, p) = (113, 0.05)0.13480.12700.11170.02560.01000.08930.08910.09550.01990.12430.11710.1042
(U, L, T) = (17, 1, 10)1.06341.06490.97950.17820.98550.66090.66030.68130.96841.04681.04840.9623
(n, p) = (105, 0.1)0.14420.14440.13160.06930.00650.09830.09820.10310.01310.13820.13850.1266
(U, L, T) = (11, 0, 6)1.02010.96890.88130.18950.99080.68880.68850.72480.97870.98210.93260.8530
(n, p) = (113, 0.05)0.14250.13460.12000.04170.00500.10250.10220.11000.00970.13570.12800.1123
(U, L, T) = (18, 1, 11)1.10781.10921.02730.36500.99430.74970.74900.77690.98581.07061.07190.9962
(n, p) = (105, 0.1)0.15020.15040.13800.14830.00360.10900.10890.11460.00720.14290.14310.1300
C ^ pmk C ^ pTk P ^ p P ^ pk
F-TChenQF-TChenQF-TChenQ
(U, L, T) = (13, 0, 5)0.87220.87150.91650.90991.10871.05910.98950.88430.88450.9396
(n, p) = (113, 0.05)0.12830.12780.13310.06630.12300.11490.10220.10240.10090.1039
(U, L, T) = (22, 1, 10)1.07161.07021.12240.90821.26871.26961.20561.08081.07931.1361
(n, p) = (105, 0.1)0.15050.15030.15910.07000.13430.13440.12480.11970.11950.1248
(U, L, T) = (10, 0, 5)0.56380.56370.58370.88830.95220.90270.81310.57110.57170.5977
(n, p) = (113, 0.05)0.09130.09120.09860.06240.10560.09800.08400.07330.07270.0785
(U, L, T) = (17, 1, 10)0.65130.65080.67010.89071.05651.05800.97480.65660.65600.6780
(n, p) = (105, 0.1)0.09920.09910.10480.06790.11180.11200.10090.08050.08040.0848
(U, L, T) = (11, 0, 6)0.66160.66130.70020.65081.00640.95690.87350.67950.68000.7184
(n, p) = (113, 0.05)0.08650.08610.09330.09480.11160.10380.09020.08290.08200.0886
(U, L, T) = (18, 1, 11)0.72320.72240.75210.70191.10061.10201.02230.74490.74410.7731
(n, p) = (105, 0.1)0.09320.09300.09840.09850.11650.11670.10580.08820.08810.0931
Table 4. The results of the simulation study for binomial processes (sample size m = 100 ).
Table 4. The results of the simulation study for binomial processes (sample size m = 100 ).
C ^ p C ^ pc C ^ py C ^ pk C ^ pyk C ^ pm
F-TChenQF-TChenQF-TChenQ
(U, L, T) = (13, 0, 5)1.11151.06080.98790.56461.00100.88710.88650.94470.99621.08651.03690.9593
(n, p) = (113, 0.05)0.10800.10240.09350.03100.00060.08870.08840.09390.00140.10070.09570.0884
(U, L, T) = (22, 1, 10)1.26401.26511.19965.76551.00241.07741.07591.13361.00211.25121.25241.1849
(n, p) = (105, 0.1)0.11950.11970.11240.60050.00010.10490.10470.11100.00020.11600.11610.1095
(U, L, T) = (10, 0, 5)0.95450.90410.81190.09380.97580.57320.57310.59690.94880.93310.88380.7883
(n, p) = (113, 0.05)0.09270.08730.07690.01880.00700.06170.06160.06610.01400.08650.08150.0727
(U, L, T) = (17, 1, 10)1.05271.05420.96990.17170.98600.65470.65420.67490.96931.04201.04360.9580
(n, p) = (105, 0.1)0.09950.09970.09090.04660.00450.06830.06820.07160.00910.09660.09680.0886
(U, L, T) = (11, 0, 6)1.00880.95830.87210.19050.99120.68180.68160.71740.97970.97580.92680.8487
(n, p) = (113, 0.05)0.09800.09250.08260.03210.00340.07070.07050.07620.00690.09450.08920.0783
(U, L, T) = (18, 1, 11)1.09661.09811.01720.35140.99460.74270.74190.76950.98641.06491.06620.9915
(n, p) = (105, 0.1)0.10370.10390.09530.10190.00250.07560.07550.07960.00490.09970.09990.0908
C ^ pmk C ^ pTk P ^ p P ^ pk
F-TChenQF-TChenQF-TChenQ
(U, L, T) = (13, 0, 5)0.86780.86720.91790.93501.09651.04800.98010.87510.87580.9372
(n, p) = (113, 0.05)0.08950.08920.09490.04760.08480.07890.07010.07050.06930.0730
(U, L, T) = (22, 1, 10)1.06691.06541.12000.93211.25591.25691.19411.07051.06901.1284
(n, p) = (105, 0.1)0.10540.10520.11230.05240.09280.09290.08630.08290.08270.0877
(U, L, T) = (10, 0, 5)0.56090.56090.58030.91110.94170.89320.80540.56540.56620.5921
(n, p) = (113, 0.05)0.06390.06380.06920.04520.07280.06720.05760.05060.05010.0543
(U, L, T) = (17, 1, 10)0.64840.64790.66690.91181.04591.04740.96550.65050.65000.6718
(n, p) = (105, 0.1)0.06970.06950.07370.05200.07730.07750.06970.05590.05580.0590
(U, L, T) = (11, 0, 6)0.65870.65840.69750.65090.99520.94670.86520.67250.67330.7117
(n, p) = (113, 0.05)0.06030.06000.06530.06830.07700.07130.06180.05710.05640.0612
(U, L, T) = (18, 1, 11)0.72050.71970.74930.70161.08961.09101.01260.73790.73720.7659
(n, p) = (105, 0.1)0.06530.06520.06910.07090.08050.08070.07310.06120.06110.0647
Table 5. The results of the simulation study for negative binomial processes (sample size m = 50 ).
Table 5. The results of the simulation study for negative binomial processes (sample size m = 50 ).
C ^ p C ^ pc C ^ py C ^ pk C ^ pyk C ^ pm
AnscB-CAnscB-CAnscB-C
(U, L, T) = (358, 5, 85)1.10261.27261.03571.00050.80510.78730.99891.0889
(r, p) = (3, 0.03)0.15570.46390.17320.00080.10900.16330.00130.15390.4484
(U, L, T) = (1428, 76, 462)1.05571.15420.95031.00000.84180.79830.99871.04291.1345
(r, p) = (5, 0.01)0.14590.79830.15440.00070.11540.15760.00150.14420.2740
(U, L, T) = (318, 5, 85)1.07021.23610.64430.99860.74030.73300.99591.05691.2138
(r, p) = (3, 0.03)0.15120.46980.18770.00200.10080.15030.00400.14940.4543
(U, L, T) = (1428, 26, 462)1.42421.81182.88631.00120.84180.82950.99981.40691.7774
(r, p) = (5, 0.01)0.19681.07152.34630.00100.11530.18570.00200.19451.0388
(U, L, T) = (328, 5, 100)1.07861.24550.73920.99920.75720.74760.99701.01751.1667
(r, p) = (3, 0.03)0.15240.46820.19150.00160.10290.15380.00310.14850.4169
(U, L, T) = (1428, 36, 400)1.31481.57762.59981.00120.84180.82470.99981.26351.5175
(r, p) = (5, 0.01)0.18170.73481.71120.00100.11540.18060.00200.15880.7202
C ^ pmk C ^ pTk P ^ p P ^ pk
AnscB-CAnscB-CAnscB-C
(U, L, T) = (358, 5, 85)0.79480.77470.90791.08791.24710.79440.7752
(r, p) = (3, 0.03)0.10490.16260.07090.12920.42750.08900.1541
(U, L, T) = (1428, 76, 462)0.83110.78570.91121.04661.14090.83450.7911
(r, p) = (5, 0.01)0.11070.15740.06750.11730.25540.09220.1465
(U, L, T) = (318, 5, 85)0.73070.72120.90641.05591.21110.73040.7217
(r, p) = (3, 0.03)0.09680.14940.07050.12540.43400.08250.1415
(U, L, T) = (1428, 26, 462)0.83110.81650.91221.41191.78540.83450.8221
(r, p) = (5, 0.01)0.11070.18510.06770.15831.01920.09220.1754
(U, L, T) = (328, 5, 100)0.71250.70420.75221.06431.22040.74710.7360
(r, p) = (3, 0.03)0.08620.14930.10780.12640.43220.08420.1449
(U, L, T) = (1428, 36, 400)0.81050.78900.78371.30341.55630.83450.8173
(r, p) = (5, 0.01)0.11310.16620.09750.14610.69450.09220.1702
Table 6. The results of the simulation study for negative binomial processes (sample size m = 100 ).
Table 6. The results of the simulation study for negative binomial processes (sample size m = 100 ).
C ^ p C ^ pc C ^ py C ^ pk C ^ pyk C ^ pm
AnscB-CAnscB-CAnscB-C
(U, L, T) = (358, 5, 85)1.09131.17341.07471.00070.79840.79900.99941.08331.1619
(r, p) = (3, 0.03)0.10620.21910.12470.00040.07440.12460.00070.10670.2134
(U, L, T) = (1428, 76, 462)1.04581.08750.98211.00000.83450.82360.99921.03811.0769
(r, p) = (5, 0.01)0.10120.15780.11390.00040.07910.12360.00080.10130.1536
(U, L, T) = (318, 5, 85)1.05921.13990.65100.99880.73420.73700.99671.05151.1287
(r, p) = (3, 0.03)0.10310.22370.14570.00130.06880.11220.00260.10360.2178
(U, L, T) = (1428, 26, 462)1.41081.57202.50481.00140.83450.84221.00001.40051.5552
(r, p) = (5, 0.01)0.13650.53401.33330.00070.07910.14560.00130.13660.5186
(U, L, T) = (328, 5, 100)1.06761.14860.75350.99950.75090.75340.99781.01031.0837
(r, p) = (3, 0.03)0.10390.22250.14860.00100.07020.11550.00190.10220.1991
(U, L, T) = (1428, 36, 400)1.30251.41322.35471.00130.83450.83981.00001.25991.3681
(r, p) = (5, 0.01)0.12600.38091.11230.00060.07910.14220.00130.11090.3815
C ^ pmk C ^ pTk P ^ p P ^ pk
AnscB-CAnscB-CAnscB-C
(U, L, T) = (358, 5, 85)0.79230.79180.93371.07611.15260.78720.7860
(r, p) = (3, 0.03)0.07240.12530.05160.08910.19970.06130.1185
(U, L, T) = (1428, 76, 462)0.82820.81610.93661.03731.07780.82770.8170
(r, p) = (5, 0.01)0.07680.12440.04830.08220.14070.06410.1161
(U, L, T) = (318, 5, 85)0.72860.73040.93201.04441.11960.72400.7250
(r, p) = (3, 0.03)0.06680.11260.05140.08650.20480.05680.1063
(U, L, T) = (1428, 26, 462)0.82820.83450.93761.39931.55640.82780.8354
(r, p) = (5, 0.01)0.07680.14620.04860.11090.51010.06410.1388
(U, L, T) = (328, 5, 100)0.70980.71250.75431.05271.12810.74050.7411
(r, p) = (3, 0.03)0.05930.11340.07710.08720.20350.05800.1096
(U, L, T) = (1428, 36, 400)0.80800.81000.78251.29191.39960.82780.8331
(r, p) = (5, 0.01)0.07810.12870.07140.10240.36050.06410.1353
Table 7. The number of nonconformities in samples of 100 printed circuit boards.
Table 7. The number of nonconformities in samples of 100 printed circuit boards.
Sample number12345678910
Nonconformities2124161215528203125
Sample number11121314151617181920
Nonconformities20241619101713221839
Sample number21222324252627282930
Nonconformities30241619171516181215
Sample number31323334353637383940
Nonconformities24212820251918211622
Sample number414243444546
Nonconformities19121491621
Table 8. Indices for the Example 1 data.
Table 8. Indices for the Example 1 data.
Technique C p C pk C pm C pmk C py C pyk C pTk P p P pk
Anscombe’s1.20851.04871.19551.03741.00261.00250.93311.14270.9916
Freeman-Tukey’s1.20741.04951.19431.0382 1.14160.9924
Q1.16681.10531.15001.0895 1.10421.0460
Table 9. The number of nonconforming cans in samples of size n = 50 cans.
Table 9. The number of nonconforming cans in samples of size n = 50 cans.
Sample number12345678910
Number of nonconforming cans8756452347
Sample number11121314151617181920
Number of nonconforming cans65537961043
Sample number21222324252627282930
Number of nonconforming cans58119735214
Sample number31323334353637383940
Number of nonconforming cans5376446856
Table 10. Indices for different transformation methods for the Example 2 data.
Table 10. Indices for different transformation methods for the Example 2 data.
Technique C p C pk C pm C pmk C py C pyk C pTk P p P pk
Freeman-Tukey’s1.05000.97731.04040.96840.99890.95820.93500.87250.8121
Chen’s1.04930.97791.03980.9691 0.87190.8125
Q1.04650.89941.03080.8859 0.85950.7387
Table 11. The number of inspected products until r = 5 nonconforming items are found.
Table 11. The number of inspected products until r = 5 nonconforming items are found.
712288118273747433945
3010533102493115381865
61593073396934552969
994338563828161410662
61244824483958204629
46303962773143361922
45352063433745366856
90147365502723602743
36772881503567194741
24282858366131296285
Table 12. Indices for different transformation methods for the Example 3 data.
Table 12. Indices for different transformation methods for the Example 3 data.
Technique C p C pk C pm C pmk C py C pyk C pTk P p P pk
Anscombe’s1.00460.79260.99630.78611.00100.99940.86311.07600.8489
Q1.51080.59031.46390.5720 1.57720.6162
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Alevizakos, V. Process Capability and Performance Indices for Discrete Data. Mathematics 2023, 11, 3457. https://doi.org/10.3390/math11163457

AMA Style

Alevizakos V. Process Capability and Performance Indices for Discrete Data. Mathematics. 2023; 11(16):3457. https://doi.org/10.3390/math11163457

Chicago/Turabian Style

Alevizakos, Vasileios. 2023. "Process Capability and Performance Indices for Discrete Data" Mathematics 11, no. 16: 3457. https://doi.org/10.3390/math11163457

APA Style

Alevizakos, V. (2023). Process Capability and Performance Indices for Discrete Data. Mathematics, 11(16), 3457. https://doi.org/10.3390/math11163457

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop