MSP
Statistiques industrielles
Management de la production et de la qualite
[email protected]
Universit
e de Caen Basse-Normandie
12 novembre 2015
[email protected] UCBN
MSP
12 novembre 2015
1 / 32
Premi`ere partie I
MSP
Capabilit
e
Capabilites des processus et des
syst`emes de mesures
Hypoth`
eses
Etude de capabilit
e
Capabilit
es
Int. Confiance
Index
Capabilite des processus
Index
[email protected] UCBN
MSP
12 novembre 2015
2 / 32
MSP
Chapitre
Capabilit
e
Hypoth`
eses
Etude de capabilit
e
Capabilit
es
Int. Confiance
Index
Capabilite des processus
Hypoth`eses
Etude de capabilite
Indices de capabilites
Intervalles de confiance
[email protected] UCBN
MSP
12 novembre 2015
3 / 32
MSP
Paragraphe
Capabilit
e
Hypoth`
eses
Etude de capabilit
e
Capabilit
es
Int. Confiance
Index
Capabilite des processus
Hypoth`eses
Etude de capabilite
Indices de capabilites
Intervalles de confiance
[email protected] UCBN
MSP
12 novembre 2015
4 / 32
MSP
Processus stable et sous-controle
Capabilit
e
Hypoth`
eses
Etude de capabilit
e
Capabilit
es
Int. Confiance
Index
On suppose que le processus est contr
ole par une statistique de
qualite quantitative Xt avec nt mesures pour lechantillon t. On
suppose que les mesures individuelles doivent etre toutes
comprises
I
entre une valeur maximale Ts (tolerance superieure) ou
USL (Upper Specification Limit) et
une valeur minimale Ti (tolerance inferieure) ou LSL
(Lower Specification Limit).
On construit une carte aux mesures avec un risque de premi`ere
esp`ece = 0.27%.
I
est stable : la loi de Xt est constante gaussienne de
moyenne 0 et decart type 0
est sous-contr
ole toutes les statistiques sont dans entre les
limites de contr
ole tYt [LCLt , UCLt ]
[email protected] UCBN
MSP
12 novembre 2015
5 / 32
MSP
Paragraphe
Capabilit
e
Hypoth`
eses
Etude de capabilit
e
Capabilit
es
Int. Confiance
Index
Capabilite des processus
Hypoth`eses
Etude de capabilite
Indices de capabilites
Intervalles de confiance
[email protected] UCBN
MSP
12 novembre 2015
6 / 32
MSP
Etude de capabilite
Capabilit
e
Hypoth`
eses
Etude de capabilit
e
Capabilit
es
Int. Confiance
Index
Etude de la capabilite dun procede
Sous les hypoth`eses precedentes, cest `a dire le procede nest
affecte que des causes communes, toute les causes speciales
pouvant affecter le processus ayant ete eliminees. Cette
situation correspond `a la meilleure performance de
fabrication possible.
I
Est ce que le processus est stable ?
Est ce que la distribution des mesures est une loi normale
constante ?
La proportion de non conformes, cest `a dire ne respectant
pas les tolerances.
Calcul dindices de capabilites mesurant la capacite du
procede `a respecter les tolerances.
[email protected] UCBN
MSP
12 novembre 2015
7 / 32
MSP
Exemple
Capabilit
e
Hypoth`
eses
Mesure dune longueur dune pi`ece
Etude de capabilit
e
Capabilit
es
Int. Confiance
Index
Une machine fabrique des pi`eces dont on mesure une
dimension. Cette axe doit respecter les specifications suivantes
78mm 12. 100 echantillons de 5 mesures ont ete effectuees.
Etudier la capabilite de cette machine.
Fichier :http://www.math.unicaen.fr/~kauffmann/data/msp-capabilite-etude.csv
1
2
3
4
5
6
X1
74.38
77.38
73.61
82.60
77.92
73.66
X2
74.40
76.86
73.33
77.28
74.28
83.24
X3
78.21
82.95
82.57
75.48
68.24
85.94
X4
80.01
72.82
83.99
75.28
82.82
82.30
X5
80.68
83.71
74.47
75.25
75.16
75.31
xbar
77.54
78.74
77.59
77.18
75.68
80.09
R
6.30
10.89
10.66
7.35
14.58
12.28
Table : Premiers echantillons
[email protected] UCBN
MSP
12 novembre 2015
8 / 32
MSP
Exercice
Capabilit
e
Hypoth`
eses
Etude de capabilit
e
Capabilit
es
Int. Confiance
Index
1. Etudier la stabilite, la maitrise statistique du processus ?
2. Calculer letendue de chaque echantillon
3. Calculer la moyenne x et letendue de chaque echantillon R.
4. Calculer la moyenne des moyennes x et letendue moyenne R.
5. Calculer une estimation de lecart type
b=s=
R
d2 (n=5)
6. Construire un histogramme de toutes les observations, ajouter la
densite dune loi normale de moyenne x et decart type s. Est ce
que la distribution des donnees est bien approximee par une loi
normale.
7. Construire un diagramme de type probability plot, doit-on
refuser lhypoth`ese que les observations suivent une loi normale ?
8. Calculer P1 = Pr ([X < Ti ]), P2 = Pr ([X > Ts ]) ou
X N (x, s). En deduire la proportion moyenne de pi`eces
defectueuses. Si on fabrique un million de pi`eces quel sera le
nombre moyen de pi`eces defectueuses.
[email protected] UCBN
MSP
12 novembre 2015
9 / 32
MSP
Exemple
Capabilit
e
Hypoth`
eses
Histogram of X.vecteur
Etude de capabilit
e
0.12
Capabilit
es
Int. Confiance
0.06
0.00
0.02
0.04
Density
0.08
0.10
Index
65
70
75
80
85
90
X.vecteur
Figure : Densite
[email protected] UCBN
MSP
12 novembre 2015
10 / 32
MSP
Exemple
Capabilit
e
Hypoth`
eses
Etude de capabilit
e
normal probability plot
Capabilit
es
Int. Confiance
0.99
probability
Index
0.95
0.90
0.80
0.70
0.60
0.50
0.40
0.30
0.20
0.10
0.05
0.01
normal
mean 76.78
sd 3.66
AD 0.33
pvalue 0.51
65 70 75 80 85 90
x
Figure : Probability plot
[email protected] UCBN
MSP
12 novembre 2015
11 / 32
MSP
Exemple
Capabilit
e
Hypoth`
eses
Etude de capabilit
e
Capabilit
es
Int. Confiance
Index
Ti
65
Ts
70
75
80
85
90
Figure : Densite
[email protected] UCBN
MSP
12 novembre 2015
12 / 32
MSP
Exemple
Capabilit
e
Hypoth`
eses
xbar Chart
for X[, variables]
Etude de capabilit
e
82
Capabilit
es
UCL
Int. Confiance
80
Index
78
CL
76
74
Group summary statistics
72
LCL
12
19
26
33
40
47
54
61
68
75
82
89
96
Group
Number of groups = 100
Center = 76.78358
StdDev = 3.65546
LCL = 71.87927
UCL = 81.68789
Number beyond limits = 1
Number violating runs = 3
Figure : Carte xbar
[email protected] UCBN
MSP
12 novembre 2015
13 / 32
MSP
Paragraphe
Capabilit
e
Hypoth`
eses
Etude de capabilit
e
Capabilit
es
Int. Confiance
Index
Capabilite des processus
Hypoth`eses
Etude de capabilite
Indices de capabilites
Intervalles de confiance
[email protected] UCBN
MSP
12 novembre 2015
14 / 32
MSP
Indice de capabilite Cp
Capabilit
e
Hypoth`
eses
Etude de capabilit
e
Capabilit
es
Int. Confiance
Index
On consid`ere un processus stable et sous-controle controle par
une mesure X avec un risque de premi`ere esp`ece de = 0.27%
On suppose que
I
la distribution de X peut etre modelisee par une normale
de moyenne et decart type .
Ces mesures doivent satisfaire les specifications suivantes
Ti X Ts .
On appelle indice de capabilite Cp du processus controle par X
Cp =
Ts Ti
6
Si lecart type est inconnu, on le remplace par une estimation
cp = Ts Ti
C
6b
[email protected] UCBN
MSP
12 novembre 2015
15 / 32
MSP
Indice Cp
Capabilit
e
Hypoth`
eses
Etude de capabilit
e
Capabilit
es
Int. Confiance
Index
Lindice de capabilite est le rapport entre
I la longueur Ts Ti de lintervalle de tol
erance [Ti , Ts ].
I la longueur UCL LCL = 6 de lintervalle de contr
ole
[LCL, UCL] dune carte aux mesures avec un risque de
premi`ere esp`ece = 0.27%.
Ts Ti
90 66
Cp =
=
= 1.28
6
87 69
LCL= 69
mu= 78
UCL= 87
TI= 66
65
[email protected] UCBN
TS= 90
70
MSP
75
80
85
90
12 novembre 2015
16 / 32
MSP
Processus Cp -capable
Capabilit
e
Hypoth`
eses
Etude de capabilit
e
Capabilit
es
Int. Confiance
Index
Definition
On dit que le processus stable et sous-controle ayant un indice
de capabilite Cp est Cp -capable si Cp 1, on peut le qualifier
de
2 Cp
4/3 < Cp 2
1.0 Cp < 4/3
Cp < 1
tr`es performant
bon
acceptable
inaceptable
Table : Qualification de la Cp -capabilite
Dans lexemple precedent la capabilite Cp est de 1.28, on peut
dire que le processus est acceptable.
[email protected] UCBN
MSP
12 novembre 2015
17 / 32
MSP
Indice de capabilite Cpk
Capabilit
e
Hypoth`
eses
Etude de capabilit
e
Avec les memes hypoth`eses que precedemment On appelle
indice de capabilite Cpk du processus contr
ole par X
Capabilit
es
Int. Confiance
Index
Cpl =
Ti
3
Cpk
Ts
3
= min(Cpl , Cpu )
Ti Ts
= min(
)
3
3
Cpu =
Si , sont inconnus on les remplace des estimations
b Ti Ts
b
bpk = min(
C
)
3b
3b
[email protected] UCBN
MSP
12 novembre 2015
18 / 32
MSP
Processus Cpk -capable
Capabilit
e
Hypoth`
eses
Etude de capabilit
e
Capabilit
es
Int. Confiance
Index
Definition
On dit que le processus stable et sous-controle ayant un indice
de capabilite Cpk est Cpk -capable si Cpk 1, on peut le
qualifier de
capabilite
2.5 < Cp
2 Cp < 2.5
4/3 Cp < 2
1.0 Cp < 4/3
Cp < 1
qualite
excellent
tr`es performant
bon
acceptable
inaceptable
intervalle
7.5
6
4
3
trop de pi`eces defecteuses
Table : Qualification de la Cpk -capabilite
[email protected] UCBN
MSP
12 novembre 2015
19 / 32
MSP
Exemple decentrage
Capabilit
e
Hypoth`
eses
Etude de capabilit
e
Capabilit
es
Int. Confiance
Index
Une specification pour une longueur dune monture est de
50mm 6mm La statistique de qualite longueur est supposee
statistiquement maitrisee et lestimation de lecart type est de
1.15mm. Par contre la moyenne va varier de 48.5 mm `a 53 mm,
le processus nest donc pas stable. On va voir que les indices de
capabilites peuvent nous aider `a detecter ce changement.
[email protected] UCBN
MSP
12 novembre 2015
20 / 32
MSP
Capabilit
e
Cp= 1.74
Cpk= 1.16
Hypoth`
eses
Etude de capabilit
e
Capabilit
es
Int. Confiance
Index
LCL= 44.6
mu= 48
TI= 44
[email protected] UCBN
UCL= 51.4
TS= 56
MSP
12 novembre 2015
21 / 32
MSP
Paragraphe
Capabilit
e
Hypoth`
eses
Etude de capabilit
e
Capabilit
es
Int. Confiance
Index
Capabilite des processus
Hypoth`eses
Etude de capabilite
Indices de capabilites
Intervalles de confiance
[email protected] UCBN
MSP
12 novembre 2015
22 / 32
MSP
Intervalle de confiance de Cp
Capabilit
e
Hypoth`
eses
Etude de capabilit
e
Capabilit
es
Int. Confiance
Index
On consid`ere que lon a processus stable controle par par la
statistique X Si la statistique de contr
ole X peut etre
modelisee par une loi normale de moyenne et decart type
alors un intervalle de confiance de niveau de Cp est :
s
s
2(1)/2,n1
21(1)/2,n1
cp
cp
Cp C
C
n1
n1
ou
I
2(1)/2,n1 est le quantile dordre (1 )/2 dune loi du
2 `a n 1 degres de libertes.
cp est estime `a partir de lecart type dechantillonage s
C
fonde sur n echantillons.
[email protected] UCBN
MSP
12 novembre 2015
23 / 32
MSP
Exemple
Capabilit
e
Hypoth`
eses
Etude de capabilit
e
Capabilit
es
Int. Confiance
Index
On dispose dune estimation du coeffient de capabilite de
cp = 1.58 calcule `a partir de n = 20 echantillons. Calculer un
C
intervalle de confiance de niveau 95% de Cp .
I
1 = 0.05. Il faut calculer les quantiles dordre
(1 )/2 = 2.5% et 1 (1 )/2 = 97.5% dune loi du
2 `a 19 degres de libertes.
20.025,n1 8.9, 20.975,n1 32.8
un
de confiance de niveau = 95% est de
q
q intervalle
8.9
32.8
[ 19 ,
19 ] 1.58 On a donc 95 chances sur cent que le
vrai coefficient de capabilite Cp se trouve dans lintervalle
[1.08, 2.076]
[email protected] UCBN
MSP
12 novembre 2015
24 / 32
MSP
cp
Valeurs minimales C
Capabilit
e
Hypoth`
eses
Etude de capabilit
e
Capabilit
es
Int. Confiance
Index
1
2
3
4
5
6
7
8
9
10
Cp (ref )
0.80
1.00
1.20
1.40
1.60
1.80
2.00
2.20
2.40
3.00
n=10
1.32
1.65
1.97
2.30
2.63
2.96
3.29
3.62
3.95
4.94
n=20
1.10
1.37
1.64
1.92
2.19
2.47
2.74
3.01
3.29
4.11
n=30
1.02
1.28
1.54
1.79
2.05
2.30
2.56
2.82
3.07
3.84
n=40
0.99
1.23
1.48
1.72
1.97
2.22
2.46
2.71
2.96
3.70
n=50
0.96
1.20
1.44
1.68
1.92
2.16
2.40
2.64
2.88
3.61
n=60
0.94
1.18
1.42
1.65
1.89
2.12
2.36
2.60
2.83
3.54
n=
0
1
1
1
1
2
2
2
2
3
cp requises pour avoir
Table : Valeurs minimales de C
H0 (Cp > Cp (ref )) avec un risque de se tromper de moins de = 5%
en fonction du nombre dechantillons.
[email protected] UCBN
MSP
12 novembre 2015
25 / 32
MSP
Intervalle de confiance de Cpk
Capabilit
e
Hypoth`
eses
Etude de capabilit
e
Capabilit
es
Int. Confiance
Index
Soit un processus stable contr
ole par par la statistique X qui
etre modelisee par r
une loi normale de moyenne et decart type
1
1
+ 2(n1)
alors une approximation
. En posant A =
d
2
9nCpk
qui de devient de plus en plus pr
ecise quand n devient
grand de l intervalle de confiance de niveau de Cpk est
d
d
C
pk (1 z1(1)/2 A) Cpk Cpk (1 + z1(1)/2 A)
z1(1)/2 est le quantile dordre 1 (1 )/2 dune loi
normale centree reduite
cp est estime `a partir de lecart type dechantillonage s
C
fonde sur n echantillons.
[email protected] UCBN
MSP
12 novembre 2015
26 / 32
MSP
Exemple
Capabilit
e
Hypoth`
eses
Etude de capabilit
e
Capabilit
es
Int. Confiance
Index
On dispose dune estimation du coeffient de capabilite de
d
C
e `a partir de n = 20 echantillons. Les
pk = 1.58 calcul
tolerances sont [740, 760], calculer un intervalle de confiance de
niveau 95% de Cpk .
I
1 = 0.05. Il faut calculer le quantile dordre
1 (1 )/2 = 97.5% dune loi du normale centree
reduite z0.975 1.98.
A 0.169
un intervalle de confiance de niveau = 95% est de
[1 1.98 0.169, 1 + 1.98 0.169] 1.58 = On a donc 95
chances sur cent que le vrai coefficient de capabilite Cpk se
trouve dans lintervalle
[1.03, 2.06]
[email protected] UCBN
MSP
12 novembre 2015
27 / 32
MSP
Capabilit
e
Hypoth`
eses
Etude de capabilit
e
Capabilit
es
Int. Confiance
Index
1
2
3
4
5
6
7
8
9
10
11
Cpk(ref)
0.80
1.00
1.21
1.40
1.60
1.80
2.00
1.34
2.50
3.00
4.00
n10
1.36
1.68
2.00
2.32
2.64
2.97
3.29
2.21
4.10
4.92
6.55
n20
1.13
1.40
1.66
1.93
2.20
2.48
2.75
1.84
3.43
4.11
5.47
n50
0.99
1.22
1.46
1.70
1.93
2.17
2.41
1.62
3.01
3.61
4.81
n100
0.93
1.15
1.37
1.60
1.82
2.05
2.28
1.52
2.84
3.41
4.54
n250
0.88
1.09
1.31
1.52
1.74
1.95
2.17
1.45
2.71
3.25
4.33
n500
0.86
1.07
1.28
1.49
1.70
1.91
2.12
1.42
2.64
3.17
4.23
n1000
0.84
1.05
1.25
1.46
1.67
1.88
2.08
1.39
2.60
3.12
4.16
d
Table : Valeurs minimales requises pour lindice C
pk permettant
daffirmer que le processus est Cpk -capable selon lindice souhaite
dans 95% des cas
[email protected] UCBN
MSP
12 novembre 2015
28 / 32
MSP
Capabilit
e
6
Hypoth`
eses
Etude de capabilit
e
Capabilit
es
Cpk.chap
Int. Confiance
Index
n10
n20
n50
n100
n250
n500
n1000
1.0
1.5
2.0
2.5
3.0
3.5
4.0
Cpk
d
Figure : C
pk minimale dans 95% des cas
[email protected] UCBN
MSP
12 novembre 2015
29 / 32
MSP
Chapitre
Capabilit
e
Hypoth`
eses
Etude de capabilit
e
Index
Capabilit
es
Int. Confiance
Index
[email protected] UCBN
MSP
12 novembre 2015
30 / 32
MSP
Index I
Capabilit
e
Hypoth`
eses
Etude de capabilit
e
Capabilit
es
Int. Confiance
Index
Capabilite
risque , 5
sous controle, 5
stabilite, 5
Tolerance, 5
inferieure, LSL , Ti, 5
superieure, USL,Ts, 5
Indice de capabilite
Cp
Cp capable, 17
definition, 16
exemple, 20
intervalle de confiance, 23
Intervalle de confiance Cp > Cp (ref ), 25
[email protected] UCBN
MSP
12 novembre 2015
31 / 32
MSP
Index II
Capabilit
e
Hypoth`
eses
Etude de capabilit
e
Capabilit
es
Int. Confiance
Index
Cpk
Cpk capable, 19
definition, 18
Intervalle de confiance, 26
Intervalle de confiance Cpk > Cpk (ref ), 28
[email protected] UCBN
MSP
12 novembre 2015
32 / 32