|
Análise
Não-Paramétrica
A análise
não-paramétrica permite que a pessoa caracterize os seus dados de vida
sem supor uma distribuição subjacente. Isto pode ser vantajoso pois
evita alguns dos perigos inerentes em uma análise onde uma distribuição
imprópria é suposta. Entretanto, a análise é limitada as estimativas da
confiabilidade para a escala do conjunto de dados, fazendo com que se
torne impossível fazer predições da confiabilidade fora desta escala.
Enquanto isto limita um pouco a utilidade deste método, a análise
não-paramétrica pode ser usada para defrontar com a acurácia da análise
paramétrica. O software
Weibull++
da ReliaSoft contem a técnica de análise não-paramétrica.
Há diversos métodos
não-paramétricos para conduzir uma análise, como o método de Kaplan-Meier,
atuarial simples e atuarial padrão. Um método para calcular os intervalos
de confiança na análise dos resultados também pode ser desenvolvido.
A base da análise não-paramétrica para dados de vida é a função empírica
da f.d.c., que é dada por:

Estimador de Kaplan-Meier
O estimador de Kaplan-Meier, também é conhecido como estimador
produto limite, pode ser usado para calcular valores da confiabilidade não-paramétrica
para o conjunto de dados com múltiplas falhas e suspensões. A equação do
estimador é dada por:

onde m
é o total de números de pontos e n é o total de unidades
analisadas. A variável ni é definida por:

onde
rj é o número de falhas no j-ésimo grupo, e sj
é o número de suspensões do j-ésimo grupo. Note que, as
confiabilidades estimadas estão sendo calculadas somente para os tempo
que ocorreram uma ou mais falhas. Isto ocorre porque o cálculo do valor
de ni é realizado para os tempos onde ocorreram falhas
e suspensões, e é assumido que as suspensões ocorrem um pouco depois das
falhas, portanto os itens suspensos são considerados em operação e
inclusos na contagem dos ni.
Exemplo do Estimador
de Kaplan-Meier
Um grupo de 20 unidades
são postos em um teste de vida, obtendo os seguintes resultados:
| Número em Estado |
Estado (F ou S) |
Tempo do Estado Final |
| 3 |
F |
9 |
| 1 |
S |
9 |
| 1 |
F |
11 |
| 1 |
S |
12 |
| 1 |
F |
13 |
| 1 |
S |
13 |
| 1 |
S |
15 |
| 1 |
F |
17 |
| 1 |
F |
21 |
| 1 |
S |
22 |
| 1 |
S |
24 |
| 1 |
S |
26 |
| 1 |
F |
28 |
| 1 |
F |
30 |
| 1 |
S |
32 |
| 2 |
S |
35 |
| 1 |
S |
39 |
| 1 |
S |
41 |
Usando o estimador de
Kaplan-Meier para os dados, a tabela com as confiabilidades estimadas pode
ser construída.
| Tempo do Estado Final |
Número de Falhas,
ri |
Números de Suspensões, si |
Unidades Disponíveis,
ni |
 |
 |
| 9 |
3 |
1 |
20 |
0.850 |
0.850 |
| 11 |
1 |
0 |
16 |
0.938 |
0.797 |
| 12 |
0 |
1 |
15 |
1.000 |
0.797 |
| 13 |
1 |
1 |
14 |
0.929 |
0.740 |
| 15 |
0 |
1 |
12 |
1.000 |
0.740 |
| 17 |
1 |
0 |
11 |
0.909 |
0.673 |
| 21 |
1 |
0 |
10 |
0.900 |
0.605 |
| 22 |
0 |
1 |
9 |
1.000 |
0.605 |
| 24 |
0 |
1 |
8 |
1.000 |
0.605 |
| 26 |
0 |
1 |
7 |
1.000 |
0.605 |
| 28 |
1 |
0 |
6 |
0.833 |
0.505 |
| 30 |
1 |
0 |
5 |
0.800 |
0.404 |
| 32 |
0 |
1 |
4 |
1.000 |
0.404 |
| 35 |
0 |
1 |
3 |
1.000 |
0.404 |
| 39 |
0 |
1 |
2 |
1.000 |
0.404 |
| 41 |
0 |
1 |
1 |
1.000 |
0.404 |
Como pode ser
determinado da tabela precedente, as confiabilidades estimadas para
os tempos de falhas são:
| Failure Time |
Reliability Estimate |
| 9 |
85.0% |
| 11 |
79.7% |
| 13 |
74.0% |
| 17 |
67.3% |
| 21 |
60.5% |
| 28 |
50.5% |
| 30 |
40.4% |
Método Atuarial
Simples O
método atuarial simples é uma forma simples e fácil de realizar a
análise não-paramétrica de dados e pode ser utilizado para dados com
múltiplas censuras que estão arranjados em intervalos. Este método é
baseado no cálculo do número de falhas em um intervalo de tempo, rj,
versus o número de itens operando no período de tempo, ni.
A equação para estimar a confiabilidade pelo método atuarial é dado por:

onde m é o número
total de intervalo e n é o número total de unidades analisadas. A
variável ni é definida por:

onde rj
é o número de falhas no j-ésimo intervalo e sj
é o número de suspensões no j-ésimo intervalo.
Exemplo do Método
Atuarial Simples
Um grupo de 55 unidades é posto em um teste de vida onde a cada 50
horas as unidades são inspecionadas. Os resultados deste teste é
apresentado na tabela seguinte:
| Start Time |
End Time |
Number of Failures,
ri |
Number of Suspensions,
si |
| 0 |
50 |
2 |
4 |
| 50 |
100 |
0 |
5 |
| 100 |
150 |
2 |
2 |
| 150 |
200 |
3 |
5 |
| 200 |
250 |
2 |
1 |
| 250 |
300 |
1 |
2 |
| 300 |
350 |
2 |
1 |
| 350 |
400 |
3 |
3 |
| 400 |
450 |
3 |
4 |
| 450 |
500 |
1 |
2 |
| 500 |
550 |
2 |
1 |
| 550 |
600 |
1 |
0 |
| 600 |
650 |
2 |
1 |
As
confiabilidades estimadas pelo método atuarial simples pode ser obtida
pela expansão da tabela de dados para incluir os termos usados no
cálculo das estimativas da confiabilidade pela equação do atuarial
simples.
| Start Time |
End Time |
Number of Failures,
ri |
Number of Suspensions,
si |
Available Units, ni |
 |
 |
| 0 |
50 |
2 |
4 |
55 |
0.964 |
0.964 |
| 50 |
100 |
0 |
5 |
49 |
1.000 |
0.964 |
| 100 |
150 |
2 |
2 |
44 |
0.955 |
0.920 |
| 150 |
200 |
3 |
5 |
40 |
0.925 |
0.851 |
| 200 |
250 |
2 |
1 |
32 |
0.938 |
0.798 |
| 250 |
300 |
1 |
2 |
29 |
0.966 |
0.770 |
| 300 |
350 |
2 |
1 |
26 |
0.923 |
0.711 |
| 350 |
400 |
3 |
3 |
23 |
0.870 |
0.618 |
| 400 |
450 |
3 |
4 |
17 |
0.824 |
0.509 |
| 450 |
500 |
1 |
2 |
10 |
0.900 |
0.458 |
| 500 |
550 |
2 |
1 |
7 |
0.714 |
0.327 |
| 550 |
600 |
1 |
0 |
4 |
0.750 |
0.245 |
| 600 |
650 |
2 |
1 |
3 |
0.333 |
0.082 |
Como pode ser determinado
da tabela precedente, a confiabilidade estimada para os tempos de falhas
são:
Failure Period End Time |
Reliability Estimate |
| 50 |
96.4% |
| 150 |
92.0% |
| 200 |
85.1% |
| 250 |
79.8% |
| 300 |
77.0% |
| 350 |
71.1% |
| 400 |
61.8% |
| 450 |
50.9% |
| 500 |
45.8% |
| 550 |
32.7% |
| 600 |
24.5% |
| 650 |
8.2% |
Método
Atuarial Padrão O método atuarial padrão é uma variação do método
atuarial simples e envolve o ajuste do valor dos números de itens
operando em um intervalo. O Kaplan-Meier e o método atuarial simples
assumem que as suspensões em um período de tempo ou intervalo, ocorrem
no final do intervalo, após as falhas terem ocorrido. O método atuarial
padrão, assume que as suspensões ocorreram no meio do intervalo, o que
gera uma redução, pela metade, do número de itens disponíveis no
intervalo (suspensões), ou seja:

Com este ajuste, os
cálculos podem ser realizados como no método contábil simples,

Exemplo do Método
Atuarial Padrão Neste exemplo, nós aplicaremos o método
atuarial padrão à série de dados de usada no exemplo do método atuarial simples.
A solução deste exemplo é similar àquela do exemplo precedente,
com exceção da inclusão do termo ni', que é usado na equação
precedente. Aplicando a equação aos dados, nós podemos gerar a seguinte
tabela:
| Start Time |
End Time |
Number of Failures,
ri |
Number of Suspensions,
si |
Adjusted Units, ni' |
 |
 |
| 0 |
50 |
2 |
4 |
53 |
0.962 |
0.962 |
| 50 |
100 |
0 |
5 |
46.5 |
1.000 |
0.962 |
| 100 |
150 |
2 |
2 |
43 |
0.953 |
0.918 |
| 150 |
200 |
3 |
5 |
37.5 |
0.920 |
0.844 |
| 200 |
250 |
2 |
1 |
31.5 |
0.937 |
0.791 |
| 250 |
300 |
1 |
2 |
28 |
0.964 |
0.762 |
| 300 |
350 |
2 |
1 |
25.5 |
0.922 |
0.702 |
| 350 |
400 |
3 |
3 |
21.5 |
0.860 |
0.604 |
| 400 |
450 |
3 |
4 |
15 |
0.800 |
0.484 |
| 450 |
500 |
1 |
2 |
9 |
0.889 |
0.430 |
| 500 |
550 |
2 |
1 |
6.5 |
0.692 |
0.298 |
| 550 |
600 |
1 |
0 |
4 |
0.750 |
0.223 |
| 600 |
650 |
2 |
1 |
2.5 |
0.200 |
0.045 |
Como pode ser determinado
da tabela precedente, a confiabilidade estimada para os tempos de falhas
são:
Failure Period End Time |
Reliability Estimate |
| 50 |
96.2% |
| 150 |
91.8% |
| 200 |
84.4% |
| 250 |
79.1% |
| 300 |
76.2% |
| 350 |
70.2% |
| 400 |
60.4% |
| 450 |
48.4% |
| 500 |
43.0% |
| 550 |
29.8% |
| 600 |
22.3% |
| 650 |
4.5% |
Intervalo de Confiança
Não-Paramétrico Os limites de confiança para estimadores
de confiabilidade não-paramétricos podem ser calculados pelo método
similar utilizado nos limites de confiança paramétricos. A dificuldade,
no caso não-paramétrico, é a estimação da variância. Para estimar a
variância de dados não paramétricos, o Weibull++ utiliza a fórmula de Greenwood:

onde m é o número total de intervalos, e n é o número
total de itens analisados. A variável ni é definida
por:

onde
rj é quantidade de falhas no j-ésimo
intervalo e sj é o número de suspensões no
j-ésimo grupo. Desde que, a variância possa ser calculada, o desvio
padrão pode ser determinado pela raiz quadrada da variância,

Esta informação pode ser utilizada para determinar os limites de
confiança,,

onde:

e a é o nível de confiança desejado para limites
de confiança unilaterais.
Exemplo com Intervalo
de Confiança Não-Paramétrico Nesse exemplo, nós iremos determinar
o intervalo de confiança unilateral para confiabilidade estimada pelo
método atuarial padrão, com um intervalo de confiança de 97,5%.
Uma vez que, este tipo de problema é resolvido mais facilmente
construindo uma tabela similar a seguinte:
| Failure Time |
Reliability Estimate |
Number of Failures |
Adjusted Units |
ri/ni' |
Variance |
Error |
w |
Lower CL |
Upper CL |
| 50 |
0.962 |
2 |
53 |
0.0377 |
0.0007 |
0.0262 |
4.108 |
0.861 |
0.991 |
| 150 |
0.918 |
2 |
43 |
0.0465 |
0.0016 |
0.0397 |
2.797 |
0.799 |
0.969 |
| 200 |
0.844 |
3 |
37.5 |
0.0800 |
0.0030 |
0.0547 |
2.257 |
0.706 |
0.924 |
| 250 |
0.791 |
2 |
31.5 |
0.0635 |
0.0040 |
0.0630 |
2.107 |
0.642 |
0.888 |
| 300 |
0.762 |
1 |
28 |
0.0357 |
0.0045 |
0.0668 |
2.059 |
0.609 |
0.868 |
| 350 |
0.702 |
2 |
25.5 |
0.0784 |
0.0054 |
0.0737 |
1.996 |
0.542 |
0.825 |
| 400 |
0.604 |
3 |
21.5 |
0.1395 |
0.0068 |
0.0823 |
1.964 |
0.438 |
0.750 |
| 450 |
0.484 |
3 |
15 |
0.2000 |
0.0082 |
0.0907 |
2.038 |
0.315 |
0.656 |
| 500 |
0.430 |
1 |
9 |
0.1111 |
0.0091 |
0.0953 |
2.142 |
0.260 |
0.618 |
| 550 |
0.298 |
2 |
6.5 |
0.3077 |
0.0104 |
0.1020 |
2.602 |
0.140 |
0.524 |
| 600 |
0.223 |
1 |
4 |
0.2500 |
0.0100 |
0.1000 |
3.098 |
0.085 |
0.471 |
| 650 |
0.045 |
2 |
2.5 |
0.8000 |
0.0036 |
0.0599 |
15.69 |
0.003 |
0.423 |
Os resultados são ilustrados pelo gráfico seguinte:

|