Begin'R
Les statistiques avec R
Navigation
[Retour au sommaire]
# Régression linéaire multiple: Affiner le modèle :Objectifs * Choisir un bon modèle :Remarque : Importation et préparation du jeu de données ? {#rq_preparation_donnees_regression, toggle=collapse, title-display=show} Dans les parties sur la régression linéaire, deux jeux de données sont employés. Le jeu de données `ozone`pour les exemples et le jeu de données `her`pour les exercices. | Jeu de données | Excel | CSV | Txt | Présentation du jeu de données| |-----------------|--------|-------------|----------------|-------------| | Ozone | [](donnees/Ozone.xlsx) | [](donnees/ozone.csv) | [](donnees/ozone.txt) |[ozone](#ozoneexp)| |her| |[](donnees/her.csv) | |[her](#herexp)| [Rappel : comment importer un jeu de données](caps_2_3_importation_csv.html) :Aide{#ozoneexp, toggle=popup, title-display=hidden} Un jeu de données sur la pollution d'ozone, issu d'une étude de surveillance de la qualité de l'air à Rennes pendant l'été 2001, est considéré. Dans cette partie, l'objectif est d'expliquer la pollution d'ozone `maxO3` (concentration en $\mu g /m^3$) d'un jour donné à l'aide de variables explicatives : * des variables de température T9, T12, T15 (les suffixes 9, 12 et 15 désignant respectivement des mesures prises à 9h, 12h et 15h, pour toutes les variables) ; * des variables de nébulosité Ne9, Ne12, Ne15 ; * des variables de vent Vx9, Vx12, Vx15 ; * la mesure du maximum de la concentration en ozone de la veille maxO3v. :Exemple Télécharger le fichier `ozone` et importer le jeu de données dans RStudio {#eximport} :Corrigé{#import, toggle=collapse} Le jeu de données sur la pollution d'ozone peut être importé à partir du fichier `ozone.csv` comme suit : ```r ozone <- read.csv2("donnees/ozone.csv", header = TRUE, stringsAsFactors = TRUE) ozone ``` ``` ## X maxO3 T9 T12 T15 Ne9 Ne12 Ne15 Vx9 Vx12 Vx15 maxO3v ## 1 20010601 87 15.6 18.5 18.4 4 4 8 0.6946 -1.7101 -0.6946 84 ## 2 20010602 82 17.0 18.4 17.7 5 5 7 -4.3301 -4.0000 -3.0000 87 ## 3 20010603 92 15.3 17.6 19.5 2 5 4 2.9544 1.8794 0.5209 82 ## 4 20010604 114 16.2 19.7 22.5 1 1 0 0.9848 0.3473 -0.1736 92 ## 5 20010605 94 17.4 20.5 20.4 8 8 7 -0.5000 -2.9544 -4.3301 114 ## 6 20010606 80 17.7 19.8 18.3 6 6 7 -5.6382 -5.0000 -6.0000 94 ## 7 20010607 79 16.8 15.6 14.9 7 8 8 -4.3301 -1.8794 -3.7588 80 ## 8 20010610 79 14.9 17.5 18.9 5 5 4 0.0000 -1.0419 -1.3892 99 ## 9 20010611 101 16.1 19.6 21.4 2 4 4 -0.7660 -1.0261 -2.2981 79 ## 10 20010612 106 18.3 21.9 22.9 5 6 8 1.2856 -2.2981 -3.9392 101 ## 11 20010613 101 17.3 19.3 20.2 7 7 3 -1.5000 -1.5000 -0.8682 106 ## 12 20010614 90 17.6 20.3 17.4 7 6 8 0.6946 -1.0419 -0.6946 101 ## 13 20010615 72 18.3 19.6 19.4 7 5 6 -0.8682 -2.7362 -6.8944 90 ## 14 20010616 70 17.1 18.2 18.0 7 7 7 -4.3301 -7.8785 -5.1962 72 ## 15 20010617 83 15.4 17.4 16.6 8 7 7 -4.3301 -2.0521 -3.0000 70 ## 16 20010618 88 15.9 19.1 21.5 6 5 4 0.5209 -2.9544 -1.0261 83 ## 17 20010620 145 21.0 24.6 26.9 0 1 1 -0.3420 -1.5321 -0.6840 121 ## 18 20010621 81 16.2 22.4 23.4 8 3 1 0.0000 0.3473 -2.5712 145 ## 19 20010622 121 19.7 24.2 26.9 2 1 0 1.5321 1.7321 2.0000 81 ## 20 20010623 146 23.6 28.6 28.4 1 1 2 1.0000 -1.9284 -1.2155 121 ## 21 20010624 121 20.4 25.2 27.7 1 0 0 0.0000 -0.5209 1.0261 146 ## 22 20010625 146 27.0 32.7 33.7 0 0 0 2.9544 6.5778 4.3301 121 ## 23 20010626 108 24.0 23.5 25.1 4 4 0 -2.5712 -3.8567 -4.6985 146 ## 24 20010627 83 19.7 22.9 24.8 7 6 6 -2.5981 -3.9392 -4.9240 108 ## 25 20010628 57 20.1 22.4 22.8 7 6 7 -5.6382 -3.8302 -4.5963 83 ## 26 20010629 81 19.6 25.1 27.2 3 4 4 -1.9284 -2.5712 -4.3301 57 ## 27 20010630 67 19.5 23.4 23.7 5 5 4 -1.5321 -3.0642 -0.8682 81 ## 28 20010701 70 18.8 22.7 24.9 5 2 1 0.6840 0.0000 1.3681 67 ## 29 20010702 106 24.1 28.4 30.1 0 0 1 2.8191 3.9392 3.4641 70 ## 30 20010703 139 26.6 30.1 31.9 0 1 4 1.8794 2.0000 1.3681 106 ## 31 20010704 79 19.5 18.8 17.8 8 8 8 0.6946 -0.8660 -1.0261 139 ## 32 20010705 93 16.8 18.2 22.0 8 8 6 0.0000 0.0000 1.2856 79 ## 33 20010706 97 20.8 23.7 25.0 2 3 4 0.0000 1.7101 -2.7362 93 ## 34 20010707 113 17.5 18.2 22.7 8 8 5 -3.7588 -3.9392 -4.6985 97 ## 35 20010708 72 18.1 21.2 23.9 7 6 4 -2.5981 -3.9392 -3.7588 113 ## 36 20010709 88 19.2 22.0 25.2 4 7 4 -1.9696 -3.0642 -4.0000 72 ## 37 20010710 77 19.4 20.7 22.5 7 8 7 -6.5778 -5.6382 -9.0000 88 ## 38 20010711 71 19.2 21.0 22.4 6 4 6 -7.8785 -6.8937 -6.8937 77 ## 39 20010712 56 13.8 17.3 18.5 8 8 6 1.5000 -3.8302 -2.0521 71 ## 40 20010713 45 14.3 14.5 15.2 8 8 8 0.6840 4.0000 2.9544 56 ## 41 20010714 67 15.6 18.6 20.3 5 7 5 -3.2139 -3.7588 -4.0000 45 ## 42 20010715 67 16.9 19.1 19.5 5 5 6 -2.2981 -3.7588 0.0000 67 ## 43 20010716 84 17.4 20.4 21.4 3 4 6 0.0000 0.3473 -2.5981 67 ## 44 20010717 63 15.1 20.5 20.6 8 6 6 2.0000 -5.3623 -6.1284 84 ## 45 20010718 69 15.1 15.6 15.9 8 8 8 -4.5963 -3.8302 -4.3301 63 ## 46 20010719 92 16.7 19.1 19.3 7 6 4 -2.0521 -4.4995 -2.7362 69 ## 47 20010720 88 16.9 20.3 20.7 6 6 5 -2.8191 -3.4641 -3.0000 92 ## 48 20010721 66 18.0 21.6 23.3 8 6 5 -3.0000 -3.5000 -3.2139 88 ## 49 20010722 72 18.6 21.9 23.6 4 7 6 0.8660 -1.9696 -1.0261 66 ## 50 20010723 81 18.8 22.5 23.9 6 3 2 0.5209 -1.0000 -2.0000 72 ## 51 20010724 83 19.0 22.5 24.1 2 4 6 0.0000 -1.0261 0.5209 81 ## 52 20010725 149 19.9 26.9 29.0 3 4 3 1.0000 -0.9397 -0.6428 83 ## 53 20010726 153 23.8 27.7 29.4 1 1 4 0.9397 1.5000 0.0000 149 ## 54 20010727 159 24.0 28.3 26.5 2 2 7 -0.3420 1.2856 -2.0000 153 ## 55 20010728 149 23.3 27.6 28.8 4 6 3 0.8660 -1.5321 -0.1736 159 ## 56 20010729 160 25.0 29.6 31.1 0 3 5 1.5321 -0.6840 2.8191 149 ## 57 20010730 156 24.9 30.5 32.2 0 1 4 -0.5000 -1.8794 -1.2856 160 ## 58 20010731 84 20.5 26.3 27.8 1 0 2 -1.3681 -0.6946 0.0000 156 ## 59 20010801 126 25.3 29.5 31.2 1 4 4 3.0000 3.7588 5.0000 84 ## 60 20010802 116 21.3 23.8 22.1 7 7 8 0.0000 -2.3941 -1.3892 126 ## 61 20010803 77 20.0 18.2 23.6 5 7 6 -3.4641 -2.5981 -3.7588 116 ## 62 20010804 63 18.7 20.6 20.3 6 7 7 -5.0000 -4.9240 -5.6382 77 ## 63 20010805 54 18.6 18.7 17.8 8 8 8 -4.6985 -2.5000 -0.8682 63 ## 64 20010806 65 19.2 23.0 22.7 8 7 7 -3.8302 -4.9240 -5.6382 54 ## 65 20010807 72 19.9 21.6 20.4 7 7 8 -3.0000 -4.5963 -5.1962 65 ## 66 20010808 60 18.7 21.4 21.7 7 7 7 -5.6382 -6.0622 -6.8937 72 ## 67 20010809 70 18.4 17.1 20.5 3 6 3 -5.9088 -3.2139 -4.4995 60 ## 68 20010810 77 17.1 20.0 20.8 4 5 4 -1.9284 -1.0261 0.5209 70 ## 69 20010811 98 17.8 22.8 24.3 1 1 0 0.0000 -1.5321 -1.0000 77 ## 70 20010812 111 20.9 25.2 26.7 1 5 2 -1.0261 -3.0000 -2.2981 98 ## 71 20010813 75 18.8 20.5 26.0 8 7 1 -0.8660 0.0000 0.0000 111 ## 72 20010814 116 23.5 29.8 31.7 1 3 5 1.8794 1.3681 0.6946 75 ## 73 20010815 109 20.8 23.7 26.6 8 5 4 -1.0261 -1.7101 -3.2139 116 ## 74 20010819 67 18.8 21.1 18.9 7 7 8 -5.3623 -5.3623 -2.5000 86 ## 75 20010820 76 17.8 21.3 24.0 7 5 5 -3.0642 -2.2981 -3.9392 67 ## 76 20010821 113 20.6 24.8 27.0 1 1 2 1.3681 0.8682 -2.2981 76 ## 77 20010822 117 21.6 26.9 28.6 6 6 4 1.5321 1.9284 1.9284 113 ## 78 20010823 131 22.7 28.4 30.1 5 3 3 0.1736 -1.9696 -1.9284 117 ## 79 20010824 166 19.8 27.2 30.8 4 0 1 0.6428 -0.8660 0.6840 131 ## 80 20010825 159 25.0 33.5 35.5 1 1 1 1.0000 0.6946 -1.7101 166 ## 81 20010826 100 20.1 22.9 27.6 8 8 6 1.2856 -1.7321 -0.6840 159 ## 82 20010827 114 21.0 26.3 26.4 7 4 5 3.0642 2.8191 1.3681 100 ## 83 20010828 112 21.0 24.4 26.8 1 6 3 4.0000 4.0000 3.7588 114 ## 84 20010829 101 16.9 17.8 20.6 7 7 7 -2.0000 -0.5209 1.8794 112 ## 85 20010830 76 17.5 18.6 18.7 7 7 7 -3.4641 -4.0000 -1.7321 101 ## 86 20010831 59 16.5 20.3 20.3 5 7 6 -4.3301 -5.3623 -4.5000 76 ## 87 20010901 78 17.7 20.2 21.5 5 5 3 0.0000 0.5209 0.0000 59 ## 88 20010902 76 17.3 22.7 24.6 4 5 6 -2.9544 -2.9544 -2.0000 78 ## 89 20010903 55 15.3 16.8 19.2 8 7 5 -1.8794 -1.8794 -2.3941 76 ## 90 20010904 71 15.9 19.2 19.5 7 5 3 -6.1284 0.0000 -1.3892 55 ## 91 20010905 66 16.2 18.9 19.3 2 5 6 -1.3681 -0.8682 1.7101 71 ## 92 20010906 59 18.3 18.3 19.0 7 7 7 -3.9392 -1.9284 -1.7101 66 ## 93 20010907 68 16.9 20.8 22.5 6 5 7 -1.5000 -3.4641 -3.0642 59 ## 94 20010908 63 17.3 19.8 19.4 7 8 8 -4.5963 -6.0622 -4.3301 68 ## 95 20010912 78 14.2 22.2 22.0 5 5 6 -0.8660 -5.0000 -5.0000 62 ## 96 20010913 74 15.8 18.7 19.1 8 7 7 -4.5963 -6.8937 -7.5175 78 ## 97 20010914 71 15.2 17.9 18.6 6 5 1 -1.0419 -1.3681 -1.0419 74 ## 98 20010915 69 17.1 17.7 17.5 6 7 8 -5.1962 -2.7362 -1.0419 71 ## 99 20010916 71 15.4 17.7 16.6 4 5 5 -3.8302 0.0000 1.3892 69 ## 100 20010917 60 13.7 14.0 15.8 4 5 4 0.0000 3.2139 0.0000 71 ## 101 20010918 42 12.7 14.3 14.9 8 7 7 -2.5000 -3.2139 -2.5000 60 ## 102 20010919 65 14.8 16.3 15.9 7 7 7 -4.3301 -6.0622 -5.1962 42 ## 103 20010920 71 15.5 18.0 17.4 7 7 6 -3.9392 -3.0642 0.0000 65 ## 104 20010921 96 11.3 19.4 20.2 3 3 3 -0.1736 3.7588 3.8302 71 ## 105 20010922 98 15.2 19.7 20.3 2 2 2 4.0000 5.0000 4.3301 96 ## 106 20010923 92 14.7 17.6 18.2 1 4 6 5.1962 5.1423 3.5000 98 ## 107 20010924 76 13.3 17.7 17.7 7 7 6 -0.9397 -0.7660 -0.5000 92 ## 108 20010925 84 13.3 17.7 17.8 3 5 6 0.0000 -1.0000 -1.2856 76 ## 109 20010927 77 16.2 20.8 22.1 6 5 5 -0.6946 -2.0000 -1.3681 71 ## 110 20010928 99 16.9 23.0 22.6 6 4 7 1.5000 0.8682 0.8682 77 ## 111 20010929 83 16.9 19.8 22.1 6 5 3 -4.0000 -3.7588 -4.0000 99 ## 112 20010930 70 15.7 18.6 20.7 7 7 7 0.0000 -1.0419 -4.0000 83 ## vent pluie ## 1 Nord Sec ## 2 Nord Sec ## 3 Est Sec ## 4 Nord Sec ## 5 Ouest Sec ## 6 Ouest Pluie ## 7 Ouest Sec ## 8 Nord Sec ## 9 Nord Sec ## 10 Ouest Sec ## 11 Nord Sec ## 12 Sud Sec ## 13 Sud Sec ## 14 Ouest Pluie ## 15 Nord Sec ## 16 Ouest Sec ## 17 Ouest Sec ## 18 Nord Sec ## 19 Est Sec ## 20 Sud Sec ## 21 Nord Sec ## 22 Est Sec ## 23 Sud Sec ## 24 Ouest Sec ## 25 Ouest Pluie ## 26 Sud Sec ## 27 Ouest Sec ## 28 Nord Sec ## 29 Est Sec ## 30 Sud Sec ## 31 Ouest Sec ## 32 Sud Pluie ## 33 Nord Sec ## 34 Ouest Pluie ## 35 Ouest Pluie ## 36 Ouest Sec ## 37 Ouest Sec ## 38 Ouest Sec ## 39 Ouest Pluie ## 40 Est Pluie ## 41 Ouest Pluie ## 42 Ouest Pluie ## 43 Sud Sec ## 44 Ouest Pluie ## 45 Ouest Pluie ## 46 Nord Sec ## 47 Ouest Pluie ## 48 Sud Sec ## 49 Ouest Sec ## 50 Nord Sec ## 51 Nord Sec ## 52 Ouest Sec ## 53 Nord Sec ## 54 Nord Sec ## 55 Ouest Sec ## 56 Sud Sec ## 57 Ouest Sec ## 58 Nord Sec ## 59 Est Sec ## 60 Sud Pluie ## 61 Ouest Pluie ## 62 Ouest Pluie ## 63 Sud Pluie ## 64 Ouest Sec ## 65 Ouest Pluie ## 66 Ouest Pluie ## 67 Nord Pluie ## 68 Nord Sec ## 69 Ouest Pluie ## 70 Ouest Sec ## 71 Nord Sec ## 72 Sud Sec ## 73 Sud Sec ## 74 Ouest Pluie ## 75 Ouest Pluie ## 76 Sud Sec ## 77 Sud Pluie ## 78 Ouest Sec ## 79 Ouest Sec ## 80 Sud Sec ## 81 Ouest Sec ## 82 Est Sec ## 83 Est Sec ## 84 Nord Pluie ## 85 Ouest Sec ## 86 Ouest Pluie ## 87 Nord Pluie ## 88 Ouest Pluie ## 89 Ouest Pluie ## 90 Nord Pluie ## 91 Nord Pluie ## 92 Nord Pluie ## 93 Ouest Pluie ## 94 Ouest Sec ## 95 Ouest Sec ## 96 Ouest Pluie ## 97 Nord Pluie ## 98 Nord Pluie ## 99 Nord Sec ## 100 Nord Pluie ## 101 Nord Pluie ## 102 Ouest Pluie ## 103 Ouest Sec ## 104 Est Pluie ## 105 Est Sec ## 106 Nord Sec ## 107 Ouest Pluie ## 108 Sud Sec ## 109 Sud Pluie ## 110 Sud Sec ## 111 Ouest Pluie ## 112 Sud Sec ``` :Aide{#herexp, toggle=popup, title-display=hidden} Le jeu de données `her` rassemble * des données d'identité (numéro d'indentification, sexe, age) * des caractéristiques du corps humain (poids, taille, pouls, pression systolique, pression diastolique, taux de cholestérol, indice de masse corporel, circonférence de jambe, coude, poignet et bras). Le critère utilisé par défaut dans R pour choisir un "bon modèle" est le critère AIC (pour “An Information Criterion”, proposé par Akaike). La “philosophie”" de ce critère est de trouver un compromis entre le modèle le plus simple possible (i.e. contenant le moins de variables explicatives) et qui ajuste le mieux les données. Ce critère étant une fonction croissante des deux grandeurs $p^*$ et $\mathrm{RSS}^*$, le “meilleur” modèle (au sens du critère AIC) sera alors le modèle ayant la valeur du AIC la plus faible. :Compléments Pour aller plus loin sur le critère AIC {#AIC, toggle=collapse} La formule du critère AIC (sur lequel la méthodologie de sélection de variables est fondée) est la suivante : $$\mathrm{AIC} = n \log \left( \frac{\mathrm{RSS}^*}{n} \right) + 2(p^*+1)$$ où $p^*$ correspond au nombre de variables explicatives considérées dans le modèle courant (i.e. celui pour lequel l'AIC est calculé), et $RSS^* = \sum_{i=1}^n (y_i − \hat{y}_i^*)^2$ est la somme des carrés des résidus du modèle courant ($\mathrm{RSS}$ pour **Residual Sum of Squares** en anglais). :Remarque Notons que la quantité $\frac{\mathrm{RSS}}{n}$ correspond à l'estimation du paramètre $\sigma^2$. L’exemple de la sélection pas à pas (stepwise selection) descendante est considéré ici. A chaque étape, une étape de sélection descendante (où on cherche à enlever une variable du modèle) est effectuée, ainsi qu’une étape de sélection ascendante (où on cherche à ajouter une variable dans le modèle). La procédure s’arrête lorsque l’ajout et la suppression d’une variable dans le modèle ont des conséquences néfastes sur le modèle au vu du critère utilisé, ici encore le critère AIC (par exemple). |Rôle|Commande R| Exemple| |---|---|---| |sélection de variables pas à pas|`step()`|[ici](#step)| :Exemple{#step, toggle=popup} Séléctionner le modèle le plus adapté à la description de la variable `maxO3` par les variables quantitatives du jeu de données ozone. :Rappel : faire la régression {#doreg, toggle=collapse} ```r ozone <- subset(ozone, select = -c(vent, pluie)) regmult <- lm(maxO3 ~ ., data = ozone) ``` ```r step(regmult) ``` ``` ## Start: AIC=608.22 ## maxO3 ~ X + T9 + T12 + T15 + Ne9 + Ne12 + Ne15 + Vx9 + Vx12 + ## Vx15 + maxO3v ## ## Df Sum of Sq RSS AIC ## - Vx12 1 5.3 20641 606.25 ## - Ne12 1 5.4 20641 606.25 ## - Ne15 1 16.2 20652 606.31 ## - T9 1 25.9 20662 606.37 ## - T15 1 76.1 20712 606.64 ## - Vx15 1 78.3 20714 606.65 ## - Vx9 1 108.6 20745 606.81 ## - X 1 191.3 20827 607.26 ## <none> 20636 608.22 ## - T12 1 573.6 21210 609.30 ## - Ne9 1 1146.4 21782 612.28 ## - maxO3v 1 6170.9 26807 635.53 ## ## Step: AIC=606.25 ## maxO3 ~ X + T9 + T12 + T15 + Ne9 + Ne12 + Ne15 + Vx9 + Vx15 + ## maxO3v ## ## Df Sum of Sq RSS AIC ## - Ne12 1 8.4 20650 604.30 ## - Ne15 1 15.8 20657 604.34 ## - T9 1 21.7 20663 604.37 ## - T15 1 77.0 20718 604.67 ## - Vx15 1 163.9 20805 605.14 ## - Vx9 1 167.5 20809 605.16 ## - X 1 186.1 20827 605.26 ## <none> 20641 606.25 ## - T12 1 577.1 21218 607.34 ## - Ne9 1 1142.8 21784 610.29 ## - maxO3v 1 6173.5 26815 633.56 ## ## Step: AIC=604.3 ## maxO3 ~ X + T9 + T12 + T15 + Ne9 + Ne15 + Vx9 + Vx15 + maxO3v ## ## Df Sum of Sq RSS AIC ## - Ne15 1 8.5 20658 602.34 ## - T9 1 32.8 20683 602.48 ## - T15 1 69.4 20719 602.67 ## - Vx9 1 164.9 20815 603.19 ## - Vx15 1 166.0 20816 603.20 ## - X 1 198.9 20849 603.37 ## <none> 20650 604.30 ## - T12 1 824.3 21474 606.68 ## - Ne9 1 1872.4 22522 612.02 ## - maxO3v 1 6239.3 26889 631.87 ## ## Step: AIC=602.34 ## maxO3 ~ X + T9 + T12 + T15 + Ne9 + Vx9 + Vx15 + maxO3v ## ## Df Sum of Sq RSS AIC ## - T9 1 26.0 20684 600.49 ## - T15 1 63.6 20722 600.69 ## - Vx15 1 157.8 20816 601.20 ## - Vx9 1 177.2 20835 601.30 ## - X 1 190.5 20849 601.37 ## <none> 20658 602.34 ## - T12 1 905.3 21564 605.15 ## - Ne9 1 1973.3 22631 610.56 ## - maxO3v 1 6241.5 26900 629.91 ## ## Step: AIC=600.49 ## maxO3 ~ X + T12 + T15 + Ne9 + Vx9 + Vx15 + maxO3v ## ## Df Sum of Sq RSS AIC ## - T15 1 55.1 20739 598.78 ## - Vx15 1 142.8 20827 599.26 ## - X 1 166.4 20851 599.38 ## - Vx9 1 266.6 20951 599.92 ## <none> 20684 600.49 ## - T12 1 1006.8 21691 603.81 ## - Ne9 1 1976.6 22661 608.71 ## - maxO3v 1 6369.3 27053 628.55 ## ## Step: AIC=598.78 ## maxO3 ~ X + T12 + Ne9 + Vx9 + Vx15 + maxO3v ## ## Df Sum of Sq RSS AIC ## - Vx15 1 135.5 20875 597.51 ## - X 1 156.9 20896 597.63 ## - Vx9 1 296.5 21036 598.37 ## <none> 20739 598.78 ## - Ne9 1 2070.2 22809 607.44 ## - maxO3v 1 6639.0 27378 627.89 ## - T12 1 6772.6 27512 628.43 ## ## Step: AIC=597.51 ## maxO3 ~ X + T12 + Ne9 + Vx9 + maxO3v ## ## Df Sum of Sq RSS AIC ## - X 1 95.5 20970 596.02 ## <none> 20875 597.51 ## - Vx9 1 919.1 21794 600.34 ## - Ne9 1 2669.5 23544 608.99 ## - T12 1 6637.6 27512 626.44 ## - maxO3v 1 6732.7 27607 626.82 ## ## Step: AIC=596.02 ## maxO3 ~ T12 + Ne9 + Vx9 + maxO3v ## ## Df Sum of Sq RSS AIC ## <none> 20970 596.02 ## - Vx9 1 903.4 21874 598.75 ## - Ne9 1 2714.8 23685 607.66 ## - T12 1 6650.4 27621 624.88 ## - maxO3v 1 7363.5 28334 627.73 ``` ``` ## ## Call: ## lm(formula = maxO3 ~ T12 + Ne9 + Vx9 + maxO3v, data = ozone) ## ## Coefficients: ## (Intercept) T12 Ne9 Vx9 maxO3v ## 12.6313 2.7641 -2.5154 1.2929 0.3548 ``` Tous les modèles intermédiaires des différentes étapes de la procédure stepwise sont conservés. Le modèle final ne contient plus qu’une variable de température, une de vent, une de nébulosité ainsi que maxO3v. Ainsi, la procédure a sélectionné une variable par groupe de variables corrélées entre elles. :Suite Régression linéaire multiple {#regmul, toggle=collapse, title-display=hidden} [Approche descriptive préliminaire](caps_9_1_regression_lineaire_multiple_description.html) [Mise en oeuvre](caps_9_2_regression_lineaire_multiple_mise_en_oeuvre.html) [Affiner le modèle](caps_9_3_regression_lineaire_multiple_affiner.html)
