MONITORING BIVARIATE PROCESSES

versão impressa ISSN 0101-7438 / versão online ISSN 1678-5142 MONITORING BIVARIATE PROCESSES Marcela A. G. Machado* Antonio F. B. Costa Fernando A. E...
0 downloads 1 Views 181KB Size
versão impressa ISSN 0101-7438 / versão online ISSN 1678-5142

MONITORING BIVARIATE PROCESSES Marcela A. G. Machado* Antonio F. B. Costa Fernando A. E. Claro Production Department São Paulo State University (UNESP) Guaratinguetá – SP, Brazil [email protected] [email protected] [email protected] * Corresponding author / autor para quem as correspondências devem ser encaminhadas

Recebido em 12/2008; aceito em 06/2009 Received December 2008; accepted June 2009

Abstract The T 2 chart and the generalized variance |S| chart are the usual tools for monitoring the mean vector and the covariance matrix of multivariate processes. The main drawback of these charts is the difficulty to obtain and to interpret the values of their monitoring statistics. In this paper, we study control charts for monitoring bivariate processes that only requires the computation of sample means (the ZMAX chart) for monitoring the mean vector, sample variances (the VMAX chart) for monitoring the covariance matrix, or both sample means and sample variances (the MCMAX chart) in the case of the joint control of the mean vector and the covariance matrix.

Keywords: control charts; bivariate processes; mean vector; covariance matrix.

Resumo Os gráficos de T 2 e da variância amostral generalizada |S| são as ferramentas usualmente utilizadas no monitoramento do vetor de médias e da matriz de covariâncias de processos multivariados. A principal desvantagem desses gráficos é a dificuldade em obter e interpretar os valores de suas estatísticas de monitoramento. Neste artigo, estudam-se gráficos de controle para o monitoramento de processos bivariados que necessitam somente do cálculo de médias amostrais (gráfico ZMAX) para o monitoramento do vetor de médias, ou das variâncias amostrais (gráfico VMAX) para o monitoramento da matriz de covariâncias, ou então das médias e variâncias amostrais (gráfico MCMAX) para o caso do monitoramento conjunto do vetor de médias e da matriz de covariâncias.

Palavras-chave: gráficos de controle; processos bivariados; vetor de médias; matriz de covariâncias.

Pesquisa Operacional, v.29, n.3, p.547-562, Setembro a Dezembro de 2009

547

Machado, Costa & Claro – Monitoring bivariate processes

1. Introduction The control charts are often used to observe whether a process is in control or not. When there is only one quality characteristic to control, the Shewhart charts are usually applied to detect process shifts. The power of the Shewhart control charts lies in its ability to separate the assignable causes of variation from the uncontrollable or inherent causes of variation. Shewhart control charts are relatively easy to construct and to interpret. As a result, they are readily implemented in manufacturing environments. However, there are many situations in which it is necessary to control two or more related quality characteristics simultaneously. Hotelling (1947) provided the first solution to this problem by suggesting the use of the T 2 statistic for monitoring the mean vector of multivariate processes. Costa & Machado (2007) studied the properties of the synthetic T 2 chart with two-stage sampling. Machado & Costa (2008a) considered the use of simultaneous X charts as an alternative to the use of the T 2 chart. Costa & Machado (2008a) considered the use of the double sampling procedure with the chart proposed by Hotelling. The first multivariate control chart for monitoring the covariance matrix Σ was based on the charting statistic obtained from the generalized likelihood ratio test (Alt, 1985). For the case of two variables, Alt (1985) proposed the generalized variance statistic S to control the covariance matrix Σ . Control charts more efficient than the S chart have been proposed. Recently, Costa & Machado (2008) considered the VMAX statistic to control the covariance matrix of bivariate processes. The points plotted on the VMAX chart correspond to the maximum of the sample variances of the two quality characteristics. Machado et al. (2008) obtained the properties of the VMAX chart with double sampling and Machado & Costa (2008) considered the EWMA chart based on the VMAX statistic. Costa & Machado (2008) extended the VMAX chart for the multivariate case, where p quality characteristics are under control. The synthetic VMAX chart for the bivariate case was also studied by Machado et al. (2009). There are quite a few recent papers dealing with the joint control of the mean vector and the covariance matrix of multivariate processes. Khoo (2005) proposed a control chart based on the T 2 and S statistics for monitoring bivariate processes. The speed with which the chart signals changes in the mean vector and/or in the covariance matrix was obtained by simulation. The results are not compelling, once the proposed chart is slow in signaling out-of-control conditions. Chen et al. (2005) proposed a single EWMA chart to control both, the mean vector and the covariance matrix. Their chart is more efficient than the joint T 2 and S in signaling small changes in the process. Zhang & Chang (2008) proposed two EWMA charts based on individual observations that are not only fast in signaling but also very efficient in informing which parameter was affected by the assignable cause; if only the mean vector or only the covariance matrix or both. Machado & Costa (2008b) proposed the joint use of two charts based on the non-central chi-square statistic for monitoring the mean vector and the covariance matrix of bivariate processes, named as the joint NCS charts. In general, the joint NCS charts are faster than the combined T 2 and S charts in signaling out-of-control conditions. Once the proposed scheme signals, the user can immediately identify the out-of-control variable. The aim of this paper is to study a set of new control charts for monitoring the mean vector and/ or the covariance matrix of bivariate processes.

548

Pesquisa Operacional, v.29, n.3, p.547-562, Setembro a Dezembro de 2009

Machado, Costa & Claro – Monitoring bivariate processes

This paper is organized as follows: the next section presents the assumptions and measure of efficiency. In Section 3, we present the ZMAX chart as an alternative to the T 2 chart. In Section 4, we present the VMAX chart as an alternative to the S chart. In Section 5, we show the gain in speed with which the VMAX chart signals when the variable sample size, double sampling, synthetic or EWMA procedures are applied. In Section 6, we present the synthetic MCMAX chart as an alternative to the joint use of the T 2 and S charts. Section 7 brings an illustrative example. Conclusions are in Section 8. 2. Assumptions and measure of efficiency

μ = μ

Throughout this article, it is assumed that the control charts are employed to monitor a bivariate process whose quality characteristics of interest (say, X and Y) are normally distributed with the mean vector µ and the covariance matrix Σ . The process is considered to start with the mean vector and the covariance matrix on target ( 0 and Σ = Σ 0 ), μ

where

' 0

⎛ σ x2 = (µ x ; µ y ) and Σ 0 = ⎜ ⎜⎜ ⎝ σ xy

σ xy ⎞

⎟ . Without losing generality we assume that ⎟ σ ⎟⎠ 2 y

μ

μ

σ x = σ y = 1.0 . The T 2 chart and the ZMAX chart are used to control processes that are subject to assignable causes that change the mean vector from

0

to

' 1

= (µ x + c; µ y + d ) ,

μ

μ

where c and d measure the magnitude of the shift on each variable. The S chart and the VMAX chart are used to control processes that are subject to assignable causes that change ⎛ a2 abσ xy ⎞ the covariance matrix from Σ 0 to Σ1 = ⎜ ⎟ , where a and b measure the ⎜ abσ b 2 ⎟⎠ xy ⎝ increase in the variability of each variable. Finally, the joint T 2 and S charts and the synthetic MCMAX are used to control processes that are subject to assignable causes that change the mean vector from 0 to 1 and/or change the covariance matrix from Σ 0 to Σ1 without changing the correlation, that is, ρ = σ xy . The numerical results presented in this paper are for the case of positive correlation, believed to be prevalent in industrial data. The number of samples that are taken from the process until the control chart produces a signal measures the efficiency of a control chart. The expected number of samples taken before the chart signals is called the average run length (ARL). During the in-control period, ARL = 1/ α and is called ARL0. The risk α is the well known Type I error. The control limits of the charts considered in this article were adjusted to assure a false alarm risk α = 0.005. 3. Control charts for monitoring the mean vector

A chart based on the sample means X = ∑ i =1 X i / n and Y = ∑ i =1Yi / n is considered to n

n

μ

control , where X and Y are two quality characteristics and n is the size of the samples. As an alternative to the use of two mean charts we consider a single chart based on the

Pesquisa Operacional, v.29, n.3, p.547-562, Setembro a Dezembro de 2009

549

Machado, Costa & Claro – Monitoring bivariate processes

statistic ZMAX = max{ n | X − µ x |, n | Y -µ y |} , shortly ZMAX chart. A signal is given if ZMAX > CL , the control limit for the ZMAX chart. The power of the ZMAX chart is given by: PM = 1 −



CL − c n



CL − d n

− CL − c n − CL − d n

(1)

f (Z x , Z y )dZ x dZ y

where f ( Z x , Z y ) is a standardized bivariate normal distribution function with correlation ρ . Table 1 – ARL values for the T 2 chart and for the ZMAX chart.

ρ 0.0

0.3 2

shifts ZMAX T c d CL 3.023 10.597 0.0 0.0 200.0 200.0 0.0 0.5 117.38 115.55 0.0 1.0 41.64 41.92 0.0 1.5 15.10 15.78 0.0 2.0 6.44 6.88 0.0 2.5 3.31 3.55 0.0 3.0 2.03 2.16

0.5 2

ZMAX T 3.021 10.597 200.0 200.0 117.43 110.45 41.63 37.97 15.08 13.85 6.43 5.97 3.30 3.10 2.03 1.93

0.7 2

ZMAX 3.015 200.0 117.47 41.47 14.98 6.38 3.28 2.02

T 10.597 200.0 99.73 30.60 10.51 4.47 2.38 1.56

ZMAX T2 2.996 10.597 200.0 200.0 115.78 77.98 40.32 18.98 14.53 5.94 6.21 2.58 3.21 1.53 1.99 1.16

0.5 0.5 0.5 0.5 0.5 0.5

0.5 1.0 1.5 2.0 2.5 3.0

83.15 36.41 14.38 6.32 3.28 2.02

76.87 32.95 13.64 6.28 3.35 2.09

84.04 36.91 14.53 6.36 3.29 2.03

91.65 40.09 15.75 6.84 3.49 2.11

85.43 37.49 14.63 6.36 3.29 2.02

99.73 41.92 15.01 6.10 3.02 1.83

87.04 37.82 14.46 6.24 3.22 1.99

106.71 38.74 11.36 4.18 2.07 1.36

1.0 1.0 1.0 1.0 1.0

1.0 1.5 2.0 2.5 3.0

23.44 11.89 5.83 3.17 1.99

18.49 9.36 4.92 2.87 1.90

24.11 12.27 5.98 3.22 2.01

25.82 13.03 6.45 3.49 2.16

24.93 12.65 6.09 3.25 2.02

30.60 15.01 6.88 3.48 2.06

25.96 13.01 6.10 3.21 1.99

35.25 15.73 6.14 2.82 1.65

1.5 1.5 1.5 1.5

1.5 2.0 2.5 3.0

8.09 4.83 2.90 1.91

5.76 3.55 2.32 1.67

8.50 5.08 3.02 1.96

8.53 5.09 3.10 2.05

8.91 5.30 3.10 1.99

10.51 6.10 3.48 2.16

9.42 5.52 3.14 1.98

12.58 6.83 3.45 1.97

2.0 2.0 2.0 2.5 2.0 3.0

3.54 2.45 1.76

2.51 1.85 1.45

3.79 2.62 1.85

3.63 2.53 1.83

4.01 2.75 1.90

4.47 3.02 2.06

4.26 2.88 1.94

5.39 3.45 2.15

2.5 2.5 2.5 3.0

1.96 1.55

1.50 1.27

2.11 1.66

1.99 1.59

2.24 1.74

2.38 1.83

2.38 1.82

2.82 2.07

3.0 3.0

1.35

1.15

1.44

1.37

1.51

1.56

1.59

1.78

550

Pesquisa Operacional, v.29, n.3, p.547-562, Setembro a Dezembro de 2009

Machado, Costa & Claro – Monitoring bivariate processes

Table 1 brings the ARLs of the ZMAX and T 2 charts. In terms of performance, the T 2 chart is highly affected by the correlation ρ , while ρ has a minor influence on the ZMAX chart. When the correlation is moderate (0.3-0.5) the ZMAX chart has a better overall performance (see ARL values in bold). Table 1 was built fixing n=5 and ARL0 = 200.0. The study with other n values led to the same conclusions.

4. Monitoring the covariance matrix

A chart based on the sample variances S x2 = ∑ i =1 ( X i − µ x ) / n and S y2 = ∑ i =1 (Yi − µ y ) / n n

2

n

2

is considered to control Σ . As an alternative to the use of two S 2 charts we consider a single chart based on the statistic VMAX = max{S x2 , S y2 } , shortly VMAX chart. A signal is given if VMAX>CL, the control limit for the VMAX chart. According to Costa & Machado (2008), the power of the VMAX chart is given by: PS = 1 −



nCL a 2 0

⎡ Pr ⎢ χ n2, t ρ 2 ⎢⎣ (

1− ρ 2

)


n1 . The size of each sample depends on the sample point position of the preceding sample. If the sample point falls in the central region, that is, VMAX < WLi , then the next sample size should be n1 , and if the sample point falls in the warning region, that is, WLi < VMAX < CLi , then the next sample size should be n2 .

552

Pesquisa Operacional, v.29, n.3, p.547-562, Setembro a Dezembro de 2009

Machado, Costa & Claro – Monitoring bivariate processes

The properties of a control chart with VSS are also determined by the number of samples that are taken from the process until the control chart produces a signal. We suppose that the process remains in control with Σ = Σ 0 for a long time until the occurrence of the assignable cause. After the assignable cause occurrence the covariance matrix changes, that is, Σ = Σ 1 . During the in-control period, the rate of inspected items per sampling, n , is given by n = n1 p0 + n2 (1 − p0 )

(4)

being p0 the probability of a sample point falling in the central region, given that it did not fall in the action region. When the process is out-of-control, the ARL is obtained using Markov chains (see Costa & Machado (2008) for details). 4.2 The VMAX chart with double sampling

When the double sampling (DS) scheme is in use samples of size n = n1 + n2 are taken from the process at regular time intervals. The sampling is performed in two stages. At the first stage n1 items are inspected and their quality characteristics (X;Y) are measured. With the sample variances of (X;Y) the VMAX1 is computed, that is, VMAX1= max{S12x , S12y } . If VMAX1 < WL , the sampling is interrupted and the process is considered in control, where WL is the warning limit for the first stage. The control chart signals an out-of-control condition when VMAX1 > CL1 , where CL1 is the control limit for the first stage, CL1 > WL > 0 . If WL < VMAX1 ≤ CL1 , the sampling goes on to the second stage, where the remaining n2 items of the sample are inspected and the VMAX 2 = max{S 22x , S 22y } statistic is computed, being S 22x and S 22y the sample variances of (X;Y) considering the whole sample of size n . The control chart also signals an out-of-control condition when VMAX 2 > CL2 , where CL2 is the control limit for the second stage. If the chart produces a signal when the process is operating properly with Σ = Σ 0 , this signal is a false alarm. At the first stage the false alarm risk is given by α1 and, at the second stage, it is given by α 2 . If the chart signals after the occurrence of the assignable cause, that is, when Σ ≠ Σ 0 , this signal is a true alarm. At the first stage the power of the chart detecting an assignable cause is given by p1 and, at the second stage, it is given by p2 . This way, when the DS procedure is in use, the false alarm risk and the power of the VMAX chart are given by α = α1 + α 2 and p = p1 + p2 , respectively. During the in-control period, the rate of inspected items per sampling, n , is given by: n = n1 + n2 (1- p0 )

Pesquisa Operacional, v.29, n.3, p.547-562, Setembro a Dezembro de 2009

(5)

553

Machado, Costa & Claro – Monitoring bivariate processes

where p0 is the probability that the DS procedure ends at the first stage, that is p0 = 1 − Pr [WL < VMAX 1 < CL 1 ] .

(6)

When the DS procedure is in use, the average run length (ARL) also measures the efficiency of a control chart in detecting a process change. The expression to obtain the ARL for the VMAX chart with DS is shown in Machado & Costa (2008). From earlier studies (Daudin, 1992), we observed that the charts with DS perform better with CL1 = ∞ , that is, α1 = 0 . In this case, a VMAX1 point in the region above WL always triggers the inspection of the whole sample.

4.3 The synthetic VMAX chart

According to the synthetic procedure (see Wu & Spedding, 2000), the signal is based on the Conforming Run Length (CRL). The CRL is the number of samples taken from the process since the previous nonconforming sample until the occurrence of the next nonconforming sample. In the case of the absence of a previous nonconforming sample, the CRL is the number of samples taken from the beginning of the monitoring until the occurrence of a nonconforming sample. The signal is given when the CRL ≤ L, where L is a specified positive integer. Figure 1 shows the synthetic VMAX chart. The sample is classified as nonconforming when VMAXi falls beyond the control limit CL. Samples 9 and 13 are nonconforming (Figure 1). In this case, CRL = 4 (13th sample – 9th sample = 4). As CRL < L (=5), the synthetic VMAX chart signals an out-of-control condition. VMAX

CL

L=5

1 2 3 4 5 6 7 8 9 10 11 12 13 sample number (i) Figure 1 – The synthetic VMAX chart.

554

Pesquisa Operacional, v.29, n.3, p.547-562, Setembro a Dezembro de 2009

Machado, Costa & Claro – Monitoring bivariate processes

For many types of control charts (for instance, the Shewhart chart), the ARL value for a given shift in the mean and/or in the variance of the process does not depend on whether the assignable cause is affecting the process since the beginning of the monitoring (zero-state ARL) or whether it only occurs after the process has been running for some time (steady-state ARL) that is, the zero-state ARL and the steady-state ARL are equal. However, for charts based on conforming run lengths, such as synthetic charts, this does not happen. We assume that the assignable cause only occurs after the process has been running for some time, so the performance measure will be the steady-state ARL. Davis & Woodall (2002) obtained the steady-state ARL (SSARL) through a Markov chain model. We adopted this approach to obtain the steady-state ARL (SSARL) for the synthetic VMAX chart, see Machado et al. (2009) for details. 4.4 The EWMA scheme based on the VMAX statistic

The EWMA chart we propose to detect changes in the covariance matrix Σ is based on the statistic Z i = λ VMAX i + (1 − λ ) Z i −1 ,

i = 1,2,…,

(7)

where VMAX i = max{S x2i , S y2i } , being S x2i and S y2i the sample variances of X and Y, and λ

is a smoothing parameter, satisfying 0 < λ ≤ 1 . The value of λ determines the weight given to the current sample value. When λ = 1 the EWMA places all of it weights on the most recent observation. The starting value Z 0 is often taken to be the expected in-control value of Z (see Lucas & Saccucci, 1990). A signal is given if Z i > CL , the control limit for the EWMA chart. For the EWMA chart, the expected length of time from the process change to the signal will depend on the value of the control statistic Z at the time the change occurs. We assume that the control statistic has reached its steady-state or stationary distribution by the random point in time that the change occurs. When the steady-state distribution is adopted the chart’s performance is measured by the steady-state ARL, shortly SSARL. Machado & Costa (2008) used the Markov chain approach described in Saccucci & Lucas (1990) to obtain the SSARL values.

5. Comparing charts

Table 4 brings the ARLs of the S chart, the VMAX chart and the VMAX chart with VSS, DS, synthetic and EWMA procedures.

Pesquisa Operacional, v.29, n.3, p.547-562, Setembro a Dezembro de 2009

555

Machado, Costa & Claro – Monitoring bivariate processes

Table 4 – The ARL values for the S chart and for the VMAX charts ( ρ = 0.5, n = 4, n1 = 3 and n2 = 12). VMAX limits

S

Fixed

VSS

DS

Synthetic*

EWMA**

CL WL CL1

6.134 -

4.094 -

-

2.680

2.964 -

2.783 -

-

-

4.758



-

-

WL1

-

-

2.425

-

-

-

shifts

CL2

-

-

2.523

2.127

-

-

a =b

WL2

2

2

1.0 1.2 1.3 1.4 1.5 2 3 5

-

-

1.685

-

-

-

200.0 112.48 89.13 73.30 60.35 30.24 13.64 6.37

200.0 107.0 82.9 66.1 54.1 25.4 10.7 4.77

200.0 106.73 81.57 63.84 51.02 20.91 7.36 3.23

200.0 84.0 59.6 44.2 34.0 13.3 5.13 2.46

200.0 93.1 68.9 53.0 42.1 18.4 7.70 3.67

200.0 91.7 67.4 51.6 40.9 17.8 7.58 3.74

*L=5; ** λ =0.5

From Table 4 we conclude that the VMAX chart with DS is the one with the best performance (see the ARL in bold). 6. Control charts for monitoring the mean vector and the covariance matrix

(

variances S x2 = ∑in=1 ( X i − µ x )2 / n and S y2 = ∑in=1 Yi − µ y

)2 / n

μ

A chart based on the sample means X = ∑in=1 X i / n and Y = ∑in=1 Yi / n , and on the sample is considered to control

and Σ . As an alternative to the use of two X charts and two S 2 charts, we consider a single chart based on the statistic MCMAX = max{|Z x |,| Z y |,Wx , W y } , shortly MCMAX chart, where Z x = n ( X − µ x ) , Z y = n (Y − µ y ) , Wx = kS x2 and W y = kS y2 . The parameter k is required to attend the imposed condition that, during the in-control period, the four statistics (| Z x |, | Z y |, Wx , W y ) have the same probability to exceed CL, the control limit of the proposed chart. The power of the MCMAX chart is given by P = PM + PS ′ − PM PS ′ .

(8)

where PS ′ = 1 −

556



( n −1) CL ka 2 0

⎡ Pr ⎢ χ n2, t ρ 2 ⎣ (

(n − 1)CL ⎤

1− ρ 2

1

− t 2 [( n −1) / 2]−1 dt ) < kb 2 (1 − ρ 2 ) ⎥⎦ 2( n −1) 2 Γ[(n − 1) / 2] e t

(9)

Pesquisa Operacional, v.29, n.3, p.547-562, Setembro a Dezembro de 2009

Machado, Costa & Claro – Monitoring bivariate processes

and PM is given by expression (1). The control limit CL is obtained by expression (1), with

c = d = 0 and PM = 1 − 1 − α . The value of k is obtained by expression (9), with a = b = 1 and PS ′ = 1 − 1 − α . We used a grid search to obtain CL and k. The sample points plotted on the MCMAX chart correspond to the largest value among (| Z x |, | Z y |, Wx , W y ) . The MCMAX chart is not always faster in signaling than the joint T 2 and S charts. However, the synthetic procedure improves the performance of the MCMAX chart. According to this procedure, a second sample point beyond the upper control limit and not far than L sampling intervals from the first one triggers the alarm. We also adopted the approach proposed by Davis & Woodall (2002) to obtain the ARL for the synthetic MCMAX chart. Table 5 shows the influence of the design parameter L on the MCMAX performance. As L increases the speed with which the synthetic MCMAX chart signals also increases. The gain in speed is more significant when L increases from 1 to 6. For instance, when a = b = 1 and c = d = 0.5, the ARL reduces 18.0 % (from 21.4 to 17.5). By other hand, as L increases from 6 to 10, the ARL reduces 1.0 % (from 17.5 to 17.3). Table 5 – The influence of the L value on the synthetic MCMAX chart’s performance (n=5 and ρ = 0.5). L shifts

1

d

3

4

5

6

7

8

9

10

a

b

1.0

1.0

0 0 0 0.5 0.5 0.5

200.0 200.0 200.0 200.0 200.0 200.0 200.0 200.0 200.0 200.0 45.3 40.9 39.1 37.9 37.3 36.8 36.3 36.1 35.9 35.7 21.4 19.2 18.4 17.9 17.7 17.5 17.4 17.3 17.3 17.3

1.25 1.25

0 0 0 0.5 0.5 0.5

13.6 8.9 6.5

12.0 7.8 5.7

11.3 7.3 5.4

11.0 7.1 5.3

10.8 7.0 5.2

10.8 7.0 5.2

10.6 6.9 5.2

10.6 6.9 5.1

10.6 6.9 5.1

10.6 6.9 5.1

1.5

0 0 0 0.5 0.5 0 0.5 0.5

10.1 7.6 7.2 5.8

8.5 6.5 6.0 5.0

7.9 6.1 5.6 4.7

7.5 5.9 5.3 4.5

7.3 5.8 5.2 4.4

7.3 5.8 5.2 4.4

7.1 5.7 5.0 4.3

7.0 5.7 5.0 4.3

7.0 5.7 5.0 4.3

7.0 5.7 4.9 4.3

1.0

c

2

CL 2.350 2.469 2.536 2.582 2.617 2.644 2.667 2.687 2.703 2.719

Table 6 presents the ARLs of both the synthetic MCMAX chart and the joint T 2 and S charts. The synthetic MCMAX chart is faster in signaling than the joint charts, except when the correlation between X and Y is high. When the correlation is high, the joint charts are, in general, faster in signaling assignable causes that only affect the mean and/or the variability of one of the two quality characteristics. We chose the usual sample size n=5. The conclusions are the same considering other values of n.

Pesquisa Operacional, v.29, n.3, p.547-562, Setembro a Dezembro de 2009

557

Machado, Costa & Claro – Monitoring bivariate processes

Table 6 – ARL values of the joint T 2 and S charts and the synthetic MCMAX chart (n=5 and L = 7). Shifts (mean vector) Shifts (covariance matrix) a b ρ 1.0 1.0 0.0

1.25

1.5

1.25

1.5

1.0

1.0

1.25

1.5

c

0

d

0.5

0.5

0

0.75

0.75

0

1.0

1.0

0 200.0 200.0

0.5 0 48.9* 49.2 36.3** 36.3

0.5 19.9 17.4

0.75 16.7 10.8

0 16.8 10.8

0.75 5.4 5.1

1.0 6.8 4.3

0 6.6 4.3

1.0 2.3 2.3

0.5

200.0 200.0

38.4* 38.5 36.3** 36.3

34.8 17.4

10.8 10.8

10.8 10.8

10.6 5.1

4.1 4.3

4.2 4.3

4.1 2.3

0.7

200.0 200.0 42.6 28.2

20.6* 36.3** 22.2 15.0

20.2 36.3 17.5 12.1

39.7 17.4 11.0 8.6

5.7 10.8 11.2 7.3

5.6 10.8 8.9 6.2

12.9 5.1 4.4 3.8

2.3 4.3 5.2 3.7

2.3 4.3 4.7 3.5

5.0 2.3 2.1 2.1

0.5

41.6 26.5

17.7 14.4

14.1 11.7

15.8 8.4

7.6 7.2

6.6 6.1

7.1 3.8

3.5 3.6

3.5 3.5

3.5 2.1

0.7

39.1 25.1 15.3 7.5

11.9 14.0 10.7 5.9

9.8 11.4 9.0 5.2

17.1 8.2 6.7 4.5

4.6 7.1 6.9 4.2

4.2 6.0 5.8 3.8

8.2 3.8 3.5 3.8

2.2 3.6 4.0 2.8

2.2 3.6 3.7 2.7

4.1 2.1 2.0 1.9

0.5

14.6 7.1

8.9 5.7

7.7 5.0

8.4 4.3

5.2 4.1

4.6 3.7

5.0 2.8

3.0 2.7

2.9 2.7

3.0 1.9

0.7

13.1 6.7 15.6 10.2

6.9 5.5 9.9 6.7

5.8 4.9 9.8 6.7

8.3 4.2 6.8 5.0

3.6 4.0 6.2 4.4

3.4 3.6 6.3 4.4

5.3 2.7 3.5 2.9

2.0 2.7 3.9 2.9

2.1 2.6 3.9 2.9

3.3 1.9 2.0 1.9

0.5

15.9 10.6

8.7 6.9

8.6 6.9

8.6 5.1

4.9 4.5

4.8 4.5

4.9 3.0

2.9 2.9

2.9 2.9

2.9 1.9

0.7

15.7 11.2 4.4 3.0

6.9 7.2 3.6 2.6

6.8 7.2 3.6 2.6

9.2 5.3 3.1 2.4

3.5 4.6 3.0 2.3

3.5 4.6 3.0 2.3

5.4 3.0 2.3 1.9

2.0 3.0 2.4 1.9

2.0 3.0 2.4 1.9

3.3 1.9 1.7 1.6

0.5

4.3 3.2

3.5 2.8

3.5 2.8

3.4 2.5

2.7 2.4

2.7 2.4

2.7 2.0

2.1 2.0

2.1 2.0

2.1 1.6

0.7

4.3 3.4

3.1 2.9

3.1 2.9

3.6 2.6

2.3 2.5

2.3 2.5

2.9 2.0

1.7 2.0

1.7 2.0

2.2 1.6

0.0

0.0

0.0

0.0

0

* Joint T 2 and S charts **synthetic MCMAX chart

558

Pesquisa Operacional, v.29, n.3, p.547-562, Setembro a Dezembro de 2009

Machado, Costa & Claro – Monitoring bivariate processes

7. Illustrative example

The purpose of the following example is to show the ability of the synthetic MCMAX chart in detecting shifts in the mean vector and/or in the covariance matrix. To this end, we considered a bivariate process whose quality characteristics of interest, X and Y, are normally distributed. When the process is in-control, the mean vector and the covariance matrix are ⎛1 0⎞ given by µ '0 = (0, 0) and Σ0 = ⎜ ⎟. ⎝0 1⎠ We initially generate 5 samples of size 5 with the process in control. The last 5 samples were simulated considering that the assignable cause changes the mean and the variability of X, that is, µ x = 1.0 and σ x = 1.25 . Table 7 presents the data of X, Y, Z, W and MCMAX. The control limit of CL = 2.667 was determined by expression (1), with c=d=0, to assure a false alarm risk α of 0.005. Figure 2 shows the synthetic MCMAX control chart. As the number of sampling intervals l (=2) is smaller than L (=7), the synthetic MCMAX control chart signals an out-of-control condition. Table 7 – Values of X, Y, Z, W and MCMAX (L=7). Sample 1

-0.284 -0.457

-0.467 -0.758

-0.067 -0.709

0.875 -0.332

0.188 -0.436

Z 0.110 1.204

W* 0.049 0.006

MCMAX 1.204

2

0.213 -2.741

1.624 0.525

-0.781 0.078

-1.26 1.96

0.273 0.036

0.031 0.064

0.221 0.518

0.518

3

-1.51 0.034

0.247 2.198

0.355 -1.369

-1.187 -0.234

1.026 0.243

0.478 0.390

0.209 0.297

0.478

4

2.275 0.987

0.179 -1.154

0.125 -0.302

-1.685 0.437

-0.917 0.857

0.010 0.369

0.397 0.142

0.397

5

-0.086 0.29

0.999 2.512

-0.741 -0.143

-0.755 -0.693

0.799 -0.581

0.097 0.619

0.123 0.305

0.619

6

-0.706 0.18

1.381 -1.093

2.096 1.35

-0.021 1.527

1.891 -1.286

2.076 0.303

0.270 0.309

2.076

7

0.117 1.027

2.042 0.75

-0.725 -0.214

2.679 -0.186

0.45 2.047

2.041 1.531

0.353 0.158

2.041

8

1.489 0.565

1.652 -0.115

1.438 -0.868

1.193 0.899

1.691 0.19

3.338 0.300

0.007 0.082

3.338

9

0.725 0.866

0.066 -0.668

1.183 -0.242

0.632 2.329

2.135 0.029

2.120 1.035

0.107 0.250

2.120

10

2.468 -1.311

-2.009 0.785

3.845 0.86

3.552 -0.562

0.918 1.005

3.924 0.347

1.023 0.191

3.924

X Y

* k = 0.1782

Pesquisa Operacional, v.29, n.3, p.547-562, Setembro a Dezembro de 2009

559

Machado, Costa & Claro – Monitoring bivariate processes

MCMAX 4 3

CL =2.667 2 1 0 1

2

3

4 5 6 7 8 Sample number (i )

9

10

Figure 2 – The synthetic MCMAX chart – example.

8. Conclusions

In this article we proposed new charts to control the mean vector and/or the covariance matrix of bivariate processes. The monitoring statistics associated to these charts are based on the sample means and sample variances. As the control chart’s users are, in general, more familiar with means and variances, they will have less difficult to deal with the proposed charts if compared with the charts proposed by Hotelling (1947) and Alt (1985). In terms of efficiency, the VMAX chart always performs better than the S chart. With the rule of the two points beyond the control limit, the synthetic MCMAX chart almost ever performs better than the joint T 2 and S charts. Acknowledgements

This work was supported by CNPq – National Council for Scientific and Technological Development, Project 307744/2006-0 – and FAPESP – The State of São Paulo Research Foundation, Project 2006/00491-0. We are in debt with the three anonymous referees who carefully read an earlier draft and made many constructive suggestions. During the celebration of the 40th SBPO event, this paper was elected the best one in the field of Operations Management. We would like to thank the award committee. References

(1) Alt, F.B. (1985). Multivariate quality control. In: Encyclopedia of Statistical Science [edited by S. Kotz and N.L. Johnson], 6, 110-122. (2) Aparisi, F. (1996). Hotelling’s T 2 control chart with adaptive sample sizes. International Journal of Production Research, 34, 2853-2862. (3) Aparisi, F.; Jabaloyes, J. & Carrión, A. (2001). Generalized variance chart design with adaptive sample sizes. The bivariate case. Communication in Statistics – Simulation and Computation, 30, 931-948.

560

Pesquisa Operacional, v.29, n.3, p.547-562, Setembro a Dezembro de 2009

Machado, Costa & Claro – Monitoring bivariate processes

(4) Chen, G.; Cheng, S.W. & Xie, H. (2005). A new multivariate control chart for monitoring both location and dispersion. Communications in Statistics-Simulation and Computation, 34, 203-217. (5) Costa, A.F.B. (1994). X control charts with variable sample size. Journal of Quality Technology, 26, 155-163. (6) Costa, A.F.B. & Machado, M.A.G. (2007). Synthetic control charts with two-stage sampling for monitoring bivariate processes. Pesquisa Operacional, 27, 117-130. (7) Costa, A.F.B. & Machado, M.A.G. (2008). A new multivariate control chart for monitoring the covariance matrix of bivariate processes. Communications in StatisticsSimulation and Computation, 37, 1453-1465. (8) Costa, A.F.B. & Machado, M.A.G. (2008a). Bivariate control charts with double sampling. Journal of Applied Statistics, 35, 809-822. (9) Costa, A.F.B. & Machado, M.A.G. (2009). A new chart based on sample variances for monitoring the covariance matrix of multivariate processes. International Journal of Advanced Manufacturing Technology, 41, 770-779. (10) Daudin, J.J. (1992). Double sampling X charts. Journal of Quality Technology, 24, 78-87. (11) Davis, R.B. & Woodall, W.H. (2002). Evaluating and improving the synthetic control chart. Journal of Quality Technology, 34, 200-208. (12) Epprecht, E.K. & Costa, A.F.B. (2001). Adaptive sample size control charts for attributes. Quality Engineering, 13, 465-473. (13) Hotelling, H. (1947). Multivariate quality control, illustrated by the air testing of sample bombsights. Techniques of Statistical Analysis, 111-184. (14) Khoo. B.C. (2005). A new bivariate control chart to monitor the multivariate process mean and variance simultaneously. Quality Engineering, 17, 109-118. (15) Lucas, J.M. & Saccucci, M.S. (1990). Exponentially weighted moving average control schemes: properties and enhancements. Technometrics, 32, 1-16. (16) Machado, M.A.G. & Costa, A.F.B. (2008). The double sampling and the EWMA charts based on the sample variances. International Journal of Production Economics, 114, 134-148. (17) Machado, M.A.G. & Costa, A.F.B. (2008a). Some comments on the use of principal components and simultaneous univariate control charts for multivariate process control. Pesquisa Operacional, 28, 173-196. (18) Machado, M.A.G. & Costa, A.F.B. (2008b). Monitoring the mean vector and the covariance matrix of bivariate processes. Brazilian Journal of Operations and Production Management, 5, 47-62. (19) Machado, M.A.G.; De Magalhães, M.S. & Costa, A.F.B. (2008). Gráfico de controle de VMAX para o monitoramento da matriz de covariâncias. Produção, 18, 222-239. (20) Machado, M.A.G.; Costa, A.F.B. & Rahim, M.A. (2009). The synthetic control chart based on two sample variances for monitoring the covariance matrix. Quality and Reliability Engineering International, 25, 595-606.

Pesquisa Operacional, v.29, n.3, p.547-562, Setembro a Dezembro de 2009

561

Machado, Costa & Claro – Monitoring bivariate processes

(21) Microsoft Fortran Power Station 4.0 (1995). Professional Edition with Microsoft IMSL Mathematical and Statistical Libraries. Microsoft Corporation, Washington, USA. (22) Saccucci, M.S. & Lucas, J.M. (1990). Average run lengths for exponentially weighted moving average control schemes using the Markov chain approach. Journal of Quality Technology, 22, 154-162. (23) Zhang, G. & Chang, S.I. (2008). Multivariate EWMA control charts using individual observations for process mean and variance monitoring and diagnosis. International Journal of Production Research, 46, 6855-6881. (24) Wu, Z. & Spedding, T.A. (2000). A synthetic control chart for detecting small shifts in the process mean. Journal of Quality Technology, 32, 32-38.

562

Pesquisa Operacional, v.29, n.3, p.547-562, Setembro a Dezembro de 2009

Suggest Documents