Exemple pour la fonction Analyser un plan de surface de réponse

Un ingénieur en techniques d'emballage doit s'assurer que les fermetures des sachets en plastique qui contiennent un produit sont suffisamment solides pour empêcher les fuites, mais pas trop, afin que les clients puissent ouvrir les sachets. Ces sachets garantissent que des instruments chirurgicaux restent secs et stériles jusqu'à ce que quelqu'un les ouvre. L'ingénieur souhaite optimiser la résistance de la fermeture pour qu'elle soit comprise entre 20 et 32 livres (limites inférieure et supérieure), avec une valeur cible de 26 livres. L'ingénieur souhaite également réduire la variabilité de la résistance de la fermeture à 1 ou moins. L'ingénieur détermine que la température de la barre de la thermoscelleuse, le temps d'arrêt et la pression de la barre de la thermoscelleuse sont des facteurs qui contribuent à la résistance de la fermeture. L'ingénieur détermine également que la température de la barre de la thermoscelleuse, le temps d'arrêt et la température du matériau sont des facteurs importants pour déterminer la variation. L'ingénieur conçoit un plan de surface de réponse composite central pour examiner les facteurs qui influent sur la résistance et la variabilité de la fermeture. Il utilise la transformation par logarithme népérien pour analyser la variabilité de la résistance de fermeture.

L'ingénieur collecte des données et analyse le plan afin de déterminer quels facteurs ont un impact sur la résistance de la fermeture.

  1. Ouvrez le fichier de données échantillons, RésistFermeture.MTW.
  2. Sélectionnez Stat > DOE (plan d'expériences) > Surface de réponse > Analyser un plan de surface de réponse.
  3. Dans la zone Réponses, saisissez Force.
  4. Cliquez sur Graphiques.
  5. Sous Graphiques des valeurs résiduelles, sélectionnez Quatre en un.
  6. Cliquez sur OK dans chaque boîte de dialogue.

Interprétation des résultats

Dans le tableau d'analyse de la variance, les valeurs de p pour TempéBar*TempéBar, TpsArrêt*TpsArrêt et TempéBar*TpsArrêt sont significatives. L'ingénieur peut envisager de réduire le modèle pour retirer les termes non significatifs. Pour plus d'informations, reportez-vous à la rubrique Réduction du modèle.

La valeur R2 indique que le modèle explique 78,58 % de la variance de la résistance, ce qui indique qu'il est bien ajusté aux données. Le R2 prévu de 0 suggère que ce modèle est trop ajusté, ce qui justifie sa réduction.

Le diagramme de Pareto vous permet de repérer visuellement les effets importants et de comparer l'importance relative des divers effets. De plus, vous voyez que l'effet maximal est TempéBar*TpsArrêt (AB) car c'est celui qui a la barre la plus grande.

Dans le graphique des valeurs résiduelles, les points ne suivent pas véritablement la droite de Henry. La réduction du modèle peut résoudre ce problème.

Coefficients codés

TermeCoeffCoef ErTValeur de TValeur de pFIV
Constante28,441,6617,090,000 
TempéBar1,6850,8991,870,0791,00
TpsArrêt-1,7190,899-1,910,0741,00
PressBar1,4810,8991,650,1191,00
TempéMat1,0600,8991,180,2551,00
TempéBar*TempéBar-2,6620,823-3,230,0051,03
TpsArrêt*TpsArrêt-2,4760,823-3,010,0081,03
PressBar*PressBar-1,3290,823-1,610,1261,03
TempéMat*TempéMat-1,1510,823-1,400,1811,03
TempéBar*TpsArrêt-5,811,10-5,280,0001,00
TempéBar*PressBar-0,091,10-0,080,9381,00
TempéBar*TempéMat-0,141,10-0,130,9021,00
TpsArrêt*PressBar0,551,100,500,6241,00
TpsArrêt*TempéMat0,241,100,220,8321,00
PressBar*TempéMat-0,101,10-0,090,9291,00

Récapitulatif du modèle

SR carréR carré
(ajust)
R carré
(prév)
4,4022878,58%59,84%0,00%

Analyse de la variance

SourceDLSomCar ajustCM ajustValeur FValeur de p
Modèle141137,5181,2514,190,004
  Linéaires4218,6554,6622,820,060
    TempéBar168,1368,1293,520,079
    TpsArrêt170,9470,9393,660,074
    PressBar152,6252,6162,710,119
    TempéMat126,9626,9631,390,255
  Carré4372,0793,0184,800,010
    TempéBar*TempéBar1202,61202,61110,450,005
    TpsArrêt*TpsArrêt1175,32175,3189,050,008
    PressBar*PressBar150,5250,5222,610,126
    TempéMat*TempéMat137,8737,8661,950,181
  Interaction à 2 facteurs6546,7991,1324,700,006
    TempéBar*TpsArrêt1540,47540,47027,890,000
    TempéBar*PressBar10,120,1210,010,938
    TempéBar*TempéMat10,300,3050,020,902
    TpsArrêt*PressBar14,844,8400,250,624
    TpsArrêt*TempéMat10,900,8990,050,832
    PressBar*TempéMat10,160,1600,010,929
Erreur16310,0819,380   
  Inadéquation de l'ajustement10308,2030,82098,510,000
  Erreur pure61,880,313   
Total301447,60     

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

Force=-289,3 + 2,287 TempéBar + 206,6 TpsArrêt + 0,124 PressBar + 0,594 TempéMat
- 0,00426 TempéBar*TempéBar - 39,6 TpsArrêt*TpsArrêt - 0,000532 PressBar*PressBar
- 0,00288 TempéMat*TempéMat - 0,930 TempéBar*TpsArrêt - 0,000070 TempéBar*PressBar
- 0,00028 TempéBar*TempéMat + 0,0440 TpsArrêt*PressBar + 0,047 TpsArrêt*TempéMat
- 0,00010 PressBar*TempéMat

Ajustements et diagnostics pour les observations aberrantes

ObservationForceValeur
ajustée
RésiduelleVal.
résid.
norm.
320,6914,436,262,20R
527,4321,725,712,01R
925,9920,165,832,05R
1921,3815,106,282,21R
R : Valeur résiduelle élevée