Research Article  Open Access
Yeneneh Tamirat, FuKwun Wang, YenChih Chen, "Sampling Schemes by Variables Inspection for the FirstOrder Autoregressive Model between Linear Profiles", Mathematical Problems in Engineering, vol. 2017, Article ID 3189412, 11 pages, 2017. https://doi.org/10.1155/2017/3189412
Sampling Schemes by Variables Inspection for the FirstOrder Autoregressive Model between Linear Profiles
Abstract
We present four new sampling schemes by variables inspection to deal with the firstorder autoregressive model between linear profiles. The first plan is based on exponentially weighted moving average (EWMA) and the rest of three plans are using the resubmitted sampling, repetitive group sampling (RGS), and multiple dependent state (MDS) sampling schemes. The nonlinear optimization problem is developed to find the number of profiles and the corresponding acceptance criteria, such that the producer’s and consumer’s risk are satisfied simultaneously. The efficiency of the proposed plans is compared with the conventional single sampling plan in terms of average sample number and the probability of acceptance. The result implies that all of the proposed sampling plans are superior to the single acceptance sampling plan by variables. In addition, the EWMA method appeared to be better than the others. The applications of proposed plans are shown with the help of industrial examples taken from calibration of an optical imaging system, and tire cornering stiffness test.
1. Introduction
Functional profiles express the relationship between one or more explanatory variables, the controllable inputs and a response variable which is the criticaltoquality characteristics. Thus, in statistical profile monitoring instead of individual or vector of quality characteristics we monitor the curve. Montgomery [1] defines a profile that is characterized by a functional relationship between a response variable and one or more independent variables. For example, profile monitoring has extensive application in the calibration of an optical imaging system; in tire cornering stiffness test the relationship between force and displacement can be modeled by linear profile. The fundamental concepts, methods, and issues regarding linear and nonlinear profiles monitoring methods can be found in [2]. A process yield index with the lower confidence bound for simple nonlinear and linear profiles can be found in [3] and [4], respectively. In the presence of autocorrelation between profiles [5] developed a process yield index with its approximate lower confidence bound. The author’s simulation study confirmed the method performs wellregarding bias and coverage rate.
Acceptance sampling plans lay down conditions for acceptance or rejection of the immediate lot inspected. Process capability is of utmost importance in acceptance sampling; in no event should the requirements of a sampling procedure exceed the producer’s process capability [6]. Moreover, in a highly competitive environment, sampling plans must be appropriately applied. In practice, the fraction of defectives a customer is willing to accept is getting very small, often measured in parts per million. That is, in order to reflect the actual lot quality, the sample size must be very large. In such circumstances, various authors [7–10] found variables sampling plans based on capability indices to be the most efficient methods of lot sentencing. In addition, variables sampling plans, based on the exponentially weighted moving average (EWMA) statistic, resubmitted sampling, repetitive group sampling (RGS), and multiple dependent state (MDS) sampling are proved to guarantee a reduction in sample size [11]. However, each sampling scheme should be applied cautiously to different scenarios. The scenarios are described as follows.
Scenario one: when the relationships between suppliers and customers have a longterm orientation, say ten or more previous lots are accepted, there would be an accumulated quality history. The EWMA sampling scheme takes into consideration past histories of the supplier to accept or reject the immediate lot inspected. The weights decline geometrically with the time of the observations. In quality control charts this EWMA statistic is known to be efficient to detect a small shift in the process [12–14]. Aslam et al. [15] applied the EWMA statistic in an acceptance sampling plan. For a process with linear profiles, [16] provided a single variable sampling plan based on EWMA model using the yield index . In the literature Aslam et al. [15] and Wang [16] suggested that a sampling plan based on EWMA yield index for a normal process and a process with linear profiles has more flexibility and economy than the single sampling plan based on yield index while providing the same protection to both suppliers and buyers. More information about EWMA sampling plan can be found in [17].
Scenario two: if either the producer or the consumer does not agree with the sample inspection result, resubmission might be necessary for reaching an appropriate decision. That is, if a lot is not accepted, then a resubmitted lot is tested and a decision is made irrespective of the previous test. In resubmitted sampling plans for a normal process, Liu et al. [19] propose variables sampling plan for twosided specification; [20] presented the operating procedure of the resubmitted sampling plan based on onesided capability indices. When the data can be considered as identically, independently, and normally distributed linear profiles, [21] proposed variables sampling plan for resubmitted lots based on the yield index . The author found that the resubmitted sampling plan has betteroperating characteristics curve (OC) than single sampling plan. Other research related to sampling inspection for resubmitted lots can be found in the literature [22].
Scenario three: in some situations testing is costly and destructive; a small sample is taken from the lot and analyzed. In the case of noncompliance, the second small sample is taken from the lot and analyzed. This process is repeated till the lot is accepted or rejected. The procedure is similar to sequential sampling. Balamurali and Jun [23] reported that the RGS plan is more efficient in terms of average sample number (ASN), cost, and time than single and double sampling plans. When the inspection is costly and destructive an RGS by variables based on yield index can be found in [24]. More information on RGS can be found in the literature [25, 26].
Scenario four: when the test result is marginal there might be a feeling the result is an inaccurate representation of the quality submitted. The inspection efficiency in such instances can be improved by utilizing the accumulated quality history of previous lots information. The MDS sampling plans are not only based on the current lot, but also on preceding lots [27]. With the absence of autocorrelation between linear profiles, [16] provided MDS sampling plan based on yield index . More information on MDS can be found in the literature [28].
The principal advantage of variables plans is the reduction in sample size. It must be remembered, however, that the superiority of the variables plan rests on the assumption of the underlying distribution of the measurements [6]. Whether the performance measure is single quality characteristic [7–10, 15] or profile data [16, 21] sampling schemes in the literature share two basic assumptions: the quality characteristic follows a normal distribution and data is independent. However, in some manufacturing processes, data is correlated or selfdependent. The effect of autocorrelations on the statistical properties of yield indices can be found in [5]. Accordingly, when dealing with sampling plans based on yield index the autocorrelation effect has to be taken into consideration. However, to the best of our knowledge, there is no work on sampling plan in the literature to deal with autocorrelation between profiles. In this study, we propose four new acceptance sampling plans based on the yield index for a firstorder autoregressive process. The first method is based on EWMA. The remaining three are based on resubmitted sampling, RGS, and MDS sampling, respectively. The rest of the paper is organized as follows. In the next section, yield index for autocorrelation between linear profiles is reviewed. Four sampling plans by variables are presented in Section 3. In the subsequent section, the performances of the proposed plans are compared and two illustrative examples are provided in Section 5. The former is based on the line widths of three photomask reference standards; the latter is an experimental example for tire cornering stiffness test on the ice road. Finally, we offer a conclusion and suggestions for future studies.
2. Yield Index for Autocorrelation between Linear Profiles
The profiles addressed in this paper are of the form of autocorrelation between linear profiles given aswhere n represents the number of levels and k is the number of profiles. For each of the k profiles, the response y is measured at the same n fixed locations called levels, where is the th level of the independent variable, denotes the firstorder autoregressive coefficient, denotes correlated random error, and are the intercept and slope of linear profiles, respectively, and .
Wang and Tamirat [5] derived the following process yield index for autocorrelation between linear profiles. This index is used to describe the relationship between manufacturing specifications and actual process yield.wherewhere and are the process mean and standard deviation of the response variable at the ith level of the independent variable, and are the upper and lower specification limits of the response variable at the ith level of the independent variable, respectively, , , , , is the cumulative distribution function of a standard random normal distribution, and is the inverse function of Ф. Furthermore, [5] derived the asymptotic normal distribution of index , which is given as follows:where where is the th lag autocorrelation and is the probability density function of a standard normal distribution.
3. Acceptance Sampling Plans
A sampling plan must satisfy two conditions: the probability of accepting a lot at the acceptable quality level (AQL) is greater than the producer’s confidence level and the probability of accepting a lot at the lot tolerance proportion defect (LTPD) is smaller than the consumer’s risk . Thus, for a particular sampling scheme its operating characteristic (OC) curve must pass through those two designated points (AQL, ) and (LTPD, ). EWMA, resubmitted, RGS, and MDS sampling schemes based on are presented and discussed in the following four subsections. The selection of a sampling scheme is based on purpose and information availability as discussed in the introduction part.
3.1. Sampling Plan Based on EWMA Model
The proposed sampling plan is given as follows.
Step 1. Choose the producer’s risk and consumer’s risk. Select the process yield requirements at two risks.
Step 2. Take a random sample of linear profiles k from the current lot and compute the process yield index .
Step 3. Collect process yield indices from previous lots and compute the EWMA statistic .
The EWMA model with yield index proposed by [29] is given bywhere is process yield index from current lot, is obtained from previous lots, and is smoothing constant and ranges from 0 to 1. The choice of the optimal EWMA parameter is based on minimizing the error sum of squares from previous lots, , where . Wang and Tamirat [5] apply the firstorder expansion of vvariate Taylor and central limit theorem to derive the asymptotic normal distribution of index ; see (4). The properties of EWMA can be obtained from the formula for the sum of a geometric series. The variance of the EWMA statistic is derived by . When the process is incontrol, , the variance of the EWMA statistic becomes [12]. In a steady state process, the mean and variance of are obtained as and , respectively. According to Montgomery [1], the EWMA can be viewed as a weighted average of all past and current observations; it is very insensitive to the normality assumption.
Step 4 (make a decision). Accept the lot if , where c_{a} is the critical value; otherwise reject it.
The lot acceptance probability, say , is established as
When , this sampling plan becomes a single sampling plan based on linear profiles. The ASN for this sampling plan is . The parameters of our proposed plan are determined through the following nonlinear optimization problem:
where, .
3.2. Resubmitted Sampling Plan
Resubmission is allowed times excluding the first inspection. The resubmitted sampling plan based on is described below.
Step 1. Choose the producer’s risk (α), the consumer’s risk (β), two process capability indices ( and ), and resubmission () times.
Step 2. Draw random linear profiles k from a submitted lot and compute the value.
Step 3 (make a decision). Accept the lot if , where is a critical value; otherwise resubmit a lot and go to Step 2. We can reject the lot if it is not accepted on the resubmission.
The lot acceptance probability, say , is established aswhereThe ASN for this resubmitted sampling plan is given as follows
The plan parameters of this scheme will be determined from the solution to the following nonlinear optimization problem:where .
3.3. Repetitive Group Sampling Plan
The RGS plan based on is described as follows.
Step 1. Choose the producer’s risk, the consumer’s risk, and two process capability indices ( and ).
Step 2. Draw random linear profiles k from a submitted lot and compute the value.
Step 3 (make a decision). Accept the lot if , and reject the lot if , where and are critical values; otherwise repeat Step 2.
The lot acceptance probability, say , is established asThe ASN for this RGS sampling plan is given as follows:When , this RGS plan becomes a traditional single sampling plan.
The parameters of the repetitive sampling plan can be determined through the following nonlinear optimization problem:where and .
3.4. Multiple Dependent State Sampling Plan
The MDS sampling plan based on is described as follows.
Step 1. Choose the producer’s risk, the consumer’s risk, two process capability indices ( and ), and the number of proceeding lots “.”
Step 2. Draw random linear profiles from a submitted lot and compute the value.
Step 3 (make a decision). Accept the lot if and reject the lot if , where and are critical values; accept the current lot if the proceeding lots were accepted under the condition of ; otherwise repeat Step 2.
The lot acceptance probability, say , is established aswhere and .
The average sample number for the MDS plan is established asWhen , this MDS plan becomes the RGS plan.
The parameters of the MDS sampling plan can be determined through the following nonlinear optimization problem:where , , , and are obtained by (19).
4. Comparison Study
In this section, the average sample numbers and probability of acceptance of the proposed plans are compared. Computations are carried out using computer programs written in the R language [30]. The programs are available from the authors. For instance, in order to determine the plan parameters for different combinations of using the sampling plan based on EWMA model, we used a simple algorithm to search the space of possible solutions. It consists of the following three steps.
Step 1. Assuming k ranges from 2 to 300 and is generated from a uniform distribution between and , 1,000 combinations of and were randomly generated.
Step 2. Evaluate 1,000 combinations of that must satisfy the two constraints (9) with the objective of minimizing the number of profiles.
Step 3. Repeat Steps 1 and 2 for 100 times to determine the optimal parameters.
As shown in Table 1 we consider some different combinations of and at a specified producer’s risk and customer’s risk . The probability of acceptance and the average number of profiles were calculated for the specified plans. From Table 1 all the proposed methods outperform the existing single variable scheme. The EWMA requires the smallest average number of profiles. Furthermore, it can be observed that the difference in ASN values between the variables single sampling plan and the proposed sampling plans increases when the gap between and narrows. This implies that, with the same protection level, the proposed plans will minimize the time and cost of inspection especially for cases when the process fraction of defectives required is very low.
 
is probability of accepting at and is probability of accepting at . 
The slope of the OC curves in Figures 1–4, respectively, displays the discriminatory power of the sampling scheme between good and bad lots. The steeper the slope of the OC curve, the greater the discriminatory power. Figures 1–4 also illustrate how different values of an input variable impact the sensitive the OC curves of each sampling scheme. Figure 1 indicates that the smaller the value of , the better the power. For example, if the quality level desired by the consumer is , the probability of acceptance at values of 0.25, 0.5, 0.75, and 1.0 became 1.0, 0.99, 0.96, and 0.91, respectively. That is, the associated producer’s risk turned out to be 0, 0.01, 0.04, and 0.09, respectively. Figures 2 and 3 indicate that increasing the number of resubmissions and repetitions improves the power. For instance, at a given = 1.33 and = 1.00, increasing the allowable resubmissions from two to four improves the probability of acceptance of a good lot from 0.96 to 0.99. For the MDS the effect of a number of preceding lots on the power is very marginal. From these comparisons, we may conclude that all of the proposed sampling plans are superior to the variable single acceptance sampling plan in terms of power. Figure 5 compared the efficiency of the proposed sampling plans with an existing single sampling plan at = 1.50, α = 0.05, and β = 0.10 for several values of = 1.0(0.05)1.4. The result indicates that the EWMA sampling plan is more efficient in terms of average sample number than other four schemes, MDS is second best, and resubmitted and RGS are very close to each other.
As stated in the introduction each scheme is designed for specific purpose and information availability. For example, EWMA is useful when ten or more previous lots are accepted. However, in some cases, data might be available for two or more methods, for such instances based on Table 1 and Figure 5 the EWMA technique appeared to be more efficient in terms of cost than the other methods. That is, it provides the smallest ASN.
5. Illustrative Examples
Example 1. We considered the calibration of an optical imaging system in [18]. The process is in statistical control, the functional profiles from phase I analysis resulted in an AR1 model, and the relationship between the response variable and the explanatory variable is modeled as , with a residual standard deviation of 0.07. Table 2 shows specification limits of the response variable. Suppose that the values of and are set to 1.50 and 1.25 and the two levels of risk are set (0.05, 0.1). Table 3 depicts sample plans that are equivalent, since they all provide the same protection to both suppliers and customers. The EWMA technique provides the smallest ASN. The plan parameters are found to be = 1.4057 and . Furthermore, the yield index of 10 previously accepted lots are assumed to be (1.4862, 1.5112, 1.5134, 1.5034, 1.4991, 1.5002, 1.4715, 1.4908, 1.4729, and 1.4862) the optimal value of is obtained as 0.1273. We collected 26 profiles and the immediate lot inspected resulted in = 1.2613. The corresponding value became 1.4637 since is greater than 1.4057 according to the EWMA scheme; the current lot can be deemed as accepted.


Example 2. This example is a real case study from the automobile parts industry. The goal was to determine whether the experimental setup is acceptable or not. The procedure is used to simulate tire cornering stiffness on the ice road. The lab measurements were done on the plunger tester equipped with the servohydraulic system. Type of tire tested is regular production, commercially available tire (255/80R17.5) with snow tread. Tests were conducted by pull force developed by 20 kg counterweights; each time add a 20 kg counterweight and record the corresponding tire displacement. The maximum load is 100 kg; that is, the measurement is taken at 20, 40, 60, 80, and 100 kg, respectively. Table 4 shows the upper and lower specification limits at each level of the independent variable. The setup is shown in Figure 6, from left to right total experiment setup, counterweights, and closeup view of contact between tire and snow in a container.

Regression analysis was used to fit the relationship between force and displacement; we found that the relationship between the response variable and the explanatory variable can be modeled as , , with a residual standard deviation of 0.06, where is the ith level of the independent variable (). Suppose that the values of and are set to 1.50 and 1.33 and the two levels of risk are set . Five equivalent sampling plans are provided in Table 5. However, since we do not have prior test information it is not possible to apply EWMA and MDS. The resubmitted scheme is more appropriate than repetitive; that is, if a test is not accepted, make an improvement on the setup and the next setup is tested and a decision is made irrespective of the previous test. After the second test we found which is greater than the critical value, ; it is possible to conclude the setup is acceptable.

6. Conclusion
In this paper, we developed four alternative acceptance sampling plans based on the process yield index to deal with lot sentencing for autocorrelation between linear profiles. The proposed sampling plans reduce significantly the average sample number as compared with the traditional single sampling plan. Under specific conditions, the proposed sampling plans become equal and reduced to a single sampling plan by variables. Since the EWMA method considers the quality history of the previous lots the number of profiles required is smaller than the other three sampling schemes. RGS and MDS schemes are equivalent; the advantage of MDS is it takes into consideration past history, which results in a reduction of cost and time of retesting. Resubmitted and RGS schemes can be used to solve disputed inspection procedures.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
This work was supported by the Ministry of Science and Technology in Taiwan under the Grant MOST 1052218E468003.
References
 D. C. Montgomery, Introduction to statistical quality control, vol. 1, Wiley, New York, NY, USA, 6th edition, 2003.
 R. Noorossana, A. Saghaei, and A. Amiri, Statistical analysis of profile monitoring, Wiley, New York, NY, USA, 2011.
 F.K. Wang and Y.C. Guo, “Measuring process yield for nonlinear profiles,” Quality and Reliability Engineering International, vol. 30, no. 8, pp. 1333–1339, 2014. View at: Publisher Site  Google Scholar
 F.K. Wang, “A process yield for simple linear profiles,” Quality Engineering, vol. 26, no. 3, pp. 311–318, 2014. View at: Publisher Site  Google Scholar
 F.K. Wang and Y. Tamirat, “Process yield analysis for autocorrelation between linear profiles,” Computers & Industrial Engineering, vol. 71, no. 1, pp. 50–56, 2014. View at: Publisher Site  Google Scholar
 E. G. Schilling and D. V. Neubauer, Acceptance sampling in quality control, Statistics: Textbooks and Monographs, CRC Press, Boca Raton, FL, USA, 2nd edition, 2009. View at: Publisher Site  MathSciNet
 C.W. Wu and S.W. Liu, “Developing a sampling plan by variables inspection for controlling lot fraction of defectives,” Applied Mathematical Modelling: Simulation and Computation for Engineering and Environmental Systems, vol. 38, no. 910, pp. 2303–2310, 2014. View at: Publisher Site  Google Scholar  MathSciNet
 L.C. Sheu, C.H. Yeh, C.H. Yen, and C.H. Chang, “Developing acceptance sampling plans based on incapability index Cpp,” Applied Mathematics & Information Sciences, vol. 8, no. 5, pp. 2509–2514, 2014. View at: Publisher Site  Google Scholar
 N. Kurniati, R.H. Yeh, and C.W. Wu, “Designing a variables twoplan sampling system of type TNTVSS(nT, nN; K) for controlling process fraction nonconforming with unilateral specification limit,” International Journal of Production Research, vol. 53, no. 7, pp. 2011–2025, 2015. View at: Publisher Site  Google Scholar
 S.W. Liu and C.W. Wu, “A quick switching sampling system by variables for controlling lot fraction nonconforming,” International Journal of Production Research, vol. 54, no. 6, pp. 1839–1849, 2016. View at: Publisher Site  Google Scholar
 M. Aslam and F.K. Wang, “Acceptance sampling plans for linear profiles with onesided specifications,” Journal of Statistical Computation and Simulation, vol. 87, no. 4, pp. 806–816, 2017. View at: Publisher Site  Google Scholar  MathSciNet
 P. Čisar and S. M. Čisar, “Optimization methods of EWMA statistics,” Acta Polytechnica Hungarica, vol. 8, no. 5, pp. 73–87, 2011. View at: Google Scholar
 J. S. Hunter, “The exponentially weighted moving average,” Journal of Quality Technology, vol. 18, no. 4, pp. 203–210, 1986. View at: Google Scholar
 J. M. Lucas and M. S. Saccucci, “Exponentially weighted moving average control schemes: properties and enhancements,” Technometrics, vol. 32, no. 1, pp. 1–29, 1990. View at: Publisher Site  Google Scholar  MathSciNet
 M. Aslam, M. Azam, and C.H. Jun, “A new lot inspection procedure based on exponentially weighted moving average,” International Journal of Systems Science, vol. 46, no. 8, pp. 1392–1400, 2015. View at: Publisher Site  Google Scholar
 F.K. Wang, “A Single Sampling Plan Based on Exponentially Weighted Moving Average Model for Linear Profiles,” Quality and Reliability Engineering International, vol. 32, no. 5, pp. 1795–1805, 2016. View at: Publisher Site  Google Scholar
 M. Aslam, M. Azam, and C.H. Jun, “Improved acceptance sampling plan based on EWMA statistic,” Sequential Analysis. Design Methods & Applications, vol. 34, no. 3, pp. 406–422, 2015. View at: Publisher Site  Google Scholar  MathSciNet
 C. Croarkin and P. Tobias, NIST/SEMATECH ehandbook of statistical methods, U.S., Department of Commerce, 2012.
 S.W. Liu, S.W. Lin, and C.W. Wu, “A resubmitted sampling scheme by variables inspection for controlling lot fraction nonconforming,” International Journal of Production Research, vol. 52, no. 12, pp. 3744–3754, 2014. View at: Publisher Site  Google Scholar
 N. Kurniati, R.H. Yeh, and C.W. Wu, “A sampling scheme for resubmitted lots based on onesided capability indices,” Quality Technology and Quantitative Management, vol. 12, no. 4, pp. 501–515, 2016. View at: Publisher Site  Google Scholar
 F.K. Wang, “Variables sampling plan for resubmitted lots in a process with linear profiles,” Quality and Reliability Engineering International, vol. 32, no. 3, pp. 1029–1040, 2016. View at: Publisher Site  Google Scholar
 K. Govindaraju and S. Ganesalingam, “Sampling inspection for resubmitted lots,” Communications in Statistics B: Simulation and Computation, vol. 26, no. 3, pp. 1163–1176, 1997. View at: Publisher Site  Google Scholar
 S. Balamurali and C.H. Jun, “Repetitive group sampling procedure for variables inspection,” Journal of Applied Statistics, vol. 33, no. 3, pp. 327–338, 2006. View at: Publisher Site  Google Scholar  MathSciNet
 C.W. Wu, M. Aslam, J. C. Chen, and C.H. Jun, “A repetitive group sampling plan by variables inspection for product acceptance determination,” European Journal of Industrial Engineering, vol. 9, no. 3, pp. 308–326, 2015. View at: Publisher Site  Google Scholar
 C.W. Wu, T.H. Wu, and T. Chen, “Developing a variables repetitive group sampling scheme by considering process yield and quality loss,” International Journal of Production Research, vol. 53, no. 7, pp. 2239–2251, 2015. View at: Publisher Site  Google Scholar
 A. Yan and S. Liu, “Designing a repetitive group sampling plan for Weibull distributed processes,” Mathematical Problems in Engineering, Article ID 5862071, Art. ID 5862071, 10 pages, 2016. View at: Publisher Site  Google Scholar  MathSciNet
 S. Balamurali and C.H. Jun, “Multiple dependent state sampling plans for lot acceptance based on measurement data,” European Journal of Operational Research, vol. 180, no. 3, pp. 1221–1230, 2007. View at: Publisher Site  Google Scholar
 M. Aslam, F. Wang, N. Khan, and C. Jun, “A multiple dependent state repetitive sampling plan for linear profiles,” Journal of the Operational Research Society. View at: Publisher Site  Google Scholar
 Y. Tamirat and F.K. Wang, “Sampling Plan based on the Exponentially Weighted Moving Average Yield Index for Autocorrelation within Linear Profiles,” Quality and Reliability Engineering International, vol. 32, no. 5, pp. 1757–1768, 2016. View at: Publisher Site  Google Scholar
 R Core Team, R: A language and environment for statistical computing, R Foundation for Statistical Computing, Vienna, Austria, 2016.
Copyright
Copyright © 2017 Yeneneh Tamirat et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.