Exemple pour Analyser une réponse binaire pour un plan de criblage définitif

Des ingénieurs qualité veulent améliorer un procédé de confection de bretzels. La couleur est un critère de qualité essentiel. Ils utilisent un plan de criblage définitif pour déterminer les facteurs pouvant avoir un effet sur la couleur des bretzels. Pour l'expérience, les contrôleurs trient rapidement de petits lots de bretzels dans deux catégories : conforme et non conforme.

  1. Ouvrez le fichier de données échantillons, CouleurBretzels.MTW.
  2. Sélectionnez Stat > DOE (plan d'expériences) > Criblage > Analyser la réponse binaire
  3. Dans la zone Nom d'événement, saisissez Evénement.
  4. Dans la zone Nombre d'événements, saisissez Couleur correcte.
  5. Dans la zone Nombre d'essais, saisissez Essais.
  6. Cliquez sur Termes.
  7. Dans Inclure les termes suivants, sélectionnez Quadratique complet. Cliquez sur OK.
  8. Cliquez sur Pas à pas.
  9. Dans Méthode, sélectionnez Critères d'informations ascendantes.
  10. Cliquez sur OK dans chaque boîte de dialogue.

Interprétation des résultats

Le diagramme de Pareto montre les barres des termes issus du meilleur modèle selon le critère AICc. Deux effets principaux se trouvent dans le modèle : Temps de cuisson (E) et Température de cuisson 2 (H). Le modèle inclut également le terme au carré pour Temps de cuisson et l'effet de l'interaction entre les deux facteurs.

Les ingénieurs conviennent que ce modèle correspond à leurs connaissances du procédé. Ils décident de l'utiliser pour organiser une expérimentation approfondie.

Méthode

Fonction de liaisonLogit
Lignes utilisées50

Sélection ascendante des termes

AICc minimal atteint = 243,23

Informations de réponse

VariableValeurDénombrementNom
d'événement
Couleur correcteEvénement4235Evénement
  Non-événement765 
EssaisTotal5000 

Coefficients codés

TermeCoeffCoef ErTFIV
Constante2,3940,145 
Temps de cuisson0,73490,05381,11
Température de cuisson 20,54510,05411,20
Temps de cuisson*Temps de cuisson-0,3840,1531,04
Temps de cuisson*Température de cuisson 2-0,51060,05621,24

Rapports des probabilités de succès pour les prédicteurs continus

IncrémentRapport des
probabilités
de succès
IC à 95
%
Temps de cuisson2*(*; *)
Température de cuisson 215*(*; *)
Les rapports des probabilités de succès ne sont pas calculés pour les prédicteurs qui sont
     inclus dans les termes d'interaction, car ces rapports dépendent des valeurs des autres
     prédicteurs dans les termes d'interaction.

Récapitulatif du modèle

R carré de
la somme des
carrés des
écarts
R carré (ajust)
de la somme des
carrés des
écarts
AICAICcBIC
95,81%95,29%241,87243,23251,43

Tests d'adéquation de l'ajustement

TestDLKhi deuxValeur de p
Somme des carrés des écarts4532,280,922
Pearson4531,930,929
Hosmer-Lemeshow87,100,526

Analyse de la variance

SourceDLSomme des
carrés des
écarts
ajustée
Moyenne
ajustée
Khi deuxValeur de p
Modèle4737,452184,363737,450,000
  Temps de cuisson1203,236203,236203,240,000
  Température de cuisson 21100,432100,432100,430,000
  Temps de cuisson*Temps de cuisson16,7706,7706,770,009
  Temps de cuisson*Température de cuisson 2180,60580,60580,610,000
Erreur4532,2760,717   
Total49769,728     

Equation de régression en unités non codées

P(Evénement)=exp(Y')/(1 + exp(Y'))
Y'=-11,984 + 3,361 Temps de cuisson + 0,08740 Température de cuisson 2
- 0,0961 Temps de cuisson*Temps de cuisson
- 0,01702 Temps de cuisson*Température de cuisson 2

Ajustements et diagnostics pour les observations aberrantes

ObservationProbabilité
observée
Valeur
ajustée
RésiduelleVal.
résid.
norm.
10,98000,93762,02982,13R
70,98000,93961,95812,00R
240,90000,9497-2,0182-2,15R
R : Valeur résiduelle élevée