Print version ISSN 0101-7438
Pesqui. Oper. vol.32 no.1 Rio de Janeiro Jan./Apr. 2012 Epub Apr 03, 2012
Linda Lee HoI,*; Roberto QuininoII
IUniversidade de São Paulo, Departamento de Engenharia de Produção Escola Politecnica. E-mail: firstname.lastname@example.org
IIUniversidade Federal de Minas Gerais, Departamento de Estatística ICEX. E-mail: email@example.com
In this paper, a procedure for the on-line process control of variables is proposed. This procedure consists of inspecting the m-th item from every m produced items and deciding, at each inspection, whether the process is out-of-control. Two sets of limits, warning (µ0 ± W) and control (µ0 ± C), are used. If the value of the monitored statistic falls beyond the control limits or if a sequence of h observations falls between the warning limits and the control limits, the production is stopped for adjustment; otherwise, production goes on. The properties of an ergodic Markov chain are used to obtain an expression for the average cost per item. The parameters (the sampling interval m, the widths of the warning, the control limits W and C(W < C), and the sequence length (h) are optimized by minimizing the cost function. A numerical example illustrates the proposed procedure.
Keywords: on-line process control, quality control by variables, Markov chain, warning limit, sequence of inspections.
Taguchi et al. (1985, 1989) proposed an economical procedure to monitor on-line process control with attributes as well as variables. This inspection system is automatic and allows for the sampling of only a single item at a time. In general, the proposed system can be implemented in high-speed electronics manufacturing facilities, where the testing equipment is connected to a central computer where data are taken automatically as the items that are produced are tested.
The procedure of on-line control for use in monitoring processes has been studied by many authors, such as Adams & Woodall (1989), Nayebpour & Woodall (1993), Srivastava & Wu (1991, 1995), Gong & Tang (1997), Borges et al. (2001), Wang & Yue (2001), Dasgupta (2003), Trindade et al. (2007), Dasgupta & Mandal (2008), Ho et al. (2007), Trindade et al. (2007) and Quinino et al. (2010). Ho & Trindade (2009) presented a revision related to this topic and proposed an economical X control chart model for short-run production.
In Taguchi's approach [Taguchi et al. (1989, 2004)], initially the process mean is µ = µ0 (state I). At a random time, the process mean may shift from µ = µ0 to µ = µ1 >µ0 (state II). To control the quality of the process, after every m items have been produced, the m-th item is examined. If the value of the quality characteristic does not satisfy the control limit, then production is stopped for adjustment. However, Taguchi et al. (1989; 2004) did not assume explicitly a probability function for the shift of the process mean from µ = µ0 to µ = µ1 > µ0, and many simplifications and approximations were used to calculate the average cost per produced item and the optimum values of the sampling interval m (hereafter denoted as mº) andthe width of the control limit C.
Adams & Woodall (1989) argued the inadequacy of Taguchi's model and presented alternative procedures for random walk models. The optimum sampling interval and control limits obtained by direct search are presented in tables. Srivastava & Wu (1991, 1995) presented a second order approximation for Taguchi's loss function and therefore presented closed expressions for the optimum sampling interval and control limits were also yielded.
Traditionally, three horizontal lines are used to draw the standard (Shewhart) control charts: a central line (representing the average value of the quality characteristic at an in-control state) and two other horizontal lines, called the control limits (upper and lower). However, some analysts may suggest two sets of limits. The outer limits are the usual action limits and are called control limits (CL) (points that plot outside of this limit lead to a stoppage for an assignable cause) and the inner limits are called warning limits (WL). Points that fall between the warning limits and the control limits may cause suspicion that the process may not be operating properly. One may also suggest for an additional rule, similar to the additional rules described in the Western Electric Handbook (1956) to decide if the production process is in-control or not. Such additional rules are usually based on the sequence of the results of last inspections.
The purpose of this paper is to present a procedure for on-line process control by adding alternative criteria that lead to an interruption of the production process not included in the previous listed papers. As in Taguchi's approach, the m-th item is inspected at every m produced items. If the value of the monitored statistic falls beyond the control limits (red zone RZ) or a sequence of h observations falls between the warning limits and the control limits (yellow zone YZ), then the production is stopped for adjustment; otherwise, (green zone GZ) production goes on. The parameters h, C, W and m are determined such that they minimize the average cost per item of the controlled system. Figure 1 represents schematically the new proposal.
This paper expands the procedure proposed by Ho et al. (2007). Their contribution can be viewed as a specific case of the new proposal, when h = 1 and W = C. The generalization proposed in the current paper, h > 1 and W < C, does not constitute a simple methodological expansion of Ho et al.'s proposal. A new development, based on the Markov chain, is required to facilitate the advancements and to provide better understanding for the reader. Additionally, restrictions related to the minimum value of the in-control average run length (ARL0) and a maximum value of the average run length during an out-of-control condition (ARLε) are included in the optimization process. A comparison of the results from the current proposal with the results obtained by Ho et al. (2007) is successful. This paper is organized as follows: Section 2 describes the probabilistic inspection model. In Section 3, an economic model is developed to determine the optimum designs. In Section 4, a numerical example illustrates the proposed model. The conclusions are presented in Section 5.
2 PROBABILISTIC MODEL
Consider a process that produces items and a quality characteristic of interest (X) that can be described by a Normal distribution with mean µ and standard deviation σ. The process starts in-control, in other words, its process mean is µ = µ0 (State I), and the process may shift to µ = µ1, µ1 ≠ µ0. The duration of the process in-control condition is usually modeled by an exponential distribution for a continuous case. The geometric distribution behaves similarly to the exponential distribution, but it is typically used for the discrete case in which the duration is measured by the number of units that are produced before the shift. Following previous papers [Nayebpour & Woodall (1993), Nandi & Sreehari (1999), Jiang & Tsui (2000), Borges et al. (2001), Ho et al. (2007), Trindade et al. (2007), Dasgupta (2008) and Ding & Gong (2008)], this study relies on a geometric distribution with parameter π ,0 < π < 1, to describe the random shift from State I to State II. One assumes that the probability of shifts from State II to State I is zero without any intervention in the process.
Both the geometric and exponential distributions are memoryless, which facilitates mathematical analysis. However, these distributions are useful for other reasons than their mathematical facility. Many researchers have recently used exponential or geometric distributions to describe shifts from in-control to out-of-control states. For example, Wang & Sheu (2003), Zhang et al. (2008), Serel & Moskowistz (2008), and Lim & Cho (2009) employed exponential distributions, but Ho et al. (2007), Trindade et al. (2007), Dasgupta (2008) and Ding & Gong (2008) used geometric distributions. Such distributions not only facilitate the development of mathematical models but also allow application to real problems, as mentioned in the cited papers.
Three horizontal lines are used to build standard (Shewhart) control charts, one drawn at µ = µ0 and two others at µ0 ± C (upper and lower control limits (CL)). However, one may suggest two sets of limits; the outer limits are the usual action limits (CL), and the inner limits are drawn at µ0 ± W , W < C. To illustrate these limits (see Fig. 2), one may use different colors to identify the different zones. Points that plot outside of the control limits lead to a stoppage for an assignable cause (red zone); points that fall between the warning limits and the control limits may reveal that the process may not be operating adequately (yellow zone); otherwise, the process must operate adequately if the points fall between the warning limits (green zone).
Let Xi , i = 1, 2, 3,..., ∞, denote the observed values of some characteristic of interest at the i-th inspection. The system of inspection may be described by a Markov chain with a set of states (s, k). The first index indicates in which state [State I: in-control (s = 0) or State II: out-of-control (s > 0)] the process is at the moment of inspection and assumes three values: s = 0, the inspected item is produced at µ = µ0; s = 1, the inspected item is produced at µ = µ1, but the process mean shifts from µ = µ0 to µ = µ1 in the current inspection (thus, a portion of the items is produced at µ = µ0, and at least the inspected item is produced at µ = µ1); s = 2, the process mean has shifted to µ = µ1, and all of the items in the current inspection are produced at µ = µ1. The second index k is related to the result of the inspection:
k = 1, the observed value of the current inspection is in the red zone (xi > µ0 + C or xi <µ0 C); the next inspection is only after a production of m items, and the process is stopped for adjustment;
k = 0, the observed value is in the green zone, (µ0 W < Xi < µ0 + W) and the process is not adjusted;
k = 1, the observed value is in the yellow zone (µ0 + W < Xi < µ0 + C or µ0 C < Xi < µ0 W), but the value of the previous inspection was not in the yellow zone and the process is not adjusted;
k = 2, the observed value is in the yellow zone (µ0 + W < Xi < µ0 + C or µ0 C < Xi < µ0 W), like the previous inspection and the process is not adjusted;
k = 3, the observed value is in the yellow zone (µ0 + W < Xi < µ0 + C or µ0 C < Xi < µ0 W), like the two previous inspections and the process is not adjusted;
k = h, the observed value is in the yellow zone (µ0 + W < Xi < µ0 + C or µ0 C < Xi < µ0 W), like the (h 1) previous inspections; the process is adjusted.
Note that there is an adjustment only if k =1 or k = h. A total of 3(h + 2) states of chain (s, k) is used to describe the inspection process. The second index assumes integer values in the interval 1 < k < h and k = t > 0, which indicates a sequence of t inspected items for which the observed values are in the yellow zone. Before the probabilistic model, let us introduce the following notations:
G0 = P(µ0 W < Xi < µ0 + W|µ = µ0)
G1 = P(µ0 W < Xi < µ0 + W|µ = µ1)
Y0 = P(µ0 + W < Xi < µ0 + C|µ = µ0) + P(µ0 C < Xi < µ0 W|µ = µ0)
Y1 = P(µ0 + W < Xi < µ0 + C|µ = µ1) + P(µ0 C < Xi < µ0 W|µ = µ1)
R0 = P(Xi > µ0 + C|µ = µ0) + P(Xi < µ0 C|µ = µ0)
R1 = P(Xi > µ0 + C|µ = µ1) + P(Xi < µ0 C|µ = µ1)
The transition probabilities from state (s, k) at the moment of inspection i to state (s*, k*) at the moment of inspection (i + 1) are elements of the matrix P
where Ass* denotes the matrix of transitory probabilities p(s, k)(s*,k*); s, s*= 0, 1, 2; k, k* = 1, 0, ..., h. The elements of A00 are the transitory probabilities from states (0, k) to states (0, k*). That is, the process is in-control at the i-th inspection, and it also stays in-control at the (i + 1)-th inspection. The non-null elements of the matrix A00 are expressed in (2).
The elements of A01 are the transition probabilities from states (0, k) to states (1, k*). That is, the process is in control at the i-th inspection, but the parameter shifts from µ = µ0 to µ = µ1. Thus, some items are produced at µ = µ0, and at least the inspected item is produced at µ = µ1. The non-null elements of the matrix A01 are expressed in (3), as follows:
for 0 < k < h, with qm = (1 π)m.
The matrix A02 is a null matrix; A10 = A20 and their non-null elements are the transitory probabilities with k =1 or h. That is, the process is adjusted in the previous inspection, and it restarts in-control in the current inspection employing the sampling interval m:
Similarly, the non-null elements of A11 = A21 are also transition probabilities with k = 1 or h. However, after the adjustment, the process restarts in-control, but the parameter µ shifts in the current inspection. Thus, the non-null elements of these matrices are
Next, we need the matrices A12 and A22, which are the transition probabilities from states (1, k) to states (2, k*) and from states (2, k) to states (2, k*), respectively. In these cases, the parameter shifted in the previous inspections, and the process has not been adjusted. A12 = A22 and the non-null probabilities of these matrices are
for 0 < k < h. Note that the Markov chain built in this way incorporates information about the process (whether it is in-control or out-of-control), the length of the sequence of observations in the yellow zone and also the length of sampling intervals m in the matrix P. Matrix P is ergodic recurrent and consequently Δ = limu → ∞ Pu exists and does not depend on the probability of the initial states of the process. Thus, all rows of Δ are equal. By denoting the first row of the matrix Δ by δ = (δ(0,1), . . . , δ(0,h), δ(1,1), . . . , δ(1,h), δ(2,1), . . . , δ(2,h)) and solving the system of linear equations δ = δP subject to the restriction δ(s, k) = 1, δ(s, k) can be determined. This solution can be interpreted as the proportion of inspections in state (s, k), s = 0, 1, 2; k =1, 0,..., h over a large number of inspections.
3 COST FUNCTION
To obtain the cost function, more assumptions are needed. Once one decides to make an adjustment, the stoppage of the process is instantaneous. After the adjustment, the process restarts at State I (µ = µ0). In this paper, the costs follow a structure similar to the usual economic designs, where ci is the cost for a single inspection and ca is the cost for adjustment. In this paper, one item is stated as non-conforming if the observed value of the quality characteristic is beyond the specification limits (µ0 ± LE). Thus, the next costs related to non-conforming and discarded items are also included in the model: cn is the cost for delivering a non-conforming item (this item is sent to the customer or to the following production stage), and cD is the cost to discard the examined item. For each cycle of inspection, the (m 1) items are sent to the customer or to the following stages of production. It is assumed that all inspected items are discarded or rectified in the posterior stages once the inspection is partially destructive.
For a sufficiently large number of inspections,
δ = (δ(0,1), . . . , δ(0,h), δ(1,1), . . . , δ(1,h), δ(2,1), . . . , δ(2,h))
is the vector of the probability of the states of the Markov chain. Let C be the random variables related to the cost at each cycle of inspection. We assume discrete values related to the states of the Markov chain. The cost of the state (s, k), s = 0, 1, 2; k =1, 0,..., h may be written as
ci , cD, respectively, are the costs to inspect and discard a single item, constants for all states (s, k), s = 0, 1, 2; k =1, 0,..., h;
a(s, k) is the cost to adjust the process;
n(s, k) is the cost to send non-conforming items to the customer or to later stages of the process.
In the next paragraphs, the different costs are detailed.
- Cost to adjust the process:
The cost for adjustment is included for states (s, k), s = 0, 1, 2; k =1, h. Thus,
a(s, k) = ca; s = 0, 1, 2; k =1, h.
- The cost to send non-conforming items to the customer or to later stages of the process:
For the states (0, k), k =1, 0,..., h, all items are produced at µ = µ0; thus,
n(0, k) = cn p1(m 1); k =1, 0,..., h,
with p1 = 1 P(µ0 LE < Xi < µ0 + LE|µ = µ0). Similarly for the states (2, k),when all items are produced at µ = µ1, we have
n(2, k) = cn p2(m 1);
for k =1, 0,..., h, with p2 = 1 P(µ0 LE < Xi < µ0 + LE|µ = µ1). For the states (1, k), k =1, 0,..., h, v< m items are produced at µ = µ0, and then, m v areproduced at µ = µ1. Taking into account all possibilities
for k =1, 0,..., h. For a large number of inspections and considering our setting as a renewal-reward process, the average cost per item CE(m) is the ratio of the expected cost per inspection cycle E(V ) by the quantity of items sent to the customer or to the following stages of production, expressed as
Montgomery (2001) claimed that the economical design of the control chart must be evaluated jointly with some statistical criteria. He recommended strongly that the optimization of the cost function should be subject to an adequate statistical restriction. In this study, a set of additional restrictions are considered: a minimum value for the in-control average run length (ARL0) isdenoted by ϕ0, and a maximum value for the average run length during an out-of-control condition (ARLε) is denoted by ϕ1. These values (ϕ0 and ϕ1) are fixed as a part of the responsibility for the control system of the productive process.
The problem consists of determining the values of
with m the sampling interval; the constants W and C are used to draw, respectively, the warning limit and the control limit, and h is the length of the sequence of values that are in the yellow zone.
The values of ARL0, ARLε can be calculated by employing the stationary vector
δ = (δ(0,1), δ(0,0), . . . , δ(0,h1), δ(0,h), δ(1,1), δ(1,0), . . . ,
δ(1,h1), δ(1,h), δ(2,1), δ(2,0), . . . , δ(2,h1), δ(2,h))
and can be expressed as:
4 NUMERICAL EXAMPLES AND DISCUSSIONS
To illustrate the proposed procedure, consider the example adapted from Taguchi et al. (1989), Taguchi et al. (2004) and Trindade et al. (2007).
A manufacturer of high-volume-integrated circuits wants to install a system to control the measurement of some dimension of interest. Historical data allows an estimation of the cost components as ci = $0.25, cn = $20, cD = $2.0, and ca = $900. The specification limits (LE) were fixed as ±1.5, and the shifts from process in-control [State I (µ0 = 0)] to out-of-control [State II (µ1 = 1)] can be described by a geometric distribution with the parameter π = 0.001. The standard deviation of the process is known and is equal to 0.5. Managers responsible for the production system adopt ϕ0 = 370 and ϕ1 = 5 (minimum and maximum values for ARL0 and ARLε, respectively). To calculate the optimum values, a program in MatLab was developed for this task. (Interested readers can request a copy of the program directly from the authors.)
Figure 3 shows the plots of the expected cost versus the sampling interval m for some cases chosen to illustrate the behavior of the optimum set. The optimal design is as follows: The sampling interval is m0 = 27. The length of the sequence of observations in the yellow zone is h0 = 3. The width of the warning limit is W 0 = 0.8, and the width of the control limit is C0 = 1.6; the results yield an average cost of $1.381. The current proposal is 4.4% less expensive when compared with the approach presented in Ho, Medeiros & Borges (2007). In that case, the sampling interval increases to m0 = 32 as the average cost also increases to $1.445 (with h = 1 and W = C = 1.4).
4.1 Sensitivity analysis
For the sensitivity analysis, each parameter is analyzed in a ±15% range of values, and other parameters are kept equal to the numerical values presented at the beginning of this section.
Table 1 summarizes the results. As expected, increases in the cost cn and the parameters µ1 and π yield a decrease in the values of the sampling interval m (and vice versa). The values of h, W and C are not sensible in the range of variation up to 15% in the parameters of the process and costs. Only increases in µ1 produce more significant changes. An increase of 6% in µ1 results in a decrease in h to 2 and increases in the values of W 0 and C0, respectively, to 2.0 and 3.4.
With these observations, it is possible to notice that the optimum policy is reasonably robust in the range of variation of approximately ±15% for the input parameters of the models. This evidence is important because the estimation of the parameters is not free of errors.
The optimum solution presented above is not influenced by the values of ϕ0 = 370 and ϕ1 = 5, which are chosen by managers responsible for the production process because the non-restricted solution of the expression (10) [that is, m0 = 27; h0 = 3; W 0 = 1.6, C0 = 3.2, $1.381] presents values of ARL0 = 408.17 and ARLε = 4.87, which satisfy the two restrictions (ϕ0 = 370 and ϕ1 = 5). However, the solution may be altered depending on the values chosen for ϕ0 and ϕ1. Table 2 shows the results for some values of ϕ0 and ϕ1. An increase in ϕ0 (keeping ϕ1 constant) tends to increase the optimum interval sampling m while the values of h, W and C remain unaltered. An increase in ϕ1 does not alter significantly the results; however, when both ϕ0 and ϕ1 increase simultaneously, a considerable increase is observed in the sampling interval m. The value of m changes from 27 to 193. As expected, the cost increases as more demanding values of ϕ0 and ϕ1 are fixed.
The results presented in Table 1 are very useful, but they do not help us to indicate easily which parameters have a greater impact on the cost function. Identification of such parameters is very important because the parameter values are not free of errors, and those parameters that produce more impact on the cost need more (special) attention. For that, we use a full factorial experiment with the cost as the response and the following factors: X1 the cost to inspect (ci), X2 the cost to send a non-conforming item (cn); X3 the cost to scrap an inspected item (cD), X4 the cost for adjustment (ca), X5 the process mean when the process is out-of-control (µ1) and X6 the probability for a shift in the process mean (π). Three levels for each factor are chosen; the value described in the beginning of this section is set as a reference level, and two other values are fixed at ±15% from the reference level. Such design results in 729 possible combinations, and a multiple regression model with dummy variables will be employed. For each factor, two dummy variables are defined, as follows:
for j = 1,..., 6. The choice of which level will be the first, second or third level will not alter the results. Estimates of the coefficients β0, βj1, βj2, j = 1, ..., 6 of the regression model are obtained by the well-known least squares method. The range of the estimates of the coefficients (known as utilities) for each factor was the metric employed to evaluate the impact of each factor on the cost function. For more details, see Green & Srinivasan (1990) and Malhotra (1999). For example, the estimate of the coefficient (utility) of the first level of the factor j is , for the second level and for the third level ( ), with and obtained by the least squares method. The impact of factor j is proportional to the range [ max(, , , )] when compared with the range of the estimates of the coefficients from the other factors. Using this method, the impact of each factor is obtained and summarized in Table 3. Note that the factors X4, the cost for adjustment (ca), and X6, the probability for a shift in the process mean (π), yield higher impacts, 31.3% and 35.1%, respectively. Thus, these factors require careful evaluation and estimation. Some descriptive statistics may help to understand these results. In Figure 4, boxplots of each factor are drawn. Observe the upward tendency of the medians of the factors pointed out as more important while the medians of the other factors are kept constant when the levels alter.
In earlier papers related to this subject, the optimum corrective policy designs for Taguchi's on-line quality monitoring procedure for variables are presented, and the sampling interval and the width of the control limits are fixed. In Ho, Medeiros & Borges (2007), the authors determine the optimal sampling interval m and the optimal control limit that minimizes the expected cost function.
In this paper, we extend the system of control by employing control zones. Three zones are stated: the green (µ0 W < Xi < µ0 + W), yellow (µ0 + W < Xi <µ0 + C or µ0 C <Xi < µ0 W) and red (X > µ0 + C or Xi < µ0 C). An observed value in the green zone results in the decision that the process is in-control and the production goes on; however, one decides in favor of an adjustment of the process if the observed value belongs to the red zone; moreover, a sequential verification takes place if an observation is in the yellow zone. If a sequence of h values in the yellow zone is observed, then an adjustment is decided on.
The aim is to find the values of m, h, W and C that minimize the average cost. According to the results presented in Section 4, the inspection control and policy presented in this paper may yield results that are at least equal to or more economical than the results obtained by Ho, Medeiros & Borges (2007).
Future research can be proposed. For example, the number of zones may increase. There is a tendency to have more automation in process control, so more complex alternatives may not produce excessive complexity for the user. The corrective adjustment based on h sequential observations in the yellow zone (µ0 + W < Xi <µ0 + C or µ0 C < Xi < µ0 W), may be modified. For example Rule A: only h sequential observations in the warning region above the centerline µ0 + W < Xi < µ0 + C; Rule B only h sequential observations in the warning region below the centerline. Fixed h and the warning limit W , Rule A (applied in this study) is more likely (higher probability) than Rule B when the process is out-of-control. Otherhand, Rule B is less unlikely (lower probability) than Rule A when the process is in-control. Preliminary simulation study conducted by the authors suggest that rules A and B present similar performance under an economical perspective once the optimum parameters h, C and W are different for each rule with similar average cost. However, a wider research may be realized to obtain better conclusions. Another possibility is to consider different sizes of sampling intervals. A longer sampling interval may be chosen for the next inspection if the current observation falls in the green zone or a shorter sampling interval may be chosen in the next inspection if the current point falls in the yellow zone. Other distributions for the process failure mechanism can be developed. A natural candidate is a discrete Weibull distribution [Nakagawa & Osaki (1975)] because it can take on a variety of shapes and mimic the behavior of other distributions. Compared with a geometric distribution, the discrete Weibull distribution is more flexible in that it includes models with an increasing or decreasing hazard rate.
The authors would like to thank the suggestions from the two anonymous referees which have greatly contributed to the improvement of this paper. They would like to acknowledge CNPq and Fapesp for their partial financial support of this research.
 ADAMS BM & WOODALL WH. 1989. An analysis of Taguchi's on-line process control procedure under a random-walk model. Technometrics, 31: 401413. [ Links ]
 BORGES W, HO LL & TURNES O. 2001. An analysis of Taguchi's on-line quality monitoring procedure for attributes with diagnosis errors. Applied Stochastic Models in Business and Industry, 17: 261276. [ Links ]
 DASGUPTA T & MANDAL A. 2008. Estimation of process parameters to determine the optimum dagnosis interval for control of defective items. Technometrics, 50: 167181. [ Links ]
 DING J & GONG L. 2008. The effect of testing equipment shift on optimal decisions in a repetitive testing process. European Journal of Operational Research, 186: 330350. [ Links ]
 GREEN PE & SRINIVASAN V. 1990. Conjoint analysis in marketing: new developments with implications for research and practice. Journal of Marketing, 54: 319. [ Links ]
 HO LL, MEDEIROS PG & BORGES W. 2007. An alternative model for on-line quality monitoring for variables. International Journal of Production Economics, 107: 202222. [ Links ]
 JIANG W & TSUI KL. 2000. An economic model for integrated APC and SPC control charts. IIE Transactions, 32: 505513. [ Links ]
 MALHOTRA NK. 1999. Marketing research: an applied orientation. Prentice-Hall. Inc. [ Links ]
 NANDI SN & SREEHARI M. 1999. Some improvements in Taguchi's economic method allowing continued quality deterioration in production process. Communications in Statistics Theory and Methods, 28(5): 11691181. [ Links ]
 NAKAGAWA T & OSAKI S. 1975. The Discrete Weibull Distribution. IEEE Transactions on Reliability, R-24 5: 300301. [ Links ]
 NAYEBPOUR MR & WOODALL WH. 1993. An analysis of Taguchi's on-line quality monitoring procedure for attributes. Technometrics, 35: 5360. [ Links ]
 SRIVASTAVA MS & WU Y. 1995. An improved version of Taguchi's on-line control procedure. Journal of Statistical Planning and Inference, 43: 133145. [ Links ]
 SRIVASTAVA MS & WU Y. 1991. A second order approximation on Taguchi's on-line procedures. Communications in Statistics Theory and Methods, 20: 21492168. [ Links ]
 TAGUCHI G. 1985. Quality Engineering in Japan. Communications in Statistics Theory and Methods, 14: 27852801. [ Links ]
 TAGUCHI G, CHOWDHURY S & WU Y. 2004. Taguchi's Quality Engineering Handbook. John Wiley & Sons, Inc. New Jersey. [ Links ]
 TAGUCHI G, ELSAYED EA & HSIANG T. 1989. Quality Engineering in Production Systems. McGraw-Hill, New York. [ Links ]
 TRINDADE A, HO LL & QUININO R. 2007. Monitoring process for attributes with quality deterioration and diagnosis errors. Applied Stochastic Models in Business and Industry, 23(4): 339358. [ Links ]
 WESTERN ELECTRIC COMPANY. 1956. Statistical Quality Control Handbook, Western Eletric Company, Inc., Indianapolis, Indiana. [ Links ]
Received March 11, 2010 / Accepted November 12, 2011
* Corresponding author