Next Article in Journal
A New Instrumental-Type Estimator for Quantile Regression Models
Next Article in Special Issue
Use of Statistical Process Control for Coking Time Monitoring
Previous Article in Journal
Numerical Investigation for Nonlinear Thermal Radiation in MHD Cu–Water Nanofluid Flow in a Channel with Convective Boundary Conditions
Previous Article in Special Issue
Control Chart T2Qv for Statistical Control of Multivariate Processes with Qualitative Variables
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Study on the Laney p′ Control Chart with Parameters Estimated from Phase I Data: Performance Evaluation and Applications

Department of Industrial Engineering and Management, Yuan Ze University, No. 135, Yuan-Tung Road, Chung-Li District, Taoyuan City 32003, Taiwan
*
Author to whom correspondence should be addressed.
Mathematics 2023, 11(15), 3411; https://doi.org/10.3390/math11153411
Submission received: 29 June 2023 / Revised: 1 August 2023 / Accepted: 3 August 2023 / Published: 4 August 2023
(This article belongs to the Special Issue Statistical Process Control and Application)

Abstract

:
The Laney p′ control chart is a new type of attribute control chart that can be applied in situations where the process exhibits either overdispersion or underdispersion. While it has gained acceptance in the industry, there is still limited knowledge about its effectiveness in detecting process variation. It is well known that before applying a control chart, understanding its performance is crucial, especially when the parameters of the control chart need to be estimated from historical data. In this study, we used simulations to investigate the ability of the Laney p′ control chart to detect process variations when the parameters are estimated. We designed appropriate experiments to assess the impact of overdispersion on the average run length (ARL) performance. In this study, we assumed that the overdispersion comes from the variation in the mean fraction nonconforming of each sample. The mean value varies according to a uniform distribution. This study evaluated the performance of the Laney p control chart using the average of the ARL (AARL) and the standard deviation of the ARL (SDARL). Additionally, real-world data were utilized to illustrate the practical applications of the Laney p control chart in the PCB and IC substrate industries. The research findings can serve as valuable guidance for practical implementation.

1. Introduction

Control charts are a statistical process control (SPC) tool used to monitor a process over time and determine whether a manufacturing or business process is in control or out of control. An in-control process refers to a process that is stable, predictable, and exhibits only random cause variations [1]. Conversely, an out-of-control process may display various nonrandom patterns on the control chart, indicating the presence of assignable causes that are responsible for variations in process performance [1]. By understanding the underlying causes, organizations can make informed decisions to improve processes, reduce variation, and enhance product or service quality.
Generally, control charts can be divided into two types: variables control charts and attribute control charts [1]. The specific type of control chart used depends on the nature of the data being monitored. A variable control chart is used to monitor and control process variability when the data being measured are quantitative and continuous. On the other hand, an attribute control chart is a tool used to monitor and control process quality characteristics that are discrete or qualitative in nature. In SPC, the p control chart is an attribute control chart used to monitor the fraction of nonconforming items in a sample [1]. The sample fraction nonconforming (also known as nonconforming rate) is defined as the ratio of the number of nonconforming units in a sample (denoted as X ) to the sample size (denoted as n ). The p control charts have been widely used in monitoring manufacturing processes. Recently, the use of p control charts is increasing in the healthcare industry [2,3,4,5].
This section provides a brief overview of the p chart. For more details, please refer to Montgomery [1]. The development of the p chart requires estimating the process fraction nonconforming, p , from a set of historical data because the value of p is typically unknown. This estimation is carried out by selecting m preliminary samples (also known as subgroups), each with a size of n . Denote p ¯ as the average of these sample fractions nonconforming. The statistic p ¯ is then used to estimate the unknown fraction nonconforming, p . Based on this estimation, the parameters of the p control chart can be expressed as follows:
UCL = p ¯ + 3 p ¯ 1 p ¯ n
CL = p ¯
LCL = p ¯ 3 p ¯ 1 p ¯ n
If the true fraction nonconforming in the production process is known or is a specified standard value, p , then the p ¯ value in Equations (1)–(3) can be replaced with p .
The above description represents the two-phase implementation of control charts. In Phase I, a set of process data is collected and used to estimate the unknown parameters when the process is in control. These estimates are then used to construct the Phase II chart, which monitors the process by plotting the chart statistics and comparing their values to the control limits. It is well known that the more accurate the estimates of the control chart parameters in Phase I, the better the Phase II control chart will perform.
In certain cases where the p control chart is used to monitor the fraction of nonconforming items, the sample consists of a complete inspection of the process output during a specific timeframe. Due to variations in the number of units produced during each period, this leads to variable control limits on the chart.
In the past, numerous studies have been conducted to enhance the traditional p control chart. One area of research involves exploring different control charts, such as MA [6], EWMA [7], and CUSUM [8], to improve the ability to detect changes in the nonconforming rate. Another line of research aims to address the issue of ARL-biased characteristics that emerge when using traditional p control charts to detect changes in the nonconforming rate [9,10,11]. A third category of research focuses on employing supplementary rules to enhance the detection capability of traditional p control charts [9,12,13]. The final research direction aims to modify the traditional p control chart to address the issue of variation in nonconforming rates between samples. This aspect is also the primary focus of our study.
It is important to note that the traditional p control chart assumes that the samples are independent and that the probability of nonconforming is constant over time. If these assumptions are not met, it may be necessary to consider alternative control charts. For instance, in the healthcare application of SPC methods, one issue that raises concern is that of overdispersion. Overdispersion commonly arises in the cases when sample sizes are very large and the parameter p (event probability) is not constant but changes over time.
The Laney p chart [14] is a new type of attribute control chart that is particularly useful in scenarios where there are large sample sizes and the data exhibit overdispersion. Overdispersion can lead to false indications of out-of-control points on a traditional p chart. The Laney p chart incorporates both the variation within samples and the variation between consecutive samples in its definition of common cause variation. When there is overdispersion, the control limits on a Laney p chart are wider in comparison to those of a traditional p chart. The wider control limits mean that only significant deviations in the process are identified as out of control.
The Laney p chart is also useful in the situation when data exhibit underdispersion. Underdispersion, which can occur with samples of any size, is often caused by a lack of randomness. Underdispersion can result in control limits that are too wide for the data. The Laney p chart corrects for underdispersion by calculating narrower control limits.
The calculations for the Laney p chart include σ Z , which is an adjustment for overdispersion or underdispersion. In order to compute σ Z , the sample statistic p i is first standardized using Equation (4):
Z i = p i p ¯ σ p i
where
σ p i = p ¯ 1 p ¯ n i
The value of σ Z can be estimated as
M R ¯ z d 2
where d 2 = 1.128 if a moving range of two observations is used.
The average of the moving range of Z i can be obtained by
M R ¯ Z = i = 1 m Z i Z i 1 m 1
where m is the number of samples. The parameters for the Laney p chart are
UCL = p ¯ + 3 σ p i × σ z
CL = p ¯
LCL = p ¯ 3 σ p i × σ z
A σ Z value of 1 indicates that no adjustment is necessary and that the Laney p control chart is the same as a traditional p chart.
Equations (8) and (10) reveal that the utilization of a p chart incorporates an extra variance component resulting from the variability in parameter p over time, in addition to the conventional calculation of binomial sampling variance.
After conducting a comprehensive literature review, it is evident that a research gap exists regarding the impact of parameter estimation on the performance of the p control chart. Previous studies have primarily focused on evaluating performance solely based on the average run length (ARL). However, when parameters are estimated, the ARL becomes a random variable, necessitating the consideration of both the average value of the average run length (AARL) and the standard deviation of the average run length (SDARL) for performance evaluation [15,16]. This is crucial because different practitioners may obtain varying estimates of the process parameters, leading to different ARL values. Consequently, practitioner-to-practitioner variability introduces randomness into the ARL.
The main objective of this study is to evaluate the performance of the p chart with estimated parameters based on the AARL and SDARL metrics. We also investigate how the number of samples in Phase I influences the performance of the control chart when the parameters are unknown. Furthermore, we evaluate the performance of the p chart when the values of sample size and nonconforming rate are allowed to vary from sample to sample.
The rest of the paper is organized as follows. Section 2 provides a review of previous studies on the Laney attribute control charts, as well as studies on the performance evaluation of the Laney p control charts, where the parameters are estimated from historical data. Section 3 outlines the experimental setup and simulation procedure employed in this study. Section 4 presents the experimental results. In Section 5, we utilize real data to explain the applications of the Laney p control chart in the PCB and IC substrate industries. Finally, Section 6 summarizes the findings and proposes potential avenues for future research.

2. Related Work

This section presents a literature review on topics that are relevant to this study, which includes the application of Laney control charts, as well as the research on evaluating the performance of the p control chart.
In the past, many researchers have studied the impact of overdispersion or underdispersion on data and proposed solutions. Jones and Govindaraju [17] pointed out that overdispersion refers to the variability in data exceeding the variability assumed by their probability distribution. They proposed a simple graphical method to verify the distributional assumption of attribute control charts. This method can also indicate whether there is excessive dispersion or underdispersion in the data. They also briefly reviewed the literature on handling process attributes that do not follow the binomial or Poisson distribution.
Laney [14] introduced the p control chart as a means to prevent negative outcomes resulting from overdispersion in process data. This new type of attribute control chart can also be utilized in cases where a process exhibits insufficient dispersion. Laney demonstrated the application of the p control chart using examples from Heimann [18] but did not conduct performance testing on this chart. Since its introduction, the Laney p control chart has been widely used in healthcare quality monitoring, particularly in cases with large sample sizes. Mohammed and Laney [19] further applied the p control chart to monitor overdispersion in healthcare performance data, highlighting its relevance in healthcare quality monitoring.
Sellers [20] proposed a generalized statistical control chart that can be used to monitor count data for overdispersion or underdispersion. The distributions considered in this study include the binomial, Poisson, and negative binomial distributions.
Vidmara and Blagus [21] pointed out that the identification of outliers in overdispersed proportion data is crucial for effective healthcare quality monitoring. Their study found that the Laney method produces the lowest rate of false alarms. However, in situations where the sample size is small and the proportion is very low, it is difficult to detect outliers, or when the proportion value is very high, outliers are difficult to detect irrespective of sample size.
Evaluating the microbiological quality of pharmaceutical products is an important criterion for determining whether they can be safely released into the drug market. Eissa [22] stated that since the data of the total viable count (TVC) in pharmaceuticals do not follow any specific distribution type, the Laney U control chart is suitable for application in such cases. Eissa [22] believes that when there is significant overdispersion or underdispersion in the data distribution, the Laney U control chart pattern is highly appropriate for monitoring the microbiological characteristics of pharmaceutical products.
Moon [23] applied the Laney U control chart to schedule performance management. In this study, planned value costs were considered as the sample size, and each dollar was assumed to represent an attribute (inspection) unit. Earned value cost, which represents the value of work completed to date, was treated as a variable. By dividing the earned value cost by the planned value cost, a schedule performance index (SPI) was obtained. The study treated SPI as a sample statistic and monitored its variation using the Laney U control chart.
Arafah [24] utilized the Laney p control chart to monitor the variations in the COVID-19 cases in Jordan. The study defined the infection rate (IR) as the number of confirmed cases divided by the number of polymerase chain reaction (PCR) tests conducted. The objective of the research was to understand the effectiveness of the government’s restrictive measures in controlling the infection rate during the COVID-19 pandemic, which included restrictions on population movement and activities. The study results demonstrated that implementing restrictive measures effectively reduced the infection rate, while relaxing these measures had the opposite effect.
Valdés-Manuel and Cogollo-Flórez [25] pointed out that overdispersion is a phenomenon that frequently occurs in large-sample data analysis. In the analysis of discrete data, it refers to a higher level of variation in the data than what is implied by the reference binomial or Poisson distribution. In clinical laboratories, there is often high variability in the proportion of nonconforming units, leading to the occurrence of overdispersion. Therefore, it is necessary to analyze and identify the most suitable control chart to overcome the limitations of traditional control charts when dealing with overdispersed data.
With the introduction and application of the Laney p control chart, some researchers have recently begun to investigate its performance. Ahsan et al. [26] used simulation to evaluate and compare the in-control and out-of-control average run lengths for the p and Laney p control charts. Their study aimed to determine whether the performance of the Laney p control chart is superior to that of the p control chart.
Hagan and Li [27] evaluated the Phase II performance of the Laney p control chart using a simulation approach. Their study allowed for variation in the sample size ( n ) and process nonconforming rate ( p ) between samples. The study assumed that n and p followed truncated normal distributions with a mean of μ and a standard deviation of μ / 5 (or μ / 3 ). When p was allowed to vary across subgroups, its simulated value of μ was either 0.1 or 0.5. When n was allowed to vary, its simulated value of μ varied from 10 to 2000. The performance metric was ARL.
The studies of [26,27] assumed that the process parameters are known and evaluated the performance of the control chart in Phase II. In other words, these studies did not take into account the impact of parameter estimation errors on the performance of the control chart. Furthermore, it is worth noting that the study conducted by Ahsan et al. [26] did not consider the variation in the fraction of nonconforming items between samples.
Recently, Goedhart and Woodall [28] proposed a new method for calculating the control limits of the p chart. They demonstrated that the Laney p chart does not perform well when there are significant variations in subgroup sizes. Laney’s method [14] only yields appropriate limits when there is no intersubgroup variation. Additionally, they pointed out that the estimated control limits for the p chart are excessively wide for smaller sample sizes and too narrow for larger sample sizes. They showed that their method is able to handle situations involving varying subgroup sizes and intersubgroup variation.

3. Methods

3.1. Performance Metrics

The detection capability of a control chart is usually evaluated using the average run length (ARL). The run length refers to the number of samples needed until the chart indicates an out-of-control condition. The ARL represents the expected value of the run length. There are two types of ARL: in-control ARL ( ARL 0 ) and out-of-control ARL ( ARL 1 ). A control chart with a smaller in-control ARL is more sensitive or responsive to detecting process shifts, but it also has a higher probability of Type I errors. Conversely, a control chart with a larger out-of-control ARL is less sensitive or slower to detect process shifts, resulting in a higher probability of Type II errors. A well-designed control chart should have a small ARL 1 when the process undergoes changes, while maintaining a large ARL 0 when there are no changes in the process.
When the control chart parameters are known, the ARL is a constant value. For example, in the case of using a 3-sigma control limit for the x ¯ chart, assuming a normal distribution, the Type I error probability α is 0.0027. In this case, the ARL 0 is 1 / α , which is equal to 370.37. This means that, on average, the control chart will signal a false alarm every 370.37 samples when the process is actually in control.
When the parameters of the control chart are estimated from the samples, the ARL becomes a random variable due to the influence of Phase I sampling [15]. When control charts are constructed using estimated parameters, the average of the ARL (AARL) is commonly used as a measure of performance evaluation. However, AARL alone does not capture the other important properties of ARL. Another significant measure is the standard deviation of the ARL (SDARL), which reflects the practitioner-to-practitioner variability [16]. Practitioner-to-practitioner variability refers to the differences in Phase I data collected by different users, resulting in different parameter estimates and consequently affecting the performance of the control chart. Jones and Steiner [29] first proposed the SDARL metric to study the estimation effect on the performance of the risk-adjusted CUSUM chart. Many other authors used it in evaluating control chart performance when the parameters are estimated [30,31,32]. In this study, we employed both the AARL and SDARL metrics when assessing the performance of the p chart with estimated parameters.
When the process parameters are estimated, in order to achieve a stable in-control ARL ( ARL 0 ), we need a sufficient amount of Phase I sample data so that the in-control AARL value can approach ARL 0 and should have a smaller value of SDARL. Zhang et al. [32] suggested that the SDARL value should be within 10% of the expected in-control ARL ( ARL 0 ). In general, for traditional variable control charts, when the process parameters are known, it is possible to obtain the desired ARL 0 by selecting appropriate control chart parameter combinations [15,32]. However, for p and p control charts, where the sample sizes and nonconforming rates may vary across samples, it becomes more challenging to estimate the expected ARL 0 accurately. Therefore, the criteria proposed by Zhang et al. [32] may not necessarily be fully applicable to the performance evaluation process of p and p control charts in this study.

3.2. Experimental Setup

To investigate the impact of estimating process parameters, we estimated the unknown parameters based on m Phase I samples of size n . The estimated parameters include the in-control process mean p and σ Z . In our study, we considered values of m ranging from 20 to 10,000. It is important to emphasize that the majority of the nonconforming rates analyzed in this study are below 0.05, aiming to reflect the typical conditions found in the general manufacturing industry. In contrast, Hagan and Li [27] examined nonconforming rates of 0.1 and 0.5, which were considered relatively high for the general manufacturing industry. After Phase I analysis, we constructed the Phase II control chart by substituting the unknown parameters with their estimates derived from the generated Phase I samples. We then evaluated the average run length (ARL) value using a simulation approach, which involves 10,000 simulation runs. Finally, we calculate the average ARL (AARL) and standard deviation of ARL (SDARL) by performing 10,000 simulation iterations, each with a unique set of Phase I samples.
In this study, we assume that the between-sample variation comes from the variation in the mean fraction nonconforming of each sample. The mean varies according to a uniform distribution, denoted as U( a ,   b ). A value sampled from the interval a to b is then used to simulate the binomial random variates. Assuming that the number of samples collected in Phase I is m , the simulation procedure used in this study includes the following steps:
  • Define the range of variation for the in-control (process under normal conditions) nonconforming rate p as a uniform distribution U p m i n ,   p m a x and the range of variation for n as a uniform distribution U n m i n ,   n m a x .
  • Using the simulated values of p i and n i , we generate m independent X i values from a binomial distribution with parameters p i and n i . Then, we calculate p i = X i / n i .
  • Calculate the average of the sample fraction nonconforming:
    p ¯ = i = 1 m X i i = 1 m n i .  
    The statistic p ¯ estimates the unknown fraction nonconforming p .
  • Calculate the control limits for the Laney p control chart using Equations (4)–(10).
The above steps describe how to generate m sets of sample data for Phase I analysis. Next, we generate Phase II data and evaluate the run length based on the control limits mentioned above. Evaluating the ARL performance of the p or p control charts is not a simple task. Both the sample size, n , and the process fraction nonconforming, p , will affect the characteristics of the ARL. Additionally, the ARL of the p control chart exhibits ARL-biased characteristics [9]. Under this characteristic, the in-control ARL is not the highest point on the entire ARL curve. Figure 1 illustrates the ARL curve for the p control chart with ARL-biased characteristics. The y -axis of the ARL curve in Figure 1 is displayed on a logarithmic scale. In this example, the in-control fraction nonconforming is 0.04, and the sample size for each group is n = 400 , with ARL 0 = 268.08 . The traditional x ¯ control chart exhibits ARL-unbiased characteristics. When the process undergoes abnormal changes, its ARL (i.e., out-of-control ARL 1 ) will be lower than the ARL when the process is in a normal state (i.e., in-control ARL 0 ).
Figure 2 illustrates the process for estimating the AARL and SDARL in this study. For each set of Phase I data, this study conducted 10,000 experiments ( j = 10,000) to estimate the ARL. We computed a specific set of control limits for each Phase I dataset. This means that throughout the 10,000 simulation runs conducted to obtain an ARL value, the same control limits were applied in each experiment. As the data collected from Phase I may vary, it reflects the variability among practitioners. This study generated 10,000 sets of Phase I data ( i = 10,000) to obtain the AARL and SDARL values.

4. Results and Discussion

The aim of this study is to assess the performance of the Laney p control chart through various experiments. In this section, we present the results obtained from the simulation experiments. We investigate the effect of between-sample variation on the ARL performance.
Firstly, we consider the performance of the traditional p control chart and the Laney p control chart under different sample sizes ( n ) when the nonconforming rate varies between samples. In this experiment, we assume p ~ U 0.045 ,   0.055 , which means the average is 0.05. When the sample sizes are 100, 150, 225, 300, and 350, the theoretical ARL 0 values of the p chart are 233.96, 277.54, 422.76, 365.86, and 279.28, respectively.
The results of the p chart are presented in Table 1, while Table 2 displays the results of the Laney p chart. In general, as the number of Phase I samples increases, the SDARL becomes smaller. The results show that the use of a small number of Phase I samples to estimate the unknown parameters leads to a huge variability in the ARL values, resulting in a very poor chart performance. In this specific case, the variation in the nonconforming rate 𝑝 is not considered significant. It can be observed that the in-control AARL of the 𝑝 control chart and Laney p control chart does not differ significantly when the number of samples 𝑚 is equal to or greater than 1000. However, when the number of samples 𝑚 is small, both the AARL and SDARL of the Laney p control chart will be higher compared to those of the traditional 𝑝 control chart.
When the sample size 𝑛 is large, even if there is not much variation in the nonconforming rate 𝑝, the Laney p control chart will still exhibit higher AARL and SDARL values compared to the traditional 𝑝 control chart. This is evident by comparing column (e) of Table 1 to column (e) of Table 2.
The results from Table 2 also highlight that when the sample size ( n ) is small ( n 300 ), the AARL of the Laney p control chart does not reach the ARL 0 value of the traditional control chart. This observation further illustrates that the Laney p control chart may not adequately widen the control limits when the sample size is relatively small. However, when the sample size ( n ) is increased to 350, the AARL of the Laney p control chart exceeds the theoretical ARL 0 value when the number of samples m is 10,000. This indicates that the control limits are excessively widened in this scenario.
Finally, there are some unusual phenomena observed in column (c) of Table 1 and Table 2 that need to be explained. By observing Table 1 and Table 2, it can be seen that under the combination of n and p in column (c) of Table 1, the SDARL is noticeably larger than other combinations. Furthermore, column (c) of Table 2 shows that the SDARL for the Laney p control chart is smaller than the traditional p control chart. After conducting an in-depth investigation, we found that in column (c), under the combination of n and p , there may or may not be a lower control limit. This causes an increase in the variability in in-control ARL, resulting in an increase in SDARL. For columns (a) and (b) of Table 1, there are no lower control limits under the combination of n and p . Furthermore, for columns (d) and (e) of Table 1, both combinations of n and p have lower control limits. These phenomena result in smaller variations in ARL for columns (a), (b), (d), and (e), indicating smaller values of SDARL.
Next, we consider the impact on the performance of traditional p control charts and Laney p control charts when the sample size n is fixed, and the nonconforming rate varies within different ranges. We consider the case where the true fraction nonconforming p = 0.05 and sample size n = 300 . In this case, the ARL 0 = 365.86 . Table 3 presents the in-control AARLs and SDARLs of the p control chart when the process parameters are estimated. On the other hand, Table 4 shows the in-control AARLs and SDARLs for the Laney p chart.
From panel (a) of Table 3, when the number of samples, m , exceeds 3000, the AARL of the p control chart approaches the theoretical ARL 0 value of 365.86. This demonstrates the high accuracy of the simulation program developed in this study.
Although the Laney p control chart is suitable for situations where the nonconforming rate of the process is excessively dispersed, one may be interested in investigating its performance when there is no variation in sample size n and process p . If m is small, one can see from panel (a) of Table 4 that the p chart may tend to overadjust the control limits, leading to the AARLs higher than the ARL 0 when the fraction nonconforming remains constant. When m is greater than 2000, the AARL is close to the true ARL value. This can be explained by the fact that the estimate of σ Z is more accurate when m increases. When there is neither overdispersion nor underdispersion in the nonconforming rate, and if the number of Phase I samples ( m ) is large, the AARL of the p control chart (panel (a) of Table 4) is close to that of the p control chart (panel (a) of Table 3). However, the p control chart exhibits a higher SDARL. This is because the p control chart requires an additional estimation of σ Z , leading to greater variability in its performance.
When the control chart parameters are estimated from Phase I samples, it is important to investigate the influence of the number of samples, m , on the size of the SDARL. When using control charts, it is generally recommended to have at least 20 or 25 Phase I samples to estimate the unknown parameters [1]. However, in this particular case, even with a significantly larger number of samples exceeding 20, the desired performance is not achieved. For instance, when using 50 samples, the AARL of the p chart is 335.73, which corresponds to approximately 91% of the desired ARL 0 . Additionally, the SDARL is 128.22, approximately 35% of the ARL 0 . To attain satisfactory performance, it is necessary to utilize 1000 samples, as this yields an AARL of 362.66 (around 99% of the desired ARL 0 ) and an SDARL of 24.61 (approximately 7% of the desired ARL 0 ). In the case of the Laney p control chart, utilizing 5000 samples is recommended, as it produces an AARL of 358.30 (about 98% of the desired ARL 0 ) and an SDARL of 34.91 (approximately 9.5% of the desired ARL 0 ). It is evident that when using the Laney p control chart, we require a larger number of samples to meet the recommended requirements of Zhang et al. [32].
Next, we investigate the impact of overdispersion in the process fraction nonconforming on the p control chart. Overdispersion refers to the variation in data that exceeds the assumed probability distribution. In this experiment, we assume that the average of the fraction nonconforming is fixed, but the variability in the fraction nonconforming is uniformly distributed.
From the results in Table 3, it can be observed that when the sample size n is fixed, the AARL of the p control chart decreases as the variation in the nonconforming rate increases. From Table 4, it can be seen that the AARL of the p control chart increases as the variability in nonconforming rate increases. The AARL value has exceeded the case when the nonconforming rate follows a binomial distribution (i.e., no overdispersion or underdispersion), indicating that the control limits are excessively widened.
The above results also show that when there is overdispersion in 𝑝, even if the sample size 𝑛 is not very large, it will result in a much lower in-control ARL than the in-control ARL value of the traditional 𝑝 control chart. In other words, using a traditional 𝑝 control chart in this scenario would lead to an elevation in Type I error.
We will now examine situations where the sample sizes differ from those in the previous experiment, while ensuring that the product of sample size and nonconforming rate remains constant. Table 5 and Table 6 display the AARLs and SDARLs of the p chart and Laney p chart with a sample size of 3000 and an average fraction nonconforming of 0.005, respectively. In this scenario, the ARL₀ of the p chart is 290.73.
Furthermore, Table 7 and Table 8 present the AARLs and SDARLs of the p chart and Laney p chart, respectively, with a larger sample size of n   = 30 , 000 and an average fraction nonconforming of p   = 0.0005 . In this case, the ARL₀ of the p chart is 284.51.
Upon comparing Table 3, Table 4, Table 5, Table 6, Table 7 and Table 8, we observe that the conclusions drawn from Table 3 and Table 4 are also applicable to Table 5 and Table 6, as well as Table 7 and Table 8. That is, when the sample size remains constant and the variability in the process nonconforming rate increases, the in-control AARLs of the traditional p control charts decrease. The p chart can address the between-sample variation. The in-control AARLs are close to the advertised values. In some cases, the AARLs may exceed the advertised value.
Next, we consider the scenario where both the sample size, n , and the nonconforming rate, p , are allowed to vary between samples. Table 9 shows the results of the traditional p control chart, while Table 10 presents the results of the Laney p control chart. By comparing Table 5 with Table 9 and Table 6 with Table 10, we can observe that the variation in sample size ( n ) does not have a significant impact on the AARL performance of both the traditional p control chart and the Laney p control chart. However, it is worth noting that when both n and p vary, it can lead to an increased variability in parameter estimation, as evidenced by the increase in SDARL.
Next, we set a fixed number of samples for the Phase I dataset and examine the influence of sample size on the in-control AARL and SDARL for the p control chart. In this experiment, we consider a wider range of variations in the sample size. Based on previous experiments, it has been found that when the number of samples is 1000 or more, the AARL and SDARL of the p control chart yield better results. Therefore, we fix the number of samples at 1000 and investigate the impact of sample size on AARL and SDARL when there are variations in the nonconforming rate between samples. Table 11 summarizes the results of the experiment.
From the table, we can observe that the in-control AARL of the traditional p control chart is very small, far below the theoretical values of the 3-sigma p control chart. This is due to the impact of overdispersion. As the sample size increases, the influence of overdispersion on the in-control AARL becomes more significant. In other words, when there is a high degree of variation in the fraction nonconforming between samples, the reduction in AARL for the p control chart becomes more substantial as n increases.
Finally, we compare the out-of-control AARL and SDARL of the p chart and Laney p chart. The magnitude of shift considered is represented by a multiple k of σ p . If the in-control nonconforming rate varies between U( a ,   b ), the range of variation for the out-of-control mean is U( a ,   b ) + k σ p . After estimating the control limits using a set of Phase I data, the estimation of the out-of-control ARL is performed. The AARL and SDARL performance of the traditional p chart and the Laney p chart are summarized in Table 12 and Table 13, respectively. In this comparison, the number of Phase I samples ( m ) for the traditional p chart and the Laney p chart is 1000 and 3000, respectively, with a fixed sample size of 3000. The selection of m values for the p chart and the Laney p chart is based on the following rationale. According to Table 5, it is evident that when the value of 𝑚 is 1000 or higher, the SDARL for the 𝑝 control chart stays within 10% of the expected in-control ARL. Similarly, Table 6 shows that the Laney p control chart exhibits a similar behavior when the 𝑚 value is 3000 or above. Table 12 and Table 13 also include the in-control AARL and SDARL values for ease of comparison.
From Table 12 and Table 13, it is observed that both the p chart and Laney p chart exhibit a decrease in AARL as the magnitude of shift increases. When the nonconforming rate is kept fixed among samples, from column (a) of Table 12 and Table 13, we can see that the AARLs of the traditional p control chart and the Laney p′ control chart are very close, but the SDARLs of the Laney p chart are still greater than those of the p chart. As explained earlier, this is because the Laney p′ control chart requires an additional estimation of σ Z , resulting in an increased variability in ARL. By comparing columns (b) and (c) of Table 12 and Table 13, we can observe significant differences in the in-control and out-of-control performance between the traditional p control chart and the Laney p′ control chart. This is due to variations in the nonconforming rates among different samples, indicating excessive process dispersion. The datasets generated based on the scenarios described in columns (b) and (c) of Table 12 have been checked using Minitab statistical software (Version 19, Minitab LLC, State College, PA, USA) and the diagnostic method described in [17]. This confirms the presence of overdispersion. Additionally, comparing column (d) of Table 12 and Table 13, we can notice that the in-control and out-of-control performance of the two control charts are very similar. The diagnostic method described in [17] confirms that overdispersion does not exist.
Based on the above comparison, this study summarizes the advantages and disadvantages of the 𝑝 control chart and the Laney p′ control chart as follows. While the calculation of the 𝑝 control chart is relatively simple, it tends to have a higher Type I error when the nonconforming rate varies between samples. On the other hand, the Laney p′ control chart has the disadvantage of requiring an estimation of the variation in the nonconforming rate between samples. This typically requires a larger number of samples to achieve an accurate parameter estimation. When the number of samples 𝑚 is small, the SDARL for the Laney p′ control chart is higher compared to the traditional 𝑝 control chart. If there is minimal variation in the nonconforming rate between samples, the control limits of the Laney p′ control chart closely resemble those of the traditional 𝑝 control chart. However, due to the need for estimating between-sample variation, the variability in its performance also increases. Additionally, the Laney p′ control chart addresses the issue of overdispersion by widening the control limits. This study found that the Laney p′ control chart can indeed reduce Type I errors caused by overdispersion. However, in some cases, this approach may result in excessively wide control limits, thereby reducing its ability to detect changes in the nonconforming rate.

5. Applications of the Laney p Chart in Manufacturing Processes

In previous studies [2,3,4,5], researchers have pointed out that the healthcare industry often experiences the phenomenon of overdispersion due to its extremely large sample sizes. However, overdispersion also frequently occurs in the general manufacturing industry, which offers numerous opportunities for the application of Laney p control charts. This section will utilize real cases and data collected from the industry to explain the reasons for having such large sample sizes and demonstrate how to apply Laney p   control charts to address the issue of overdispersion and obtain reasonable results. The examples are from the PCB and IC subtractive manufacturing processes. In the process of circuit board and chip substrate manufacturing, a panel can be divided into multiple strips, with each strip containing numerous units. Inspection can be conducted either at the strip level or the unit level in these processes. The size of each sample represents the total number of strips (or units) inspected within a week (or a batch). When calculating the nonconforming rate, a large sample size will be involved.

5.1. Example 1

During the ball placement process of a ball grid array (BGA), the operator first uses an automated optical inspection (AVI) machine to scan the bump area of the BGA. The AVI machine compares the scanned data against an image file or predefined criteria to determine if the bumps meet the required specifications. Any captured images showing abnormal bumps are further examined by operators using the verify and recheck system (VRS) machine. The VRS operator views the magnified images displayed on the VRS screen (at a 32× magnification).
If the bump size of the solder balls within the location marked by the AVI is higher or lower than the adjacent solder balls, they are classified as either large-sized balls or small-sized balls. Considering the nature and severity of the product, if the VRS machine confirms the presence of any large or small balls, regardless of their quantity, the unit is considered nonconforming.
The quality engineer responsible for process improvement decided to set up a control chart in order to improve the fraction of nonconforming units produced by this process. In this application, the 100% inspection approach is adopted. The nonconforming rate is calculated by dividing the number of nonconforming units found per lot by the total number of units inspected in each lot. In this specific case, the sample sizes represent the total number of units inspected in each lot. However, it is important to note that this value is not constant. The number of units in each lots varies depending on the part number, typically ranging from 1800 to 7000 units. In the following presentations, the p chart and Laney p chart were produced using Minitab statistical software.
The quality control engineer gathered data from 20 lots to establish a p   control chart. The data are shown in Table 14. Figure 3a presents a traditional p control chart for the 20 lots. The sample fraction nonconforming from each sample is plotted on this chart. The control limits of the 𝑝 control chart vary as a result of employing the 100% inspection method and the variability in daily production output. The data demonstrate a significant amount of overdispersion. The p control chart indicates many points beyond the control limits. However, the analysis of these points does not reveal any reasonable assignable causes for them. The diagnostic method described in [17] confirmed the presence of overdispersion. As a result, the engineer decided to use the Laney p control chart. Figure 3b displays the Laney p control chart, which shows that all points are within the control limits. This suggests a stable process with no assignable causes present. Thus, the engineer has concluded Phase I of the control usage.
Next, the engineer collected Phase II data to monitor the production process. From Figure 3b, it is evident that Sample 25 plots outside the upper control limit, indicating a significant increase in the process nonconforming rate. Subsequently, a process diagnosis must be conducted to identify the assignable causes and initiate improvements. There are several reasons for the generation of large and small balls in this process. The investigation of Sample 25 revealed that the assignable cause is the variation in incoming material. During ball placement, when the substrate magnification matches the steel plate magnification, μ -balls can smoothly fit into the solder mask openings. However, if there is a variation in the magnification of incoming materials (with some substrates being larger or smaller), it can lead to an improper alignment between the steel plate and the substrate. This misalignment causes some balls to fit properly while others do not. As a consequence, the balls cannot accurately land in the solder mask opening, leading to ball misalignment issues.

5.2. Example 2

A PCB lithography process consists of several steps, including preprocessing, lamination, exposure, etching, and finishing. In this process, the nonconforming rate caused by dirt particles is an important quality indicator. The process utilizes a negative-type dry film that polymerizes and hardens when exposed to ultraviolet (UV) light. After lamination is completed, if there are dirt particles adhering to the substrate, the image will be masked during exposure. The unexposed areas will not polymerize or harden, and during development, the dry film is removed (developed), exposing the underlying copper. During the etching stage, the copper layer will be etched away by the etching solution, resulting in nonconformities such as open circuits, nicks, and pinholes.
The quality control engineer collected daily data to create a p control chart with the aim of improving the production process. In this application, the inspection method used is 100% inspection. The current quality control specification states that if any nonconformities are confirmed on a strip, that strip is considered nonconforming. The nonconforming rate is calculated by dividing the number of nonconforming strips found per day by the total number of strips inspected each day. The number of strips in each panel varies depending on the part number, typically ranging from 8 to 12 strips. The production capacity is approximately 28,000 to 35,000 panels per day. To establish a Phase I control chart, preliminary data were collected over a period of 25 days. The data are shown in Table 15, and the traditional p control chart is plotted in Figure 4a. The 𝑝 control chart exhibits variable control limits due to the implementation of 100% inspection and the fluctuations in daily production output. In this specific case, the data demonstrate a significant amount of overdispersion due to the large sample size (200,000 to 400,000 strips per day).
An examination of the p control chart reveals that there are many points outside the control limits. Out of the 25 points, 23 points exceed the control limits. This can lead to a misconception that the process is influenced by many assignable causes and is out of control. The occurrence of data points falling outside the control limits can imply the potential existence of assignable causes that warrant investigation and corrective action. However, after conducting a comprehensive analysis of the process, it was determined that the process remains in control, with no identifiable assignable causes present. This phenomenon has been pointed out by Mohammed and Laney [19], who noted that overdispersion caused by large sample sizes may lead to the identification of an inappropriately large number of data points indicating assignable cause variation.
The engineer confirmed the necessity of using the Laney p’ chart by utilizing the diagnostic method described in [17]. The Laney p control chart is shown in Figure 4b, with control limits adjusted for overdispersion. This control chart reveals that all points fall within the control limits, indicating a stable process without any assignable causes. Based on this observation, the engineer adopted the Laney p chart for Phase II usage.

6. Conclusions

In this study, we employed simulation to investigate the performance of the traditional p and Laney p control charts when the parameters are estimated. Both charts were evaluated and compared in terms of AARL and SDARL. Overdispersion was accounted for by allowing for between-sample variation in the nonconforming rate ( p ) . Additionally, variations between samples in both the sample size ( n ) and the nonconforming rate of the process were also considered.
This study found that the Laney p control chart can indeed reduce Type I errors caused by overdispersion, but it may also result in excessively wide control limits and an increased risk of Type II errors. Our results show that the p chart requires a much larger amount of Phase I data than previously recommended in the literature in order to sufficiently reduce the variation in the chart performance. These research findings can serve as a valuable reference for practitioners.
Future research can be conducted in several directions. Firstly, a broader range of combinations of n and p should be investigated to achieve a more comprehensive understanding. Furthermore, the methodology presented in this study can be extended to evaluate the effectiveness of the Laney U control chart. A more comprehensive assessment of its performance will facilitate its practical implementation and promotion in various industries.

Author Contributions

Conceptualization, C.-S.C.; methodology, C.-S.C. and P.-W.C.; software, C.-S.C., P.-W.C. and C.-W.W.; formal Analysis, C.-S.C., P.-W.C. and C.-W.W.; investigation, C.-S.C. and P.-W.C.; data curation, P.-W.C. and C.-W.W.; writing—original draft preparation, C.-S.C., P.-W.C. and C.-W.W.; writing—review and editing, C.-S.C., P.-W.C. and C.-W.W.; supervision, C.-S.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was partially funded by the National Science and Technology Council, R.O.C. (grant number NSTC 112-2222-E-155-001).

Data Availability Statement

All data generated in this study are available upon request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Montgomery, D.C. Introduction to Statistical Quality Control, 8th ed.; John Wiley & Sons: New York, NY, USA, 2020. [Google Scholar]
  2. Woodall, W.H. The use of control charts in health-care and public-health surveillance. J. Qual. Technol. 2006, 38, 89–104. [Google Scholar] [CrossRef]
  3. Duclos, A.; Voirin, N. The p-control chart: A tool for care improvement. Int. J. Qual. Health Care 2010, 22, 402–407. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Suman, G.; Prajapati, D. Control chart applications in healthcare: A literature review. Int. J. Metrol. Qual. Eng. 2018, 9, 5. [Google Scholar] [CrossRef] [Green Version]
  5. Kam, Y.; Hee, V. A study on the application of control chart in healthcare. ITM Web Conf. 2021, 36, 01001. [Google Scholar]
  6. Khoo, M.B. A moving average control chart for monitoring the fraction non-conforming. Qual. Reliab. Eng. Int. 2004, 20, 617–635. [Google Scholar] [CrossRef]
  7. Gan, F.F. Monitoring observations generated from a binomial distribution using modified exponentially weighted moving average control chart. J. Stat. Comput. Simul. 1990, 37, 45–60. [Google Scholar] [CrossRef]
  8. Gan, F.F. An optimal design of CUSUM control charts for binomial counts. J. Appl. Stat. 1993, 20, 445–460. [Google Scholar] [CrossRef]
  9. Acosta-Mejia, C.A. Improved p charts to monitor process quality. IIE Trans. 1999, 31, 509–516. [Google Scholar] [CrossRef]
  10. Morais, M.C. An ARL-unbiased np-chart. Econ. Qual. Control 2016, 31, 11–21. [Google Scholar] [CrossRef]
  11. Argoti, M.A.; Carrión-García, A. A heuristic method for obtaining quasi ARL-unbiased p-charts. Qual. Reliab. Eng. Int. 2019, 35, 47–61. [Google Scholar] [CrossRef] [Green Version]
  12. Wludyka, P.S.; Jacobs, S.L. Runs rules and p-charts for multistream binomial processes. Commun. Stat. Simul. Comput. 2002, 31, 97–142. [Google Scholar] [CrossRef]
  13. Khilare, S.K.; Shirke, D.T. Fraction nonconforming control charts with m-of-m runs rules. Int. J. Adv. Manuf. Technol. 2015, 78, 1305–1314. [Google Scholar] [CrossRef]
  14. Laney, D.B. Improved control charts for attributes. Qual. Eng. 2002, 14, 531–537. [Google Scholar] [CrossRef]
  15. Hany, M.; Mahmoud, M.A. An evaluation of the Crosier’s CUSUM control chart with estimated parameters. Qual. Reliab. Eng. Int. 2016, 32, 1825–1835. [Google Scholar] [CrossRef]
  16. Saleh, N.A.; Mahmoud, M.A.; Jones-Farmer, L.A.; Zwetsloot, I.; Woodall, W.H. Another look at the EWMA control chart with estimated parameters. J. Qual. Technol. 2015, 47, 363–382. [Google Scholar] [CrossRef]
  17. Jones, G.; Govindaraju, K. A graphical method for checking attribute control chart assumptions. Qual. Eng. 2001, 13, 19–26. [Google Scholar] [CrossRef]
  18. Heimann, P.A. Attributes control charts with large sample sizes. J. Qual. Technol. 1996, 28, 451–459. [Google Scholar] [CrossRef]
  19. Mohammed, M.A.; Laney, D. Overdispersion in health care performance data: Laney’s approach. BMJ Qual. Saf. 2006, 15, 383–384. [Google Scholar] [CrossRef] [Green Version]
  20. Sellers, K.F. A generalized statistical control chart for over-or under-dispersed data. Qual. Reliab. Eng. Int. 2012, 28, 59–65. [Google Scholar] [CrossRef]
  21. Vidmar, G.; Blagus, R. Outlier detection for healthcare quality monitoring–a comparison of four approaches to over-dispersed proportions. Qual. Reliab. Eng. Int. 2014, 30, 347–362. [Google Scholar] [CrossRef]
  22. Eissa, M.E. Application of Laney control chart in assessment of microbiological quality of oral pharmaceutical filterable products. Bangladesh J. Sci. Ind. Res. 2017, 52, 239–246. [Google Scholar] [CrossRef] [Green Version]
  23. Moon, J. An investigation into the use of Laney u chart as a visual schedule tracker to graphically monitor the schedule performance index. J. Eng. Proj. Prod. Manag. 2020, 10, 35–41. [Google Scholar]
  24. Arafah, M. Using the Laney p’ control chart for monitoring COVID-19 cases in Jordan. J. Healthc. Eng. 2022, 2022, 6711592. [Google Scholar] [CrossRef] [PubMed]
  25. Valdés-Manuel, J.I.; Cogollo-Flórez, J.M. Monitoring overdispersed process in clinical laboratories using control charts. DYNA 2022, 89, 28–33. [Google Scholar] [CrossRef]
  26. Ahsan, M.; Mashuri, M.; Khusna, H. Evaluation of Laney p’ chart performance. Int. J. Appl. Eng. Res. 2017, 12, 14208–14217. [Google Scholar]
  27. Hagan, J.; Li, B. Phase II performance of P-charts and P’-charts. J. Med. Stat. Inform. 2018, 6, 3. [Google Scholar] [CrossRef] [Green Version]
  28. Goedhart, R.; Woodall, W.H. Monitoring proportions with two components of common cause variation. J. Qual. Technol. 2021, 54, 324–337. [Google Scholar] [CrossRef]
  29. Jones, M.A.; Steiner, S.H. Assessing the effect of estimation error on risk-adjusted CUSUM chart performance. Int. J. Qual. Health Care 2012, 24, 176–181. [Google Scholar] [CrossRef]
  30. Lee, J.; Wang, N.; Xu, L.; Schuh, A.; Woodall, W.H. The effect of parameter estimation on upper-sided Bernoulli Cumulative Sum charts. Qual. Reliab. Eng. Int. 2013, 29, 639–651. [Google Scholar] [CrossRef]
  31. Zhang, M.; Peng, Y.; Schuh, A.; Megahed, F.M.; Woodall, W.H. Geometric charts with estimated control limits. Qual. Reliab. Eng. Int. 2013, 29, 209–223. [Google Scholar] [CrossRef]
  32. Zhang, M.; Megahed, F.M.; Woodall, W.H. Exponential CUSUM charts with estimated control limits. Qual. Reliab. Eng. Int. 2014, 30, 275–286. [Google Scholar] [CrossRef]
Figure 1. The ARL curve for the p control chart with ARL-biased characteristics.
Figure 1. The ARL curve for the p control chart with ARL-biased characteristics.
Mathematics 11 03411 g001
Figure 2. Illustration of the process for estimating AARL and SDARL in this study.
Figure 2. Illustration of the process for estimating AARL and SDARL in this study.
Mathematics 11 03411 g002
Figure 3. Control charts for Example 1: (a) the traditional p control chart; (b) the Laney p control chart.
Figure 3. Control charts for Example 1: (a) the traditional p control chart; (b) the Laney p control chart.
Mathematics 11 03411 g003
Figure 4. Control charts for Example 2: (a) the traditional p control chart; (b) the Laney p control chart.
Figure 4. Control charts for Example 2: (a) the traditional p control chart; (b) the Laney p control chart.
Mathematics 11 03411 g004
Table 1. The effect of dispersion in p on the in-control AARL of the p chart.
Table 1. The effect of dispersion in p on the in-control AARL of the p chart.
( a )   n   = 100
p ~ U 0.045 ,   0.055
( b )   n   = 150
p ~ U 0.045 ,   0.055
( c )   n   = 225
p ~ U 0.045 ,   0.055
( d )   n   = 300
p ~ U 0.045 ,   0.055
( e )   n   = 350
p ~ U 0.045 ,   0.055
m AARLSDARLAARLSDARLAARLSDARLAARLSDARLAARLSDARL
20327.32345.69317.31225.48292.50177.66273.62137.39256.47121.50
50273.16183.91285.29163.35282.15140.03266.4198.79258.9397.56
100237.1573.63270.68124.10275.80108.30261.0773.98257.0481.71
1000215.752.17246.062.47289.1183.16284.5520.20225.4038.15
2000215.752.17246.062.46299.8978.63287.227.66217.5919.91
3000215.752.17246.062.46304.1176.64287.593.64215.464.76
4000215.752.17246.062.46312.2271.12287.612.87215.243.77
5000215.752.17246.062.46315.8068.64287.412.88215.002.16
6000215.752.17246.052.46319.4765.36287.412.88215.002.16
8000215.752.17246.032.45324.3860.31287.412.88215.002.16
10,000215.752.17246.032.45328.2356.46287.412.87214.982.16
Table 2. The effect of dispersion in p on the in-control AARL of the p chart.
Table 2. The effect of dispersion in p on the in-control AARL of the p chart.
( a )   n   = 100
p ~ U 0.045 ,   0.055
( b )   n   = 150
p ~ U 0.045 ,   0.055
( c )   n   = 225
p ~ U 0.045 ,   0.055
( d )   n   = 300
p ~ U 0.045 ,   0.055
( e )   n   = 350
p ~ U 0.045 ,   0.055
m AARLSDARLAARLSDARLAARLSDARLAARLSDARLAARLSDARL
20557.71773.40615.71817.45637.49848.38633.91843.58648.33856.97
50378.23476.36447.62547.64513.17607.95521.04623.57534.77641.76
100298.40280.70367.56363.29386.91369.49420.47414.75429.11422.95
1000222.7757.94262.7681.59306.1093.02313.45100.04327.46104.59
2000216.5016.25249.3434.35314.3269.83297.9858.23330.5693.62
3000215.894.83246.3312.92323.6362.56293.3340.41331.9293.53
4000215.914.80246.226.43327.5457.79290.8430.57334.5192.04
5000215.904.80246.062.47334.1349.20289.4622.29335.8792.32
6000215.902.17246.052.57337.4443.66288.1113.28339.6191.78
8000215.892.17246.032.52341.4435.82287.658.17346.5189.92
10,000215.892.17246.032.50343.8529.35287.412.87346.8390.15
Table 3. The effect of dispersion in p on the in-control AARL of the p chart with n   = 300 .
Table 3. The effect of dispersion in p on the in-control AARL of the p chart with n   = 300 .
( a )   n   = 300
p   = 0.05
( b )   n   = 300
p ~ U 0.04 ,   0.06
( c )   n   = 300
p ~ U 0.025 ,   0.075
m AARLSDARLAARLSDARLAARLSDARL
20344.24183.03156.4174.2425.588.01
50335.73128.22153.8752.3626.415.98
100334.7999.70153.2340.2126.694.64
1000362.6624.61164.6113.3528.592.37
2000365.1513.49166.656.5229.121.48
3000365.853.63166.962.8429.261.07
4000365.863.63167.052.3029.310.88
5000365.853.63167.062.2929.340.67
6000365.853.63167.061.6829.350.68
8000365.853.63167.061.6829.350.67
10,000365.853.63167.061.6829.350.67
Table 4. The effect of dispersion in p on the in-control AARL of the p chart with n   = 300 .
Table 4. The effect of dispersion in p on the in-control AARL of the p chart with n   = 300 .
( a )   n   = 300
p   = 0.05
( b )   n   = 300
p ~ U 0.04 ,   0.06
( c )   n   = 300
p ~ U 0.025 ,   0.075
m AARLSDARLAARLSDARLAARLSDARL
20437.20457.00453.62430.54887.39811.39
50453.10410.19451.46383.06877.27749.50
100426.80334.35382.48286.55769.08605.44
1000330.5597.61310.7893.00530.47154.53
2000340.7867.87309.3448.73521.18107.97
3000349.0952.39310.3133.40521.8088.73
4000353.9742.88312.7921.22527.6970.79
5000358.3034.91312.9117.90531.2362.25
6000360.3730.12313.8512.83534.8056.26
8000362.8522.74314.179.28541.0746.54
10,000364.3117.57314.516.47546.0837.61
Table 5. The effect of dispersion in p on the in-control AARL of the p chart with n   = 3000 .
Table 5. The effect of dispersion in p on the in-control AARL of the p chart with n   = 3000 .
( a )   n   = 3000
p   = 0.005
( b )   n   = 3000
p ~ U 0.004 ,   0.006
( c )   n   = 3000
p ~ U 0.0025 ,   0.0075
m AARLSDARLAARLSDARLAARLSDARL
20344.11182.17159.4677.6027.198.93
50332.13140.89157.2059.3728.207.14
100326.68119.57154.3749.1328.466.10
1000292.1719.08141.9411.2227.692.37
2000290.892.89140.821.9627.380.99
3000290.722.89140.711.9527.280.60
4000290.882.91140.721.9527.250.28
5000290.722.91140.691.4127.250.28
6000290.722.91140.691.4127.250.27
8000290.722.89140.691.4127.250.27
10,000290.722.89140.691.4127.250.27
Table 6. The effect of dispersion in p on the in-control AARL of the p chart with n   = 3000 .
Table 6. The effect of dispersion in p on the in-control AARL of the p chart with n   = 3000 .
( a )   n   = 3000
p   = 0.005
( b )   n   = 3000
p ~ U 0.004 ,   0.006
( c )   n   = 3000
p ~ U 0.0025 ,   0.0075
m AARLSDARLAARLSDARLAARLSDARL
20434.30462.21423.47416.02515.65422.09
50427.05403.96423.08369.48583.49394.31
100391.80322.65372.33274.48548.70344.05
1000310.4494.60294.0083.64458.03128.31
2000301.0258.22279.7462.55443.3089.45
3000293.3332.16271.7050.38441.0672.52
4000292.9725.08267.0941.90439.2054.13
5000291.5814.71263.1733.06437.7943.28
6000291.1111.27260.3127.67439.1135.13
8000290.795.32258.1617.17437.5421.14
10,000290.774.44257.5813.55436.5016.86
Table 7. The effect of dispersion in p on the in-control AARL of the p chart with n   = 30 , 000 .
Table 7. The effect of dispersion in p on the in-control AARL of the p chart with n   = 30 , 000 .
( a )   n   = 30 , 000
p   = 0.0005
( b )   n   = 30 , 000
p ~ U 0.0004 ,   0.0006
( c )   n   = 30 , 000
p ~ U 0.00025 ,   0.00075
m AARLSDARLAARLSDARLAARLSDARL
20341.54177.11157.7177.9327.618.93
50332.47141.49157.1058.9628.387.19
100288.50120.07154.7048.7328.526.09
1000288.0624.65140.3213.3527.632.46
2000284.856.15138.652.2327.221.40
3000284.632.86138.471.3927.110.90
4000284.622.85138.471.3827.070.67
5000284.632.83138.481.3927.050.60
6000284.632.85138.481.3827.070.60
8000284.592.86138.461.4027.060.60
10,000284.642.82138.481.3827.070.59
Table 8. The effect of dispersion in p on the in-control AARL of the p chart with n   = 30 , 000 .
Table 8. The effect of dispersion in p on the in-control AARL of the p chart with n   = 30 , 000 .
( a )   n   = 30 , 000
p   = 0.0005
( b )   n   = 30 , 000
p ~ U 0.0004 ,   0.0006
( c )   n   = 30 , 000
p ~ U 0.00025 ,   0.00075
m AARLSDARLAARLSDARLAARLSDARL
20443.59472.94433.71448.74523.35494.67
50432.40403.23407.10394.53566.80458.88
100398.18324.28361.43301.76548.09397.36
1000313.8095.58290.3888.48450.77149.62
2000297.0861.10276.3668.14438.84106.54
3000290.8436.12268.4855.53434.3479.31
4000288.1525.19263.7648.42431.6565.04
5000285.8316.01258.4536.11429.3149.07
6000285.4611.96257.0732.07427.9837.20
8000284.636.60254.2723.19428.1723.08
10,000284.544.35253.4118.89427.7218.38
Table 9. The effect of dispersion in n and p on the in-control AARL of the p chart.
Table 9. The effect of dispersion in n and p on the in-control AARL of the p chart.
( a )   n ~ U 2400 ,   3600
p   = 0.005
( b )   n ~ U 2400 ,   3600
p ~ U 0.004 ,   0.006
( c )   n ~ U 2400 ,   3600
p ~ U 0.0025 ,   0.0075
m AARLSDARLAARLSDARLAARLSDARL
20330.38167.73154.3672.0926.408.17
50323.90123.73151.3051.9627.566.44
100312.6788.46151.4540.0027.684.77
1000303.1927.92146.1312.3327.881.74
2000302.5821.15146.129.7427.881.42
3000302.5516.79146.017.0427.861.10
4000302.5116.64146.016.4427.881.10
5000302.4914.32145.925.5827.890.93
6000302.3411.55145.805.0827.920.87
8000302.3010.53145.764.6427.920.82
10,000302.279.72145.734.4327.900.80
Table 10. The effect of dispersion in n and p on the in-control AARL of the p chart.
Table 10. The effect of dispersion in n and p on the in-control AARL of the p chart.
( a )   n ~ U 2400 ,   3600
p   = 0.005
( b )   n ~ U 2400 ,   3600
p ~ U 0.004 ,   0.006
( c )   n ~ U 2400 ,   3600
p ~ U 0.0025 ,   0.0075
m AARLSDARLAARLSDARLAARLSDARL
20471.28478.12429.25455.93503.63488.42
50429.90400.32408.95393.01568.92453.98
100395.38315.21373.32293.14551.73402.27
1000309.6778.91291.2172.93446.84127.65
2000304.4652.31285.8649.15440.3889.36
3000304.0443.16284.7239.19435.1771.80
4000303.2337.83283.1234.93430.6661.52
5000302.4534.38282.5229.67430.4657.64
6000301.0330.64282.4328.04429.9251.35
8000301.0226.24281.8522.92429.4344.11
10,000300.1024.17280.9921.08430.9241.88
Table 11. The effect of sample size n on the in-control AARL of the p chart.
Table 11. The effect of sample size n on the in-control AARL of the p chart.
p ~ U 0.25 ,   0.75 p ~ U 0.05 ,   0.15 p ~ U 0.03 ,   0.09 p ~ U 0.025 ,   0.075 p ~ U 0.015 ,   0.045
n AARLSDARLAARLSDARLAARLSDARLAARLSDARLAARLSDARL
504.820.2483.762.04109.6342.97127.8427.85143.2131.89
1002.460.0241.896.7777.8516.9770.471.16113.0726.67
3001.530.018.060.4221.371.5528.592.3745.6910.74
10001.230.012.310.033.680.094.600.159.890.70
15001.180.011.870.022.560.043.000.065.400.21
30001.120.001.490.011.760.011.910.022.620.04
15,0001.050.001.170.001.240.011.270.011.390.01
30,0001.030.001.120.001.150.001.180.001.250.01
Table 12. AARL and SDARL performance of the p chart with m =   1000 .
Table 12. AARL and SDARL performance of the p chart with m =   1000 .
(a)
p   = 0.005
(b)
p ~ U 0.004 ,   0.006
(c)
p ~ U 0.0025 ,   0.0075
(d)
p ~ U 0.0005 ,   0.0015
ShiftAARLSDARLAARLSDARLAARLSDARLAARLSDARL
0.0292.1719.08141.9411.2227.692.37107.2913.82
1.022.351.7217.400.888.710.6114.981.44
1.59.300.588.150.365.400.317.700.66
2.04.780.244.510.173.670.184.600.34
2.52.900.102.860.082.670.103.090.19
3.02.020.062.040.042.070.072.280.11
Table 13. AARL and SDARL performance of the p chart with m =   3000 .
Table 13. AARL and SDARL performance of the p chart with m =   3000 .
(a)
p   = 0.005
(b)
p ~ U 0.004 ,   0.006
(c)
p ~ U 0.0025 ,   0.0075
(d)
p ~ U 0.0005 ,   0.0015
ShiftAARLSDARLAARLSDARLAARLSDARLAARLSDARL
0.0293.3332.16271.7050.38441.0672.52110.709.54
1.022.521.7726.543.4256.257.4815.320.81
1.59.340.6111.551.2525.983.027.850.28
2.04.800.255.970.5213.711.344.680.16
2.52.910.113.580.248.060.673.140.09
3.02.030.062.420.125.190.362.310.06
Table 14. Data for Example 1.
Table 14. Data for Example 1.
Sample NumberPhaseSample Size
(Number of Inspected Units)
Number of Nonconforming UnitsSample Fraction Nonconforming
1I5400540.01000
2I5255620.01180
3I56751280.02256
4I52351220.02330
5I2355140.00594
6I4315580.01344
7I5155480.00931
8I5560600.01079
9I5720820.01434
10I2500680.02720
11I5095780.01531
12I2870360.01254
13I6000560.00933
14I3845600.01560
15I1995320.01604
16I50701240.02446
17I2995750.02504
18I4705380.00808
19I2480660.02661
20I65001480.02277
21II5255620.01180
22II2995720.02404
23II4705380.00808
24II2480360.01452
25II65002320.03569
Table 15. Data for Example 2.
Table 15. Data for Example 2.
Sample NumberPhaseSample Size
(Number of Inspected Units)
Number of Nonconforming UnitsSample Fraction Nonconforming
1I210,00014800.00705
2I268,50019330.00720
3I272,50019340.00710
4I274,70025540.00930
5I364,50025600.00702
6I279,00018410.00660
7I280,00025760.00920
8I276,25025130.00910
9I269,75027780.01030
10I374,35023580.00630
11I264,95022550.00851
12I276,80016510.00596
13I228,70018290.00800
14I277,55023860.00860
15I220,00014750.00670
16I272,50019340.00710
17I274,50025320.00922
18I364,50026250.00720
19I268,60023540.00876
20I364,50024600.00675
21I279,00018410.00660
22I280,00025760.00920
23I268,50018330.00683
24I269,75025810.00957
25I328,25024230.00738
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

Chen, P.-W.; Cheng, C.-S.; Wang, C.-W. A Study on the Laney p′ Control Chart with Parameters Estimated from Phase I Data: Performance Evaluation and Applications. Mathematics 2023, 11, 3411. https://doi.org/10.3390/math11153411

AMA Style

Chen P-W, Cheng C-S, Wang C-W. A Study on the Laney p′ Control Chart with Parameters Estimated from Phase I Data: Performance Evaluation and Applications. Mathematics. 2023; 11(15):3411. https://doi.org/10.3390/math11153411

Chicago/Turabian Style

Chen, Pei-Wen, Chuen-Sheng Cheng, and Ching-Wen Wang. 2023. "A Study on the Laney p′ Control Chart with Parameters Estimated from Phase I Data: Performance Evaluation and Applications" Mathematics 11, no. 15: 3411. https://doi.org/10.3390/math11153411

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