http://dx.doi.org/10.1590/S0103-65132008000200003
The VMAX control chart for monitoring the covariance matrix
Marcela Aparecida Guerreiro MachadoI; Maysa Sacramento de MagalhãesII; Antônio Fernando Branco CostaI
IFEG/UNESP
IIENCE/IBGE
Neste artigo é proposto, para o monitoramento de processos normais bivariados, um gráfico de controle baseado nas variâncias amostrais de duas características de qualidade. Os pontos plotados no gráfico correspondem ao valor da maior variância amostral. O gráfico proposto, denominado gráfico de VMAX, tem um desempenho superior ao do gráfico da variância amostral generalizada |S| e, além disso, tem uma melhor capacidade de diagnóstico, ou seja, com ele é mais fácil identificar a variável que teve sua variabilidade alterada pela ocorrência da causa especial. Quando a amostragem dupla está em uso o gráfico proposto também tem um desempenho superior ao do gráfico de |S|, exceto em alguns casos em que o tamanho da segunda amostra é muito grande.
Palavras-chave: Gráficos de controle de VMAX, matriz de covariâncias, processos bivariados, amostragem dupla.
In this paper a control chart based on sample variances from two quality characteristics for monitoring bivariate normal processes is proposed. The points plotted on the chart are the maximum of the values of these two statistics. The proposed chart (VMAX chart) detects process disturbances faster than the generalized variance |S| chart and has a better diagnostic feature, that is, with the VMAX chart it is easier to relate an out-of-control signal to the variable whose variability has moved away from its in-control state. When the double sampling scheme is used the proposed chart has also a better performance, except in a few cases in which the size of the samples at the second stage is very large.
Key words: VMAX control chart, covariance matrix, bivariate processes, double sampling.
O gráfico de controle tem sido a principal ferramenta estatística utilizada no monitoramento de processos, graças a sua simplicidade operacional. De acordo com os fundamentos de Shewhart, sempre que um ponto é plotado na região de ação do gráfico, deve-se ir à procura de possíveis causas especiais que estejam afetando a qualidade do processo. Entretanto, a simplicidade operacional faz do gráfico de controle de Shewhart uma ferramenta lenta na detecção de pequenas a moderadas alterações no parâmetro do processo que está sendo monitorado. Diante disso, não têm sido poucas as inovações sugeridas ao trabalho original de Shewhart, como por exemplo, gráficos de somas acumuladas (PAGE, 1954; HAWKINS; OLWELL, 1998; REYNOLDS; STOUMBOS, 2005), gráficos de médias móveis ponderadas exponencialmente (LUCAS; SACCUCCI, 1990; DOMANGUE; PATCH, 1991; CHEN et al., 2001; COSTA; RAHIM, 2006), gráficos adaptativos (COSTA, 1994, 1999; De MAGALHÃES et al., 2002, 2006) e gráficos de controle com amostragem dupla (CROASDALE, 1974; DAUDIN, 1992; IRIANTO; SHINOZAKI, 1998; CAROT et al., 2002; HE et al., 2002; HE; GRIGORYAN, 2002, 2006; COSTA et al., 2005a).
A adoção de gráficos de controle com amostragem dupla para o monitoramento do vetor de médias de processos multivariados tem sido objeto de pesquisa recente. He e Grigoryan (2005) avaliaram o desempenho dos gráficos multivariados com amostragem dupla por meio de simulação. Alternativamente, Machado (2006) e Costa e Machado (2007) desenvolveram expressões para obtenção das propriedades dos gráficos bivariados com amostragem dupla.
Assim como é importante monitorar o vetor de médias de um processo, é também importante monitorar a sua matriz de covariâncias (YEH et al., 2003, 2004, 2005; SURTIHADI et al., 2004; GRIGORYAN; HE, 2005).
O primeiro gráfico de controle utilizado no monitoramento da matriz de covariâncias Σ se baseou na estatística obtida do teste da razão de máxima verossimilhança generalizada (ALT, 1985). Para o caso bivariado, Alt (1985) propôs o uso da variância amostral generalizada |S| para controlar a matriz de covariâncias Σ. É importante salientar que os testes estatísticos para igualdade de matrizes de covariâncias foram propostos bem antes de serem utilizados no monitoramento de processos. Em Anderson (2003), os leitores interessados encontrarão detalhes destes testes.
Aparisi et al. (1999) estenderam a aplicação da estatística |S| para o caso em que o número de variáveis do processo sob monitoramento é maior do que dois (p>2) e, considerando esta mesma estatística, propuseram um gráfico de controle bivariado com tamanho de amostra variável (APARISI et al., 2001). Grigoryan e He (2005) propuseram o gráfico de controle multivariado com amostragem dupla para o monitoramento da matriz de covariâncias. Uma revisão recente de gráficos de controle para o monitoramento da matriz de covariâncias foi elaborada por Yeh et al. (2006). Chang e Zhang (2007) apresentaram uma metodologia para o monitoramento da variabilidade de processos multivariados autocorrelacionados.
Neste artigo é proposto, para o monitoramento de processos bivariados, um gráfico de controle baseado nas variâncias amostrais de duas características de qualidade. Os pontos plotados no gráfico correspondem ao valor da maior variância. O gráfico proposto, denominado gráfico de VMAX, ou seja, da maior variância, tem um desempenho superior ao do gráfico da variância amostral generalizada |S|. Além do gráfico de VMAX, também é proposto o gráfico de VMAX com amostragem dupla. Este gráfico tem também um desempenho superior ao do gráfico da variância amostral generalizada |S|, exceto em alguns casos em que o tamanho da segunda amostra é muito grande. Nos Apêndices I e II são apresentadas as expressões para a obtenção das propriedades dos gráficos de VMAX e de VMAX com amostragem dupla, respectivamente. A principal vantagem do gráfico de VMAX em relação aos gráficos de controle multivariados é que ele possui uma melhor capacidade de diagnóstico, ou seja, com ele é mais fácil identificar a variável que teve sua variabilidade alterada pela ocorrência de uma causa especial.
2. GRÁFICOS UNIVARIADOS PARA O MONITORAMENTO DE PROCESSOS MULTIVARIADOS
Uma dificuldade encontrada ao se lidar com qualquer gráfico de controle multivariado é a interpretação prática de um sinal de fora de controle. Especificamente, não se sabe ao certo qual das variáveis (suponha p variáveis) ou qual subconjunto delas é responsável pelo sinal. Na literatura encontramos muitas maneiras para resolver este problema, desde métodos gráficos até técnicas mais analíticas. A técnica de decomposição é uma das de maior sucesso (MASON, YOUNG, 2002). Uma outra alternativa consiste em plotar gráficos univariados para controlar cada característica de qualidade de um processo normal multivariado (SEREL et al., 2000).
Como um exemplo de como controlar diferentes variáveis com gráficos de S2, considere duas características de qualidade, representadas pelas variáveis normalmente distribuídas X e Y. Primeiramente considere o caso em que as duas características são independentes e monitoradas separadamente por dois gráficos de S2 univariados com limites de controle de 4 desvios-padrão. Para amostras de tamanho cinco, cada ponto amostral tem uma probabilidade do erro do tipo I igual a 0,004 (ou seja, α =0,004) de exceder os limites de controle de 4 desvios-padrão. A probabilidade de que as variâncias amostrais estejam dentro dos limites quando o processo está em controle é igual a (1-0,004)(1-0,004). Então a probabilidade do erro total do tipo I (denotada por α'), para este caso, é dada por α' = 1-(1-0,004)(1-0,004)=1-0,992=0,008. No caso de p características de qualidade estatisticamente independentes, a probabilidade do erro total do tipo I será dada por,
Caso se deseje manter controlada a probabilidade do erro do tipo I para um processo p-variado com variáveis aleatórias independentes, a equação (1) pode ser utilizada para calcular a probabilidade do erro do tipo I de cada gráfico (MONTGOMERY, 2004). Conseqüentemente, obtém-se o coeficiente de abertura dos limites de controle de cada gráfico. É importante ressaltar que quando se trabalha com um número grande de variáveis, a probabilidade do erro do tipo I de cada gráfico separadamente é muito pequena e o fator de abertura dos limites de controle é muito grande, reduzindo, assim, a sensibilidade do gráfico de controle em detectar desajustes do processo.
Se as variáveis são dependentes, que é o caso mais comum, a obtenção dos limites de controle de cada gráfico deixa de ser trivial. Para simplificar, considere o caso bivariado em que as variáveis seguem uma distribuição normal e que a matriz de covariâncias com o processo em controle é dada
por , onde ρ é o coeficiente de correlação entre as variáveis X e Y. Se gráficos de S2 univariados são utilizados para controlar cada característica de qualidade, então a' será dado por:
Para o cálculo do fator de abertura do limite de controle, LC, foi necessário obter a expressão (A1) do Apêndice I.
3. GRÁFICO DE CONTROLE DA MAIOR VARIÂNCIA (GRÁFICO DE VMAX)
Nesta seção apresenta-se o gráfico de controle proposto para o monitoramento da matriz de covariâncias Σ de um processo bivariado. Supõe-se que o vetor de médias µ não sofre alterações.
A estatística de monitoramento corresponde ao maior valor das variâncias amostrais e , sendo X e Y duas características de qualidade que seguem uma distribuição normal bivariada. Sem perda de generalidade, supõe-se que com o processo em controle a matriz de covariâncias é dada por , sendo Σxy e Σyx as covariâncias entre X e Y. Na prática, o mais comum é encontrar variáveis com variâncias diferentes entre si e não necessariamente iguais a 1. Contudo, a análise de sensibilidade apresentada neste trabalho permanece válida, uma vez que o efeito da causa especial, que resulta em aumento da variabilidade, é sempre expresso em unidades das variâncias iniciais de X e Y, isto é, variâncias de X e Y estimadas com o processo em controle.
Outro ponto importante sobre a matriz de covariâncias Σ0 que merece atenção é a escolha do estimador a ser utilizado em sua estimação. De acordo com Djauhari (2005) e, segundo estudos que datam da década de 50 (ANDERSON, 2003), o estimador é viesado, sendo a média das matrizes de covariâncias amostrais Si, i = 1,2,...,m, das m amostras de tamanho n utilizadas na estimação de 0. Segue que ~ {|Σ0|/[m(n-1)]p} Zk, onde Zk são independentes e Zk ~ , k = 1,2,..., p. Para eliminar o viés, o estimador deve ser dividido por uma constante b3 que é função de m, n e p. É interessante mencionar que Djauhari (2005) usou a distribuição normal como aproximação da distribuição do determinante da matriz de covariâncias. Neste artigo, porém, o monitoramento é realizado com base nas variâncias amostrais, as quais têm distribuição de Qui-quadrado.
Analogamente a Surtihadi et al. (2004), assume-se duas possíveis maneiras de uma causa especial alterar a matriz de covariâncias, resultando na matriz .
Uma possibilidade (caso I) supõe que a causa especial afeta somente a variância da variável aleatória X, isto é, a = c e b = 1, ou somente a variância da variável aleatória Y, isto é, b = c e a = 1. A outra possibilidade (caso II) supõe que a causa especial altera tanto a variância de X como a de Y e, neste caso, a = b = √c, sendo c a magnitude da perturbação (c > 1). Em ambos os casos, a correlação entre X e Y, dada por , não é afetada pela causa especial. Da maneira como Σ0 foi definida, tem-se que Σ0x = Σ0y = 1 e, conseqüentemente, ρ = σxy = σyx.
Quando o gráfico de VMAX está em uso, amostras de tamanho n são retiradas do processo em intervalos de tempo regulares. Duas características de qualidade (Xi, Yi) das n unidades da amostra são medidas e a estatística é calculada, sendo as variâncias amostrais de (Xi, Yi).
Se a estatística VMAX for maior do que o limite de controle LC, o gráfico sinaliza um desajuste do processo. Após a ocorrência do sinal, o usuário pode imediatamente examinar as variâncias amostrais de X e Y para descobrir quais variáveis são responsáveis pelo desajuste do processo, ou seja, aquelas cujas variâncias amostrais são maiores que LC. Este diagnóstico deve ser feito com cautela, pois uma causa especial pode aumentar a variabilidade de X e Y e, para a amostra que sinaliza, não necessariamente, as variâncias amostrais de X e Y serão maiores que o limite de controle. O limite de controle LC do gráfico de VMAX é obtido utilizando a expressão (A1) do Apêndice I. A Figura 1 apresenta o gráfico proposto.
A eficiência de um gráfico de controle em detectar uma alteração no processo pode ser medida pelo número médio de amostras até o sinal (NMA), ver Costa et al. (2005). Durante o período em controle o NMA=1/α e é denominado NMA0. Quando um processo está sob controle é desejável que o número médio de amostras retiradas do processo desde o início do monitoramento até o sinal (NMA0) seja grande, de modo a garantir poucos alarmes falsos. Quando um processo está fora de controle é desejável que o número médio de amostras retiradas desde a ocorrência de uma causa especial até o sinal (NMA) seja pequeno, de modo a garantir uma rápida detecção de alterações no processo. No Apêndice I é desenvolvida a expressão para obtenção do NMA para o gráfico de VMAX.
O desempenho do gráfico proposto é comparado com o do gráfico da variância amostral generalizada |S|, onde tem distribuição de qui-quadrado com 2n - 4 graus de liberdade (ALT, 1985). A matriz de covariâncias S é dada por , sendo sxx e syy as variâncias amostrais de X e Y respectivamente e sxy e syx as covariâncias amostrais. O limite de controle para o gráfico de |S| é dado por:
A magnitude da perturbação é dada por c2, onde c2 é a razão entre |Σ1| e |Σ0|. Sendo assim, o gráfico de |S| terá o mesmo poder de detecção para qualquer perturbação na matriz de covariâncias que resulte no mesmo valor de c2, independente da causa especial afetar apenas a variável X ou apenas a variável Y (caso I), ou ambas (caso II). De acordo com a expressão (4) do NMA para o gráfico de |S|, o desempenho deste gráfico independe da correlação (APARISI et al., 2001):
O coeficiente de correlação ρ tem uma certa influência nas propriedades do gráfico de VMAX, ver Tabela 1. Desta forma, escolheu-se um valor de ρ intermediário para a obtenção das propriedades do gráfico de VMAX. De acordo com a Tabela 2, construída para ρ=0,5, o gráfico de VMAX é sempre mais eficiente que o gráfico da variância amostral generalizada |S|. Para uma mesma magnitude de perturbação (mesmo valor de c2), o gráfico de VMAX é mais sensível se a perturbação afetar apenas uma das variáveis (caso I). Isto se explica pelo fato da estatística de VMAX ser o valor da maior variância amostral. Para o caso I, o gráfico de VMAX tem um desempenho bem superior ao do gráfico de controle proposto por Alt (1985). Os gráficos de controle em todas as tabelas que seguem foram planejados com base em um erro do tipo I da ordem de 0,5%, ou seja, um NMA0 = 200,0.
4. DESCRIÇÃO DO GRÁFICO DE VMAX COM AMOSTRAGEM DUPLA
Suponha que o gráfico de controle de VMAX com amostragem dupla seja empregado no monitoramento da matriz de covariâncias Σ de duas características de qualidade (Xi , Yi), descritas por uma distribuição normal bivariada.
Similarmente ao gráfico de controle de Shewhart, amostras de tamanho n = n1 + n2 são retiradas do processo em intervalos de tempo regulares. A amostragem é realizada em dois estágios. No primeiro estágio n1 unidades são inspecionadas, as duas características de qualidade são medidas e a estatística é calculada, sendo e as variâncias amostrais de (Xi , Yi) da subamostra de tamanho n1.
Se VMAX1 for menor do que o limite de advertência, LA, a amostragem é interrompida. Se VMAX1 for maior do que o limite de controle do primeiro estágio, LC1, o gráfico de controle sinaliza um desajuste do processo. Se LA < VMAX1 < LC1, a amostragem vai para o segundo estágio, onde as n2 unidades restantes são inspecionadas e a estatística é calculada, sendo e as variâncias amostrais de (Xi , Yi) levando em consideração a amostra de tamanho n. O gráfico também sinaliza um desajuste do processo quando VMAX2 > LC2, sendo LC2 o limite de controle do segundo estágio.
A região dada por (LA, LC1) corresponde a região de advertência. A região acima de LC1 é a região de ação do primeiro estágio e a região acima de LC2 é a região de ação do segundo estágio.
Se um ponto cai na região de ação com o processo em controle, isto é, quando Σ = Σ0, tem-se um alarme falso. Em outras palavras, o gráfico de controle sinaliza erroneamente a existência de uma causa especial. No primeiro estágio a probabilidade de se ter um alarme falso é dada por α1 e, no segundo estágio, por α2.
Se um ponto cai na região de ação após a ocorrência da causa especial, isto é, quando Σ ≠ Σ0, tem-se um alarme verdadeiro. No primeiro estágio a probabilidade de se ter um alarme verdadeiro é dada por p1 e, no segundo estágio, por p2.
Deste modo, α = α1 + α2 é a probabilidade de se ter um alarme falso quando a amostragem dupla está em uso e p = p1 + p2 é o poder de detecção do gráfico de VMAX.
Durante o período em controle, o número médio de itens inspecionado por amostragem, , é dado por
onde p0 é a probabilidade de que a amostragem seja interrompida no primeiro estágio. Os valores de , n1 e n2 são especificados pelo usuário, onde n1 < < n2. Por exemplo, se = 5, n1 = 2 e n2 = 10, fazendo-se uso da expressão (5) tem-se que p0 = 0,7, ou seja, com o processo em controle, 70% das vezes a amostragem é interrompida no primeiro estágio.
A Figura 2 mostra um gráfico de VMAX com os valores de VMAX1 e VMAX2 plotados. O gráfico de controle dá um sinal quando VMAX1 > LC1 ou quando VMAX2 > LC2. Os pontos amostrais VMAX1, ao caírem na região de advertência, disparam a inspeção de toda a amostra, e, nestes casos, a estatística VMAX2 é calculada e comparada com LC2. Para evitar o uso de dois gráficos, um para o primeiro estágio e outro para o segundo estágio, o usuário pode construir o gráfico de VMAX com duas escalas, uma do lado esquerdo e uma do lado direito, como ilustrado na Figura 2. Os valores VMAX1 correspondem aos pontos pretos que são plotados considerando a escala da esquerda e os valores VMAX2 correspondem aos pontos brancos que são plotados considerando a escala da direita.
5. PROPRIEDADES DO GRÁFICO DE VMAX COM AMOSTRAGEM DUPLA
Com a amostragem dupla, a eficiência do gráfico de controle em detectar uma alteração no processo é medida também pelo NMA, número médio de amostras até o sinal. A expressão do NMA para o gráfico de VMAX com amostragem dupla encontra-se no Apêndice II.
As tabelas desta seção foram obtidas fixando LC1 = ∞, ou seja, α1 = 0. Neste caso, o gráfico de controle dá um sinal somente quando VMAX2 > LC2. Os pontos amostrais VMAX1, ao caírem na região de advertência, que para o caso em que LC1 = ∞ corresponde à região acima de LA, disparam a inspeção de toda a amostra, e, nestes casos, a estatística VMAX2 é calculada e comparada com LC2. Em geral, o desempenho do gráfico de VMAX tende a piorar sempre que se adotam valores de α1 > 0, ver Tabela 3.
A Tabela 4 compara os gráficos de VMAX e de |S| com amostragem dupla. Os valores da Tabela 4 para o gráfico de |S| com amostragem dupla foram obtidos por Grigoryan e He (2005) por meio de simulações, ver Apêndice III. De acordo com a Tabela 4, observa-se que quando a amostragem dupla está em uso o gráfico proposto tem, em geral, um desempenho superior ao do gráfico de |S|.
As Tabelas de 5 a 7 contêm os valores do NMA para o gráfico de VMAX com amostragem dupla. Nestas tabelas n1 é fixo enquanto n2 assume diferentes valores. Quando n2 aumenta, o gráfico de VMAX se torna mais ágil na detecção de pequenas alterações; porém, se torna lento na detecção de grandes alterações. Por exemplo, da Tabela 5 observa-se que para c2 = 1,3 o NMA diminui de 46,1 para 37,4 quando n2 aumenta de 8 para 16 (caso I). Por outro lado, para c2 = 5, o NMA aumenta de 1,53 para 1,67 quando n2 passa de 8 para 16 (caso I).
O efeito da escolha de n1 na velocidade com que o gráfico de VMAX sinaliza o desajuste pode ser observado comparando os NMAs das Tabelas 6 e 7. Quando n2 é pequeno e n1 aumenta de 3 para 4 o gráfico de VMAX se mostra mais ágil. Por exemplo, observa-se que para c2 = 1,1 e n2 =10, o NMA diminui de 65,6 para 65,0 quando n1 passa de 3 para 4 (caso I). Para a mesma perturbação, porém para n2 grande (n2 =20), o NMA aumenta de 73,8 para 75,3 quando n1 passa de 3 para 4. As Tabelas de 5 a 7 apresentam também os NMAs para o gráfico de |S| com amostragem dupla. Quando a amostragem dupla está em uso e apenas umas das variáveis é afetada pela causa especial (caso I), o gráfico de VMAX é sempre mais eficiente na detecção da perturbação do que o gráfico de |S|. Porém, se as duas variáveis são afetadas pela causa especial (caso II) e o valor de n2 é grande, o gráfico de VMAX tem desempenho inferior ao do gráfico de |S|, exceto para grandes perturbações. Vale salientar que nesses casos, os desempenhos dos dois gráficos não são substancialmente diferentes e que o gráfico de VMAX é bem mais simples de ser utilizado, pois o usuário está muito mais familiarizado com o cálculo das variâncias amostrais do que com o do determinante de uma matriz de covariâncias.
6. EXEMPLO ILUSTRATIVO
Nesta seção, considera-se o exemplo citado por Aparisi et al. (2001) para explicar o uso do gráfico de controle com amostragem dupla baseado na estatística VMAX ou na estatística |S|.
A peça da Figura 3 tem duas dimensões a serem controladas, X (distância entre centros de dois furos interiores) e Y (distância entre dois centros de furos laterais). As duas características de qualidade são correlacionadas, sendo ρ = 0,7
Quando o processo está em controle, o vetor de médias e a matriz de covariâncias do processo são conhecidos:
Inicialmente foram geradas 10 amostras com o processo em controle, onde = 4, n1 = 3 e n2 = 8. A última amostra foi gerada sob a condição do caso I ter ocorrido, em que a causa especial aumenta somente a variabilidade da primeira variável X, sendo c2 = 3,0. Portanto, . Na Tabela 8 estão os dados de (Xi , Yi) e na Tabela 9 estão as variâncias amostrais , , e e as estatísticas VMAX1, VMAX2, |S1| e |S2|. É importante salientar que |S1| é o determinante de S considerando uma amostra de tamanho n1 e |S2| de tamanho n, e que os valores plotados no gráfico de |S| são os das estatísticas |S1| e |S2| e não seus valores padronizados. Para a obtenção do limite de controle LC2 =1,103 e do limite de advertência LA =1,063 do gráfico de VMAX foi utilizada a expressão (5) e as expressões (A2) e (A5) do Apêndice II, para a = b = 1,0. O limite de controle LC2 =0,230 e o limite de advertência LA =0,099 do gráfico de |S| foram obtidos por simulação, de acordo com o esquema de amostragem dupla. Neste exemplo, adotou-se α =0,005 e LC1 = ∞.
As Figuras 4 e 5 apresentam os gráficos de VMAX e de |S| com amostragem dupla. O gráfico de VMAX sinaliza já na 1ª amostra após a ocorrência da causa especial. Entretanto, o mesmo não ocorre com o gráfico de |S|.
Neste artigo, considerou-se um gráfico de controle baseado nas variâncias amostrais de duas características de qualidade. Os pontos plotados no gráfico correspondem ao valor da maior variância. A principal razão em considerar o gráfico proposto é que ele possui uma melhor capacidade de diagnóstico, ou seja, com ele é mais fácil identificar a variável que teve sua variabilidade alterada pela ocorrência da causa especial.
O gráfico proposto é sempre mais eficiente que o gráfico da variância amostral generalizada |S|. Quando a amostragem dupla está em uso, observa-se que o gráfico proposto também é mais rápido do que o gráfico de |S|, exceto em alguns casos em que o tamanho da amostra no segundo estágio é muito grande.
Com a crescente informatização dos controles "on-line", tornou-se usual acompanhar as condições do processo por meio do monitoramento simultâneo de muitos parâmetros de qualidade. Sendo assim, a extensão natural do presente trabalho é a criação de um gráfico de VMAX para o monitoramento de processos com p variáveis de qualidade, sendo p>2. Os autores já estão trabalhando nesta direção. Uma alternativa para o controle de processos multivariados consiste em se reduzir o número de variáveis a serem monitoradas, através, por exemplo, da aplicação da técnica de análise de componentes principais (SOUZA; RIGÃO, 2005).
Os autores agradecem aos dois revisores anônimos pelos úteis comentários e sugestões. Esta pesquisa contou com o apoio da FAPESP e do CNPq.
ALT, F. B. Multivariate quality control. In: Kotz, S., Johnson, N. L., ed., Encyclopedia of Statistical Sciences. Wiley. 1985. [ Links ]
ANDERSON, T. W. An Introduction to Multivariate Statistical Analysis. Stanford: John Wiley. 2003. [ Links ]
APARISI, F.; JABALOYES, J.; CARRIÓN, A. Statistical properties of the |S| multivariate control chart. Communications in Statistics - Theory and Methods, v. 28, p. 2671-2686, 1999. [ Links ]
APARISI, F.; JABALOYES, J.; CARRIÓN, A. Generalized variance chart design with adaptive sample sizes. The bivariate case. Communications in Statistics - Simulation and Computation, v. 30, p. 931-948, 2001. [ Links ]
CAROT, V.; JABALOYES, J.M.; CAROT, T. Combined double sampling and variable sampling interval chart. International Journal of Production Research, v. 40, p. 2175-2186, 2002. [ Links ]
CHANG, S. I.; ZHANG, K. Statistical process control for variance shift detections of multivariate autocorrelated processes. Quality Technology of Quantitative Management, v. 4, p. 413-435, 2007. [ Links ]
CHEN, G.; CHENG, S. W.; XIE, H. Monitoring process mean and variability with one EWMA chart. Journal of Quality Technology, v. 33, p. 223-233, 2001. [ Links ]
COSTA, A. F. B. Joint and R charts with variable sample sizes and sampling intervals. Journal of Quality Technology, v. 31, p. 387-397, 1999. [ Links ]
COSTA, A. F. B.; MACHADO, M. A. G. Synthetic control chart with two-stage sampling for monitoring bivariate processes. Pesquisa Operacional, v. 27 , p. 117-130 , 2007. [ Links ]
COSTA, A. F. B.; RAHIM, M. A. A single EWMA chart for monitoring process mean and process variance. Quality Technology & Quantitative Management, v. 3, p. 295-305, 2006. [ Links ]
COSTA, A.F.B.; EPPRECHT E.K.; CARPINETTI, L.C.R. Controle Estatístico de Qualidade. 2. ed. São Paulo: Editora Atlas, 2005. 334 p. [ Links ]
COSTA, A.F.B.; De MAGALHÃES, M.S.; EPPRECHT E.K. The non-central chi-square chart with double sampling. Brazilian Journal of Operation and Production Management, v. 2, p. 57-80, 2005a. [ Links ]
CROASDALE, R. Control charts for a double sampling scheme based on average production run length. International Journal of Production Research, v. 12, p. 585-592, 1974. [ Links ]
DAUDIN, J. J. Double sampling charts. Journal of Quality Technology, v. 24, p. 78-87, 1992. [ Links ]
De MAGALHÃES, M. S.; COSTA, A. F. B.; EPPRECHT, E. K. Constrained optimization model for the design of an adaptive chart. International Journal of Production Research, v. 40, p. 3199-3218, 2002. [ Links ]
De MAGALHÃES, M. S.; COSTA, A. F. B.; MOURA NETO, F. D. Adaptive control charts: a Markovian approach for processes subject to independent out-of-control disturbances. International Journal of Production Economics, v. 99, p. 236-246, 2006. [ Links ]
DJAUHARI, M. A. Improved monitoring of multivariate process variability. Journal of Quality Technology, v. 37, p. 32-39, 2005. [ Links ]
DOMANGUE, R.; PATCH, SC. Some omnibus exponentially weighted moving average statistical process monitoring schemes. Technometrics, v. 33, p. 299-313, 1991. [ Links ]
GRIGORYAN, A.; HE, D. Multivariate double sampling |S| charts for controlling process variability. International Journal of Production Research, v. 43, p. 715-730, 2005. [ Links ]
HAWKINS, D. M.; OLWELL, D. H. Cumulative sum charts and charting for quality improvement. New York: Springer-Verlag. 1998. [ Links ]
HE, D; GRIGORYAN, A. Construction of double sampling S- control charts for agile manufacturing. Quality and Reliability Engineering International, v. 18, p. 343-355, 2002. [ Links ]
HE, D; GRIGORYAN, A. Multivariate multiple sampling charts. IIE Transactions, v. 37, p. 509-521, 2005. [ Links ]
HE, D; GRIGORYAN, A. Joint statistical design of double sampling and S charts. European Journal of Operational Research, v. 168, p. 122-142, 2006. [ Links ]
HE, D; GRIGORYAN, A; SIGH, M. Design of double and triple sampling control charts using genetic algorithms. International Journal of Production Research, v. 40, p. 1387-1404, 2002. [ Links ]
IRIANTO, D.; SHINOZAKI, N. An optimal double sampling control chart. International Journal of Industrial Engineering - Applications and Practice, v. 5, p. 226-234, 1998. [ Links ]
LUCAS, J. M.; SACCUCCI, M. S. Exponentially weighted moving average control schemes: properties and enhancements. Technometrics, v. 32, p. 1-16, 1990. [ Links ]
MACHADO, M. A. G. Estudo das propriedades dos gráficos de controle bivariados com amostragem dupla. Dissertação (Mestrado em Engenharia Mecânica) - Faculdade de Engenharia do Campus de Guaratinguetá, Universidade Estadual Paulista, Guaratinguetá, 129f., 2006. [ Links ]
MASON, R. L.; YOUNG, J. C. Multivariate Statistical Process Control with Industrial Applications. Philadelphia: ASA - SIAM. 2002. [ Links ]
MONTGOMERY, D. C. Introduction to Statistical Quality Control. New York: John Wiley & Sons, Inc. [ Links ] 2004.
PAGE, E. S. Continuous inspection schemes. Biometrika, v. 41, p. 100-115, 1954. [ Links ]
REYNOLDS, M. R., Jr.; STOUMBOS, Z. G. Should exponentially weighted moving average and cumulative sum charts be used with Shewhart limits? Technometrics, v. 47, p. 409-424, 2005. [ Links ]
SEREL, D. A.; MOSKOWITZ, H.; TANG, J. Univariate control charts for individual characteristics in a multinormal model. IIE Transactions, v. 32, p. 1115-1125, 2000. [ Links ]
SOUZA, A. M.; RIGÃO, M. H. Identificação de variáveis fora de controle em processos produtivos multivariados. Revista Produção, v. 15, p. 74-86, 2005. [ Links ]
SURTIHADI, J.; RAGHAVACHARI, M.; RUNGER, G. Multivariate control charts for process dispersion. International Journal of Production Research, v. 42, p. 2993-3009, 2004. [ Links ]
YEH, A. B.; LIN, D. K. J.; ZHOU, H.; VENKATARAMANI, C. A multivariate exponentially weighted moving average control chart for monitoring process variability. Journal of Applied Statistics, v. 30, p. 507-536, 2003. [ Links ]
YEH, A. B.; HUWANG, L.; WU, Y. F. A likelihood-ratio-based EWMA control chart for monitoring variability of multivariate normal processes. II Transactions, v. 36, p. 865-879, 2004. [ Links ]
YEH, A. B.; HUWANG, L.; WU, C. W. A multivariate EWMA control chart for monitoring process variability with individual observations. II Transactions, v. 37, p. 1023-1035, 2005. [ Links ]
YEH, A. B.; LIN, D. K. J.; MCGRATH, R. N. Multivariate control charts for monitoring covariance matrix: a review. Quality Technology and Quantitative Management, v. 3, p. 415-436, 2006. [ Links ]
Artigo recebido em 15/07/2007
Departamento de Produção/FEG/UNESP,
End.: Avenida Ariberto Pereira da Cunha, 333 - CEP 12516-410 - Bairro Pedregulho - Guaratinguetá - SP.