Fonctions > Plan d'expériences > Analyse de régression > Exemple : Qualité de l'ajustement
  
Exemple : Qualité de l'ajustement
Utilisez les fonctions polyfit et polyfitstat pour effectuer une régression linéaire et une analyse de la variance afin de tester la qualité de l'ajustement.
1. Définissez un tableau de données expérimentales relatif à un procédé de fabrication d'un polymère. La température de la réaction t et la vitesse d'alimentation en catalyseur fr affectent la viscosité vy du polymère.
Cliquez pour copier cette expression
2. Appelez la fonction polyfit afin de modéliser les données sous la forme d'une régression linéaire.
Cliquez pour copier cette expression
Cliquez pour copier cette expression
3. Calculez la viscosité prédite pour chaque valeur de la température et de la vitesse d'alimentation.
Cliquez pour copier cette expression
Cliquez pour copier cette expression
Cliquez pour copier cette expression
4. Calculez les résiduels (différence entre les valeurs calculées par le modèle et les valeurs mesurées).
Cliquez pour copier cette expression
5. Tracez les résiduels par rapport à la viscosité, la température et la vitesse d'alimentation observées.
Cliquez pour copier cette expression
Cliquez pour copier cette expression
Cliquez pour copier cette expression
Les tracés des résiduels indiquent que les variances de la viscosité et de la température observées augmentent selon l'accroissement, respectivement, de la viscosité et de la température.
6. Appelez polyfitstat pour calculer les différentes statistiques du modèle linéaire. Affichez la matrice ANOVA renvoyée par la fonction polyfitstat à la ligne 8.
Cliquez pour copier cette expression
Cliquez pour copier cette expression
Dans la matrice ANOVA, les sources de la variance sont divisées entre les composants de régression et de résiduels. Le composant de régression est par ailleurs divisé entre chaque coefficient de régression. Vous ne pouvez cependant pas faire la différence entre le manque d'ajustement et l'erreur pure pour le résiduel car les résultats de l'expérience vy ne comportent pas de répliques.
Calcul et utilisation de la table ANOVA pour la régression
1. Calculez la somme des carrés due à l'erreur (SSE).
Cliquez pour copier cette expression
La SSE est égale à χ2, qui est une mesure métrique générale d'adéquation de l'ajustement. Elle correspond à la quantité réduite lors du calcul d'une solution à l'aide de la méthode des moindres carrés. L'erreur correspond à une mesure qui indique la façon dont le modèle s'ajuste aux données. Elle indique dans quelle mesure la déviation n'est pas expliquée par la régression.
2. Définissez les degrés de liberté de l'erreur df_error par rapport aux degrés de liberté totaux df_total et aux degrés de liberté des paramètres df_param. Les degrés de liberté correspondent à la longueur des données, moins le nombre de paramètres d'ajustement.
Cliquez pour copier cette expression
Cliquez pour copier cette expression
Cliquez pour copier cette expression
3. Définissez la somme des carrés due à la régression (SSR) par rapport à la somme totale des carrés (SST).
Cliquez pour copier cette expression
Cliquez pour copier cette expression
4. Définissez l'erreur de carré moyenne (MSE) et l'erreur moyenne de régression (MSR). Divisez l'erreur par les degrés de liberté appropriés.
Cliquez pour copier cette expression
Cliquez pour copier cette expression
5. Réalisez une analyse du tableau de la variance afin de caractériser l'ajustement.
Somme des carrés
DF
Carré moyen
Facteur F
Régression
Cliquez pour copier cette expression
Cliquez pour copier cette expression
Cliquez pour copier cette expression
Cliquez pour copier cette expression
Erreur
Cliquez pour copier cette expression
Cliquez pour copier cette expression
Cliquez pour copier cette expression
Total
Cliquez pour copier cette expression
Cliquez pour copier cette expression
Vous pouvez comparer le tableau ci-dessus avec la matrice ANOVA polyfitstat.
Cliquez pour copier cette expression
6. Estimez dans quelle mesure le modèle s'ajuste aux données :
Cliquez pour copier cette expression
Cliquez pour copier cette expression
Ceci indique que 92.7 % de la variabilité de la viscosité s'explique par le modèle de régression linéaire appliqué.
7. Définissez le niveau de signification d'un test d'hypothèse afin de tester l'ajustement du modèle aux données.
Cliquez pour copier cette expression
8. Calculez la valeur critique F.
Cliquez pour copier cette expression
9. Testez l'hypothèse selon laquelle le modèle est ajusté aux données.
Cliquez pour copier cette expression
Acceptez l'hypothèse. Vous pouvez prévoir la viscosité du polymère grâce à ce modèle de régression linéaire.
Référence
Montgomery, D.C., Design and Analysis of Experiments, 5th ed., John Wiley & Sons, New York, 2001, pp. 398