{"id":6447,"date":"2026-10-14T18:02:07","date_gmt":"2026-10-14T18:02:07","guid":{"rendered":"http:\/\/french.arabpsychology.com\/statistics\/comment-tester-la-normalite-dans-stata\/"},"modified":"2026-10-14T18:02:07","modified_gmt":"2026-10-14T18:02:07","slug":"comment-tester-la-normalite-dans-stata","status":"publish","type":"psych_stat","link":"https:\/\/french.arabpsychology.com\/?psych_stat=comment-tester-la-normalite-dans-stata","title":{"rendered":"Comment tester la normalit\u00e9 dans Stata"},"content":{"rendered":"<p>L&rsquo;analyse statistique appliqu\u00e9e repose sur un ensemble d&rsquo;axiomes fondamentaux dont la v\u00e9rification conditionne la robustesse, la validit\u00e9 interne et la port\u00e9e inf\u00e9rentielle des conclusions scientifiques. Parmi ces postulats, l&rsquo;hypoth\u00e8se selon laquelle les variables al\u00e9atoires observ\u00e9es \u2014 ou plus pr\u00e9cis\u00e9ment les perturbations stochastiques sous-jacentes \u2014 suivent une loi normale occupe une position pr\u00e9\u00e9minente dans la mod\u00e9lisation param\u00e9trique. Que l&rsquo;on mobilise l&rsquo;\u00e9conom\u00e9trie, la biom\u00e9trie, la psychom\u00e9trie ou la sociologie quantitative, la loi de Gauss constitue le socle th\u00e9orique de l&rsquo;estimation par les moindres carr\u00e9s ordinaires, de l&rsquo;analyse de la variance (<a rel=\"noopener noreferrer\" href=\"https:\/\/www.statmethods.net\/stats\/anova.html\">ANOVA<\/a>) et des m\u00e9thodes multivari\u00e9es canoniques. N\u00e9anmoins, l&rsquo;acceptation dogmatique de cette r\u00e9gularit\u00e9 distributionnelle expose le chercheur \u00e0 des distorsions diagnostiques s\u00e9v\u00e8res, pouvant compromettre la significativit\u00e9 statistique des tests d&rsquo;hypoth\u00e8ses et invalider les pr\u00e9dictions d&rsquo;un mod\u00e8le.<\/p>\n<p>Dans l&rsquo;environnement de calcul statistique <a rel=\"noopener noreferrer\" href=\"https:\/\/www.stata.com\/\">Stata<\/a>, l&rsquo;\u00e9valuation de la normalit\u00e9 ne saurait se r\u00e9duire \u00e0 une routine m\u00e9canique ou \u00e0 l&rsquo;application aveugle d&rsquo;un indicateur unique. Le logiciel fournit un \u00e9ventail complet d&rsquo;outils analytiques allant de repr\u00e9sentations graphiques fines \u00e0 des tests formels d&rsquo;ad\u00e9quation bas\u00e9s sur l&rsquo;alignement des quantiles, les d\u00e9compositions de moments ou la distance empirique globale. Comprendre les fondements th\u00e9oriques, les \u00e9quations sous-jacentes, les propri\u00e9t\u00e9s d&rsquo;\u00e9chantillonnage et les sensibilit\u00e9s relatives de ces proc\u00e9dures s&rsquo;av\u00e8re indispensable pour tout praticien soucieux d&rsquo;int\u00e9grit\u00e9 m\u00e9thodologique. Une telle d\u00e9marche permet de d\u00e9passer l&rsquo;opposition simpliste entre rejet et acceptation d&rsquo;une distribution pour adopter une approche diagnostique int\u00e9gr\u00e9e et gradu\u00e9e.<\/p>\n<p>Ce trait\u00e9 m\u00e9thodologique propose une exploration exhaustive de l&rsquo;ensemble des techniques diagnostiques de la normalit\u00e9 univari\u00e9e, r\u00e9siduelle et multivari\u00e9e impl\u00e9ment\u00e9es sous Stata. \u00c0 travers un expos\u00e9 m\u00e9thodique associant la d\u00e9rivation math\u00e9matique des statistiques de test, l&rsquo;explication minutieuse de la syntaxe de programmation, l&rsquo;analyse d\u00e9taill\u00e9e des sorties num\u00e9riques et l&rsquo;interpr\u00e9tation des structures graphiques, ce document fournit aux analystes de donn\u00e9es, doctorants et chercheurs chevronn\u00e9s une feuille de route rigoureuse pour auditer la conformit\u00e9 de leurs donn\u00e9es empiriques aux exigences distributionnelles de la th\u00e9orie statistique classique.<\/p>\n<h2>1. Introduction aux hypoth\u00e8ses distributionnelles et au postulat de normalit\u00e9 dans Stata<\/h2>\n<h3>1.1 Importance de la distribution gaussienne en recherche empirique<\/h3>\n<p>La primaut\u00e9 de la distribution normale dans l&rsquo;inf\u00e9rence statistique contemporaine d\u00e9coule de sa structure math\u00e9matique \u00e9l\u00e9gante, caract\u00e9ris\u00e9e par une sym\u00e9trie parfaite autour de son esp\u00e9rance math\u00e9matique et une d\u00e9croissance exponentielle de ses densit\u00e9s marginales. Dans le cadre de l&rsquo;inf\u00e9rence param\u00e9trique classique, la d\u00e9rivation exacte des distributions d&rsquo;\u00e9chantillonnage des estimateurs \u2014 notamment les lois de Student pour les coefficients de r\u00e9gression et les lois de Fisher-Snedecor pour les contrastes lin\u00e9aires simultan\u00e9s \u2014 d\u00e9pend de l&rsquo;hypoth\u00e8se de normalit\u00e9 des termes d&rsquo;erreur. Lorsque cette condition d&rsquo;ellipticit\u00e9 est enfreinte dans des \u00e9chantillons de dimension restreinte, la distribution empirique des statistiques de test s&rsquo;\u00e9carte substantiellement de leurs lois th\u00e9oriques tabul\u00e9es. Cela entra\u00eene une alt\u00e9ration de l&rsquo;erreur de premi\u00e8re esp\u00e8ce, augmentant le risque de faux positifs, ainsi qu&rsquo;une d\u00e9t\u00e9rioration de la puissance statistique des estimateurs.<\/p>\n<p>Le th\u00e9or\u00e8me central limite (<a rel=\"noopener noreferrer\" href=\"https:\/\/mathworld.wolfram.com\/CentralLimitTheorem.html\">Central Limit Theorem<\/a>) est fr\u00e9quemment invoqu\u00e9 dans la pratique \u00e9conom\u00e9trique pour att\u00e9nuer la port\u00e9e des non-normalit\u00e9s. Selon ce principe fondamental formul\u00e9 sous les conditions de Lindeberg-L\u00e9vy ou de Liapounov, la somme de variables al\u00e9atoires ind\u00e9pendantes et identiquement distribu\u00e9es, dot\u00e9es d&rsquo;une variance finie, converge asymptotiquement vers une distribution gaussienne standardis\u00e9e, quelle que soit la loi m\u00e8re sous-jacente. Si ce m\u00e9canisme de convergence justifie l&rsquo;application asymptotique des tests de Wald ou du multiplicateur de Lagrange sur des macro-\u00e9chantillons, il ne saurait exon\u00e9rer l&rsquo;analyste d&rsquo;une \u00e9valuation rigoureuse sur des cohortes limit\u00e9es ou interm\u00e9diaires. Dans les \u00e9chantillons finis, les \u00e9carts marqu\u00e9s par rapport \u00e0 la normalit\u00e9 \u2014 en particulier les distributions bimodalit\u00e9s prononc\u00e9es, les asym\u00e9tries fortes ou les queues de distribution \u00e9paisses \u2014 d\u00e9gradent l&rsquo;efficacit\u00e9 relative des estimateurs lin\u00e9aires sans biais et rendent impr\u00e9cis le calcul des intervalles de confiance.<\/p>\n<p>Au-del\u00e0 de la simple inf\u00e9rence sur les moyennes conditionnelles, le maintien du postulat de normalit\u00e9 s&rsquo;av\u00e8re critique lors de la construction d&rsquo;intervalles de pr\u00e9diction individuels. Contrairement aux intervalles de confiance encadrant l&rsquo;esp\u00e9rance de la variable d\u00e9pendante, dont la pr\u00e9cision b\u00e9n\u00e9ficie directement de la loi des grands nombres, la d\u00e9limitation d&rsquo;un intervalle de pr\u00e9vision pour une observation singuli\u00e8re incorpore directement la variance du terme r\u00e9siduel individuel. Si la distribution des erreurs n&rsquo;est pas strictement normale, les bornes de pr\u00e9diction calcul\u00e9es sur la base de la loi de Student standard sous-estimeront ou surestimeront la probabilit\u00e9 r\u00e9elle de couverture, induisant des erreurs de ciblage d\u00e9cisionnel dans les applications cliniques, actuarielles ou \u00e9conom\u00e9triques.<\/p>\n<h3>1.2 Distinction entre normalit\u00e9 univari\u00e9e, r\u00e9siduelle et multivari\u00e9e<\/h3>\n<p>Une confusion m\u00e9thodologique r\u00e9currente consiste \u00e0 exiger la normalit\u00e9 marginale des variables d\u00e9pendantes ou ind\u00e9pendantes observ\u00e9es avant toute d\u00e9marche de mod\u00e9lisation. En r\u00e9alit\u00e9, le cadre du mod\u00e8le lin\u00e9aire g\u00e9n\u00e9ralis\u00e9 n&rsquo;impose aucune contrainte de normalit\u00e9 sur la distribution brute des r\u00e9gresseurs, ni m\u00eame sur la distribution marginale inconditionnelle de la variable expliqu\u00e9e. Les pr\u00e9dicteurs peuvent \u00eatre binaires, nominaux, polytomiques ou asym\u00e9triques sans enfreindre les propri\u00e9t\u00e9s des moindres carr\u00e9s ordinaires. Ce qui importe pour la validit\u00e9 du mod\u00e8le param\u00e9trique n&rsquo;est pas la normalit\u00e9 marginale d&rsquo;une variable prise isol\u00e9ment, mais la conformit\u00e9 de la distribution conditionnelle de la variable d\u00e9pendante sachant les pr\u00e9dicteurs, ce qui \u00e9quivaut \u00e0 exiger la normalit\u00e9 de la composante stochastique r\u00e9siduelle non capt\u00e9e par le mod\u00e8le.<\/p>\n<p>La normalit\u00e9 r\u00e9siduelle renvoie sp\u00e9cifiquement \u00e0 la distribution du vecteur des perturbations inobservables dans la population parente. Dans une r\u00e9gression lin\u00e9aire multiple, ces perturbations sont suppos\u00e9es \u00eatre g\u00e9n\u00e9r\u00e9es par un processus stochastique sph\u00e9rique, caract\u00e9ris\u00e9 par une esp\u00e9rance math\u00e9matique conditionnelle nulle, une homosc\u00e9dasticit\u00e9 et une ind\u00e9pendance s\u00e9rielle, le tout sous-tendu par une distribution normale. L&rsquo;audit de normalit\u00e9 doit par cons\u00e9quent s&rsquo;orienter en priorit\u00e9 vers les r\u00e9sidus estim\u00e9s apr\u00e8s convergence de l&rsquo;algorithme d&rsquo;optimisation, plut\u00f4t que vers la variable d\u00e9pendante brute dont l&rsquo;allure multimodale peut simplement refl\u00e9ter l&rsquo;effet de stratification de covariables discr\u00e8tes sous-jacentes.<\/p>\n<p>Enfin, la normalit\u00e9 multivari\u00e9e pose des exigences conceptuelles beaucoup plus restrictives que la juxtaposition de normalit\u00e9s marginales univari\u00e9es. Un vecteur al\u00e9atoire multidimensionnel suit une loi normale multivari\u00e9e si et seulement si toute combinaison lin\u00e9aire scalaire de ses composantes suit une loi normale univari\u00e9e. Il est math\u00e9matiquement d\u00e9montr\u00e9 que l&rsquo;existence de distributions marginales univari\u00e9es parfaitement gaussiennes pour plusieurs variables distinctes ne garantit en rien que leur distribution conjointe respecte la normalit\u00e9 multivari\u00e9e. Cette distinction est cruciale dans les architectures analytiques fond\u00e9es sur l&rsquo;estimation du maximum de vraisemblance \u00e0 information compl\u00e8te, telles que l&rsquo;analyse factorielle confirmatoire, la mod\u00e9lisation par \u00e9quations structurelles (<a rel=\"noopener noreferrer\" href=\"https:\/\/en.wikipedia.org\/wiki\/Structural_equation_modeling\">SEM<\/a>) ou l&rsquo;analyse discriminante lin\u00e9aire de Fisher.<\/p>\n<h3>1.3 Pr\u00e9sentation du jeu de donn\u00e9es de travail et pr\u00e9paration de l&rsquo;environnement<\/h3>\n<p>Afin de concr\u00e9tiser l&rsquo;ensemble des d\u00e9marches diagnostiques pr\u00e9sent\u00e9es au sein de ce trait\u00e9, nous mobiliserons le jeu de donn\u00e9es classique int\u00e9gr\u00e9 \u00e0 l&rsquo;environnement logiciel Stata, accessible via l&rsquo;archive interne <code>sysuse auto.dta<\/code>. Ce corpus recueille les caract\u00e9ristiques techniques, dimensionnelles et \u00e9conomiques d&rsquo;une flotte de 74 mod\u00e8les d&rsquo;automobiles commercialis\u00e9s en 1978. La composition de ce fichier combine des m\u00e9triques continues affichant des niveaux variables d&rsquo;asym\u00e9trie et d&rsquo;aplatissement, ce qui permet d&rsquo;illustrer la r\u00e9ponse empirique des diff\u00e9rents algorithmes de test face \u00e0 des configurations distributionnelles diversifi\u00e9es.<\/p>\n<p>L&rsquo;initialisation d&rsquo;un environnement d&rsquo;analyse rigoureux dans Stata n\u00e9cessite la configuration pr\u00e9alable d&rsquo;un fichier journal permettant la tra\u00e7abilit\u00e9 int\u00e9grale des manipulations et l&rsquo;ex\u00e9cution d&rsquo;une inspection pr\u00e9liminaire des observations manquantes. Les variables principales qui feront l&rsquo;objet de nos exp\u00e9rimentations sont le prix (<code>price<\/code>), la consommation de carburant exprim\u00e9e en miles par gallon (<code>mpg<\/code>), le poids du v\u00e9hicule (<code>weight<\/code>) et la cylindr\u00e9e du moteur (<code>displacement<\/code>). Les commandes fondamentales d&rsquo;amorce se structurent comme suit :<\/p>\n<p><code>sysuse auto, clear<\/code><br \/>\n<code>describe price mpg weight displacement<\/code><br \/>\n<code>summarize price mpg weight displacement<\/code><br \/>\n<code>misstable summarize price mpg weight displacement<\/code><\/p>\n<p>L&rsquo;instruction <code>misstable summarize<\/code> garantit l&rsquo;exhaustivit\u00e9 des sous-\u00e9chantillons d&rsquo;analyse en r\u00e9pertoriant les valeurs manquantes structurelles. La pr\u00e9sence de valeurs absentes non al\u00e9atoires peut biaiser les distributions marginales en \u00e9liminant s\u00e9lectivement les extr\u00eames. Dans le jeu de donn\u00e9es <code>auto.dta<\/code>, l&rsquo;int\u00e9grit\u00e9 des 74 observations est pr\u00e9serv\u00e9e pour ces vecteurs continus, fournissant un cadre de calibration stable pour l&rsquo;ensemble des estimateurs univari\u00e9s et multivari\u00e9s.<\/p>\n<h2>2. \u00c9valuation graphique pr\u00e9liminaire : L&rsquo;histogramme et la superposition de densit\u00e9 normale<\/h2>\n<h3>2.1 Syntaxe fondamentale et structure de la commande histogram<\/h3>\n<p>L&rsquo;histogramme constitue la m\u00e9thode exploratoire non param\u00e9trique la plus intuitive pour visualiser la fonction de masse ou de densit\u00e9 empirique d&rsquo;un vecteur continu. Dans l&rsquo;\u00e9cosyst\u00e8me Stata, la commande <code>histogram<\/code> (ou son abr\u00e9viation autoris\u00e9e <code>hist<\/code>) partitionne le support de la variable \u00e9tudi\u00e9e en une s\u00e9rie d&rsquo;intervalles contigus, appel\u00e9s classes ou r\u00e9ceptacles statistiques, puis comptabilise la fr\u00e9quence des r\u00e9alisations observ\u00e9es au sein de chaque segment. Cette discr\u00e9tisation permet d&rsquo;identifier imm\u00e9diatement les propri\u00e9t\u00e9s qualitatives de la distribution : la localisation du mode central, l&rsquo;\u00e9ventuelle bimodalit\u00e9 r\u00e9v\u00e9latrice d&rsquo;h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9s sous-jacentes, le sens de l&rsquo;asym\u00e9trie lat\u00e9rale ainsi que l&rsquo;\u00e9paisseur apparente des queues de distribution.<\/p>\n<p>Le param\u00e9trage par d\u00e9faut de Stata calcule une ordonn\u00e9e calibr\u00e9e en densit\u00e9 statistique, assurant que la somme des aires des rectangles \u00e9l\u00e9mentaires soit rigoureusement \u00e9gale \u00e0 l&rsquo;unit\u00e9. N\u00e9anmoins, l&rsquo;analyste peut modifier l&rsquo;\u00e9chelle verticale pour adapter la restitution aux conventions de son champ disciplinaire en mobilisant des options sp\u00e9cifiques :<\/p>\n<ul>\n<li><code>density<\/code> : option implicite tra\u00e7ant la densit\u00e9 empirique, o\u00f9 la surface totale correspond \u00e0 1.<\/li>\n<li><code>fraction<\/code> : \u00e9chelle normalis\u00e9e repr\u00e9sentant la proportion relative d&rsquo;observations par classe par rapport \u00e0 l&rsquo;\u00e9chantillon global.<\/li>\n<li><code>frequency<\/code> : d\u00e9nombrement absolu des observations contenues dans chaque intervalle.<\/li>\n<li><code>percent<\/code> : conversion de la fraction en pourcentage cumulatif par classe.<\/li>\n<\/ul>\n<p>L&rsquo;ex\u00e9cution de la commande sous sa forme univari\u00e9e \u00e9l\u00e9mentaire pour la variable du prix des v\u00e9hicules se formule ainsi :<\/p>\n<p><code>histogram price, frequency title(\"Distribution des fr\u00e9quences du prix\")<\/code><\/p>\n<figure class=\"wpt-preserved-image\"><img fetchpriority=\"high\" decoding=\"async\" src=\"http:\/\/french.arabpsychology.com\/wp-content\/uploads\/2026\/09\/normalityTestStata1.png\" alt=\"Histogram example in Stata\" \/><figcaption>Histogram example in Stata<\/figcaption><\/figure>\n<p>Ce premier trac\u00e9 met en \u00e9vidence la forte d\u00e9viation de la variable <code>price<\/code> par rapport aux crit\u00e8res de sym\u00e9trie th\u00e9orique : la distribution est \u00e9tir\u00e9e vers la droite avec une concentration majeure d&rsquo;observations dans les tranches tarifaires basses et m\u00e9dianes, compl\u00e9t\u00e9e par une tra\u00eene d&rsquo;observations extr\u00eames atteignant des valeurs sup\u00e9rieures \u00e0 15 000 dollars.<\/p>\n<h3>2.2 Superposition analytique de la courbe th\u00e9orique avec l&rsquo;option normal<\/h3>\n<p>L&rsquo;appr\u00e9ciation visuelle d&rsquo;un histogramme gagne en pr\u00e9cision par la confrontation directe des fr\u00e9quences observ\u00e9es avec la fonction de densit\u00e9 gaussienne param\u00e9trique \u00e9quivalente. L&rsquo;option <code>normal<\/code> ordonne \u00e0 Stata de calculer la moyenne arithm\u00e9tique empirique et la variance corrig\u00e9e de l&rsquo;\u00e9chantillon, puis de g\u00e9n\u00e9rer une courbe de densit\u00e9 th\u00e9orique continue calibr\u00e9e sur ces moments du premier et du second ordre :<\/p>\n<p><code>histogram price, normal title(\"Histogramme de price avec densit\u00e9 normale superpos\u00e9e\")<\/code><\/p>\n<figure class=\"wpt-preserved-image\"><img loading=\"lazy\" decoding=\"async\" src=\"http:\/\/french.arabpsychology.com\/wp-content\/uploads\/2026\/09\/normalityTestStata2.png\" alt=\"Histogram with normal curve in Stata\" \/><figcaption>Histogram with normal curve in Stata<\/figcaption><\/figure>\n<p>La pr\u00e9sence de cette r\u00e9f\u00e9rence analytique r\u00e9v\u00e8le instantan\u00e9ment l&rsquo;inad\u00e9quation de la loi normale pour mod\u00e9liser cette variable brute. Alors que la courbe th\u00e9orique impose un profil en cloche rigoureusement sym\u00e9trique, les barres empiriques s&rsquo;accumulent nettement \u00e0 gauche du sommet de la cloche th\u00e9orique, tout en laissant un d\u00e9ficit syst\u00e9matique d&rsquo;observations dans la zone centrale droite et un exc\u00e9dent dans l&rsquo;extr\u00e9mit\u00e9 sup\u00e9rieure droite. Ce d\u00e9calage illustre une asym\u00e9trie positive prononc\u00e9e coupl\u00e9e \u00e0 une structure leptokurtique.<\/p>\n<p>\u00c0 l&rsquo;inverse, l&rsquo;examen de la variable relative au poids du v\u00e9hicule (<code>weight<\/code>) avec la m\u00eame option :<\/p>\n<p><code>histogram weight, normal title(\"Histogramme de weight avec densit\u00e9 normale superpos\u00e9e\")<\/code><\/p>\n<p>met en \u00e9vidence un ajustement beaucoup plus r\u00e9gulier entre les barres d&rsquo;\u00e9chantillonnage et la fonction de densit\u00e9 gaussienne. L&rsquo;absence de d\u00e9crochages aberrants sugg\u00e8re une distribution empirique bien plus proche du mod\u00e8le th\u00e9orique de Gauss, bien que seule l&rsquo;inf\u00e9rence formelle permette d&rsquo;en certifier la validit\u00e9 statistique.<\/p>\n<h3>2.3 Optimisation du d\u00e9coupage en classes et limites de l&rsquo;inspection visuelle<\/h3>\n<p>La morphologie d&rsquo;un histogramme d\u00e9pend \u00e9troitement de la strat\u00e9gie de discr\u00e9tisation adopt\u00e9e. Le nombre et la largeur des intervalles conditionnent l&rsquo;apparence du profil distributionnel : un nombre insuffisant de classes agr\u00e8ge des caract\u00e9ristiques distinctes et occulte des bimodalit\u00e9s locales, tandis qu&rsquo;un d\u00e9coupage excessif fragmente l&rsquo;\u00e9chantillon, g\u00e9n\u00e9rant des artefacts visuels caus\u00e9s par les fluctuations d&rsquo;\u00e9chantillonnage.<\/p>\n<p>Stata met en \u0153uvre par d\u00e9faut une formule de discr\u00e9tisation optimis\u00e9e inspir\u00e9e des travaux de Scott ou de Sturges, mais permet \u00e0 l&rsquo;utilisateur de prendre le contr\u00f4le param\u00e9trique complet via deux options exclusives : <code>bin(#)<\/code>, qui sp\u00e9cifie le nombre total d&rsquo;intervalles verticaux, et <code>width(#)<\/code>, qui d\u00e9finit la largeur absolue de chaque r\u00e9ceptacle sur l&rsquo;axe des abscisses. Consid\u00e9rons la sensibilit\u00e9 de l&rsquo;interpr\u00e9tation de la variable <code>mpg<\/code> face \u00e0 ces variations de granularit\u00e9 :<\/p>\n<p><code>histogram mpg, bin(5) normal title(\"Sur-lissage : 5 classes\")<\/code><br \/>\n<code>histogram mpg, bin(30) normal title(\"Sous-lissage : 30 classes\")<\/code><br \/>\n<code>histogram mpg, width(2) start(10) normal title(\"Largeur fix\u00e9e \u00e0 2 unit\u00e9s\")<\/code><\/p>\n<p>Le sur-lissage induit par <code>bin(5)<\/code> masque les micro-d\u00e9viations structurelles, donnant une illusion de normalit\u00e9. \u00c0 l&rsquo;oppos\u00e9, le d\u00e9coupage \u00e0 30 classes accentue artificiellement les discontinuit\u00e9s d&rsquo;\u00e9chantillonnage. Il est donc recommand\u00e9 d&rsquo;adopter une largeur de bande coh\u00e9rente avec la formule de Freedman-Diaconis, bas\u00e9e sur l&rsquo;\u00e9cart interquartile :<\/p>\n<p><code>h = 2 * IQR * n^(-1\/3)<\/code><\/p>\n<p>Malgr\u00e9 ces raffinements d&rsquo;\u00e9chelle, l&rsquo;\u00e9valuation par histogramme reste une m\u00e9thode purement qualitative, hautement d\u00e9pendante de la subjectivit\u00e9 de l&rsquo;observateur. Elle sert d&rsquo;\u00e9tape exploratoire pr\u00e9liminaire, mais doit imp\u00e9rativement \u00eatre compl\u00e9t\u00e9e par des repr\u00e9sentations comparatives de rang et des tests inf\u00e9rentiels rigoureux.<\/p>\n<h2>3. Inspection visuelle avanc\u00e9e : Les diagrammes quantile-quantile et probabilit\u00e9-probabilit\u00e9<\/h2>\n<h3>3.1 Principes g\u00e9om\u00e9triques et interpr\u00e9tation du diagramme Q-Q via qnorm<\/h3>\n<p>Le diagramme quantile-quantile normal, impl\u00e9ment\u00e9 dans Stata via la commande <code>qnorm<\/code>, offre une r\u00e9solution graphique bien sup\u00e9rieure \u00e0 celle de l&rsquo;histogramme pour identifier les anomalies aux extr\u00e9mit\u00e9s d&rsquo;une distribution. Le principe sous-jacent repose sur la confrontation directe entre les quantiles empiriques observ\u00e9s dans l&rsquo;\u00e9chantillon et les quantiles th\u00e9oriques issus d&rsquo;une distribution normale standardis\u00e9e dot\u00e9e de param\u00e8tres de moyenne et de dispersion identiques.<\/p>\n<p>Sur le plan g\u00e9om\u00e9trique, si le vecteur \u00e9tudi\u00e9 d\u00e9coule fid\u00e8lement d&rsquo;un processus g\u00e9n\u00e9rateur gaussien, les points de coordonn\u00e9es repr\u00e9sentant chaque quantile se superposent le long d&rsquo;une droite diagonale de r\u00e9f\u00e9rence \u00e0 45 degr\u00e9s. Tout \u00e9cart syst\u00e9matique par rapport \u00e0 cette trajectoire lin\u00e9aire traduit une non-conformit\u00e9 distributionnelle sp\u00e9cifique :<\/p>\n<ul>\n<li><strong>Queues lourdes (leptokurticit\u00e9) :<\/strong> Les quantiles observ\u00e9s inf\u00e9rieurs se placent sous la droite de r\u00e9f\u00e9rence \u00e0 gauche, tandis que les quantiles sup\u00e9rieurs se placent au-dessus de la droite \u00e0 droite, formant une courbure caract\u00e9ristique en \u00ab S \u00bb invers\u00e9.<\/li>\n<li><strong>Queues l\u00e9g\u00e8res (platykurticit\u00e9) :<\/strong> La trajectoire des quantiles croise la droite diagonale en sens oppos\u00e9, g\u00e9n\u00e9rant une ondulation en \u00ab S \u00bb.<\/li>\n<li><strong>Asym\u00e9trie positive (\u00e9talement \u00e0 droite) :<\/strong> La trajectoire pr\u00e9sente une courbure globalement convexe, restant sous la diagonale dans la r\u00e9gion centrale avant de diverger fortement vers le haut pour les quantiles sup\u00e9rieurs.<\/li>\n<li><strong>Asym\u00e9trie n\u00e9gative (\u00e9talement \u00e0 gauche) :<\/strong> La s\u00e9rie de points dessine une courbure concave au-dessus de la droite de r\u00e9f\u00e9rence sur la majorit\u00e9 du support.<\/li>\n<\/ul>\n<p>L&rsquo;ex\u00e9cution de la proc\u00e9dure pour la variable <code>price<\/code> met en \u00e9vidence cette convexit\u00e9 :<\/p>\n<p><code>qnorm price, title(\"Diagramme Q-Q de la variable Price\")<\/code><\/p>\n<p>Le trac\u00e9 g\u00e9n\u00e9r\u00e9 montre un alignement acceptable pour les valeurs centrales modestes, suivi d&rsquo;une divergence exponentielle au-del\u00e0 du quantile th\u00e9orique d&rsquo;ordre 1, ce qui confirme que la queue sup\u00e9rieure de la distribution des prix s&rsquo;\u00e9tend bien au-del\u00e0 des bornes pr\u00e9vues par la th\u00e9orie normale.<\/p>\n<h3>3.2 \u00c9valuation de la r\u00e9gion centrale par le diagramme P-P avec pnorm<\/h3>\n<p>Alors que le diagramme quantile-quantile accentue le comportement des valeurs extr\u00eames situ\u00e9es dans les queues, le diagramme probabilit\u00e9-probabilit\u00e9, pilot\u00e9 par l&rsquo;instruction <code>pnorm<\/code>, concentre sa sensibilit\u00e9 analytique sur la zone centrale de la distribution. Cette technique juxtapose la fonction de r\u00e9partition empirique standardis\u00e9e d&rsquo;un \u00e9chantillon \u00e0 la fonction de r\u00e9partition th\u00e9orique cumulative de la loi normale.<\/p>\n<p>Le graphique repr\u00e9sente en ordonn\u00e9e la probabilit\u00e9 empirique cumul\u00e9e associ\u00e9e \u00e0 chaque observation ordonn\u00e9e, et en abscisse la probabilit\u00e9 cumulative attendue sous l&rsquo;hypoth\u00e8se nulle de normalit\u00e9 :<\/p>\n<p><code>pnorm price, title(\"Diagramme P-P de la variable Price\")<\/code><\/p>\n<p>Dans ce plan unitaire born\u00e9 sur l&rsquo;intervalle [0, 1], les points extr\u00eames convergent m\u00e9caniquement vers les coordonn\u00e9es (0,0) et (1,1). Les d\u00e9viations les plus prononc\u00e9es s&rsquo;observent au centre de la distribution :<\/p>\n<ul>\n<li>Une trajectoire s&rsquo;\u00e9levant nettement au-dessus de la bissectrice au centre indique que les probabilit\u00e9s cumul\u00e9es observ\u00e9es progressent plus vite que pr\u00e9vu par le mod\u00e8le gaussien, r\u00e9v\u00e9lant une forte concentration de donn\u00e9es sur des valeurs inf\u00e9rieures \u00e0 la moyenne.<\/li>\n<li>Une trajectoire sinuso\u00efdale coupant la diagonale r\u00e9v\u00e8le des anomalies d&rsquo;aplatissement ou des distributions multimodales au c\u0153ur de la distribution.<\/li>\n<\/ul>\n<p>La compl\u00e9mentarit\u00e9 entre <code>qnorm<\/code> et <code>pnorm<\/code> constitue une m\u00e9thode diagnostique visuelle robuste : le diagramme Q-Q s\u00e9curise l&rsquo;analyste contre le risque d&rsquo;ignorer des valeurs aberrantes ou des queues \u00e9paisses influentes, tandis que le diagramme P-P permet d&rsquo;\u00e9valuer la r\u00e9gularit\u00e9 et la sym\u00e9trie de la masse centrale de l&rsquo;\u00e9chantillon.<\/p>\n<h3>3.3 Personnalisation des trac\u00e9s graphiques pour les publications acad\u00e9miques<\/h3>\n<p>L&rsquo;int\u00e9gration de diagnostics distributionnels dans des manuscrits soumis \u00e0 \u00e9valuation par les pairs n\u00e9cessite une mise en forme graphique soign\u00e9e. Stata dispose d&rsquo;options graphiques avanc\u00e9es permettant d&rsquo;ajuster les contrastes visuels, d&rsquo;appliquer des chartes graphiques adapt\u00e9es \u00e0 l&rsquo;impression monochrome et de juxtaposer des diagnostics compl\u00e9mentaires au sein d&rsquo;une m\u00eame figure exportable en haute r\u00e9solution vectorielle.<\/p>\n<p>L&rsquo;encha\u00eenement m\u00e9thodique ci-dessous combine l&rsquo;histogramme, le diagramme Q-Q et le diagramme P-P au sein d&rsquo;un canevas unifi\u00e9, en s&rsquo;appuyant sur le sch\u00e9ma de couleurs <code>cleanplots<\/code> ou le th\u00e8me institutionnel classique <code>sj<\/code> (Stata Journal) :<\/p>\n<p><code>set scheme s2mono<\/code><br \/>\n<code>histogram price, normal frequency title(\"A: Histogramme &amp; Densit\u00e9\", size(medium)) name(g1, replace)<\/code><br \/>\n<code>qnorm price, grid title(\"B: Alignement des Quantiles\", size(medium)) name(g2, replace)<\/code><br \/>\n<code>pnorm price, grid title(\"C: Probabilit\u00e9s Cumul\u00e9es\", size(medium)) name(g3, replace)<\/code><br \/>\n<code>graph combine g1 g2 g3, col(3) iscale(0.8) title(\"Audit Distributionnel Global : Variable Price\", size(large))<\/code><br \/>\n<code>graph export \"diagnostic_normalite_price.eps\", as(eps) preview(on) replace<\/code><\/p>\n<p>L&rsquo;utilisation de formats vectoriels tels que l&rsquo;EPS (<a rel=\"noopener noreferrer\" href=\"https:\/\/en.wikipedia.org\/wiki\/Encapsulated_PostScript\">Encapsulated PostScript<\/a>) ou le PDF pr\u00e9serve la nettet\u00e9 des courbes et des points de dispersion lors des r\u00e9ductions d&rsquo;\u00e9chelle impos\u00e9es par les maquettes \u00e9ditoriales. Cette combinaison graphique offre aux \u00e9valuateurs une vue d&rsquo;ensemble du profil distributionnel de la variable \u00e9tudi\u00e9e.<\/p>\n<h2>4. Le test de normalit\u00e9 de Shapiro-Wilk sous Stata<\/h2>\n<h3>4.1 Fondements m\u00e9thodologiques et propri\u00e9t\u00e9s du coefficient W<\/h3>\n<p>Le test formalis\u00e9 par <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1093\/biomet\/52.3-4.591\">Samuel Sanford Shapiro et Martin Wilk en 1965<\/a> est largement consid\u00e9r\u00e9 comme l&rsquo;un des tests d&rsquo;ad\u00e9quation univari\u00e9s les plus puissants pour d\u00e9tecter les d\u00e9viations par rapport \u00e0 la loi normale. Le principe du test repose sur l&rsquo;\u00e9valuation de la corr\u00e9lation lin\u00e9aire entre les statistiques d&rsquo;ordre d&rsquo;un \u00e9chantillon observ\u00e9 et les esp\u00e9rances math\u00e9matiques des statistiques d&rsquo;ordre d&rsquo;un \u00e9chantillon th\u00e9orique de m\u00eame taille issu d&rsquo;une loi normale standardis\u00e9e.<\/p>\n<p>La statistique de test, d\u00e9sign\u00e9e par la notation standard <span class=\"arabpsy-math-inline\" data-latex=\"W\">$W$<\/span>, s&rsquo;exprime comme le ratio de deux estimateurs distincts de la variance de la population :<\/p>\n<p><code>W = (sum(a_i * x_(i)))^2 \/ sum((x_i - mean(x))^2)<\/code><\/p>\n<p>Dans cette formulation math\u00e9matique, <span class=\"arabpsy-math-inline\" data-latex=\"x_{(i)}\">$x_{(i)}$<\/span> d\u00e9signe la <span class=\"arabpsy-math-inline\" data-latex=\"i\">$i$<\/span>-\u00e8me observation ordonn\u00e9e de l&rsquo;\u00e9chantillon, <span class=\"arabpsy-math-inline\" data-latex=\"\\bar{x}\">$\\bar{x}$<\/span> correspond \u00e0 la moyenne empirique, et le vecteur des poids <span class=\"arabpsy-math-inline\" data-latex=\"a = (a_1, ..., a_n)&#039;\">$a = (a_1, &#8230;, a_n)&rsquo;$<\/span> est d\u00e9riv\u00e9 de la matrice de variance-covariance des statistiques d&rsquo;ordre d&rsquo;une distribution normale univari\u00e9e selon l&rsquo;\u00e9quation :<\/p>\n<p><code>a' = (m' * V^(-1)) \/ (m' * V^(-1) * V^(-1) * m)^(1\/2)<\/code><\/p>\n<p>o\u00f9 <span class=\"arabpsy-math-inline\" data-latex=\"m\">$m$<\/span> repr\u00e9sente le vecteur des esp\u00e9rances des statistiques d&rsquo;ordre gaussiennes et <span class=\"arabpsy-math-inline\" data-latex=\"V\">$V$<\/span> symbolise leur matrice de covariance. Le coefficient <span class=\"arabpsy-math-inline\" data-latex=\"W\">$W$<\/span> prend obligatoirement des valeurs contenues dans l&rsquo;intervalle semi-ouvert ]0, 1]. Une valeur de <span class=\"arabpsy-math-inline\" data-latex=\"W\">$W$<\/span> tr\u00e8s proche de l&rsquo;unit\u00e9 indique un alignement presque parfait entre les quantiles empiriques et th\u00e9oriques, apportant un soutien empirique \u00e0 l&rsquo;hypoth\u00e8se nulle <span class=\"arabpsy-math-inline\" data-latex=\"H_0\">$H_0$<\/span> de normalit\u00e9. \u00c0 l&rsquo;inverse, des valeurs de <span class=\"arabpsy-math-inline\" data-latex=\"W\">$W$<\/span> significativement inf\u00e9rieures \u00e0 1 indiquent une d\u00e9formation distributionnelle, conduisant au rejet de <span class=\"arabpsy-math-inline\" data-latex=\"H_0\">$H_0$<\/span>. L&rsquo;algorithme d\u00e9velopp\u00e9 par Patrick Royston a permis d&rsquo;\u00e9tendre la port\u00e9e op\u00e9rationnelle du test jusqu&rsquo;\u00e0 des \u00e9chantillons atteignant 2 000 unit\u00e9s, tout en fournissant une transformation normalis\u00e9e facilitant le calcul exact des p-valeurs associ\u00e9es.<\/p>\n<h3>4.2 Syntaxe op\u00e9rationnelle et lecture des sorties de la commande swilk<\/h3>\n<p>L&rsquo;impl\u00e9mentation de la proc\u00e9dure sous Stata s&rsquo;effectue via l&rsquo;appel direct de la commande <code>swilk<\/code>, suivie du nom de la variable inspect\u00e9e ou d&rsquo;une liste exhaustive de variables. Appliquons ce test \u00e0 notre jeu de donn\u00e9es :<\/p>\n<p><code>swilk price<\/code><\/p>\n<figure class=\"wpt-preserved-image\"><img loading=\"lazy\" decoding=\"async\" src=\"http:\/\/french.arabpsychology.com\/wp-content\/uploads\/2026\/09\/normalityTestStata3.png\" alt=\"Shapiro Wilk Test output in Stata\" \/><figcaption>Shapiro Wilk Test output in Stata<\/figcaption><\/figure>\n<p>La sortie standard g\u00e9n\u00e9r\u00e9e par Stata s&rsquo;articule autour d&rsquo;un tableau \u00e0 colonnes multiples contenant des m\u00e9triques fondamentales :<\/p>\n<ul>\n<li><code>Obs<\/code> : Le volume d&rsquo;observations valides prises en compte dans le calcul (ici, <span class=\"arabpsy-math-inline\" data-latex=\"N = 74\">$N = 74$<\/span>).<\/li>\n<li><code>W<\/code> : La statistique de Shapiro-Wilk calcul\u00e9e (<span class=\"arabpsy-math-inline\" data-latex=\"W = 0.758\">$W = 0.758$<\/span> pour <code>price<\/code>).<\/li>\n<li><code>V<\/code> : L&rsquo;indice standardis\u00e9 d\u00e9velopp\u00e9 par Royston, repr\u00e9sentant le ratio entre la variance attendue sous <span class=\"arabpsy-math-inline\" data-latex=\"H_0\">$H_0$<\/span> et la variance observ\u00e9e. Une valeur de <span class=\"arabpsy-math-inline\" data-latex=\"V = 1\">$V = 1$<\/span> correspond exactement au comportement th\u00e9orique d&rsquo;une distribution normale. Des valeurs de <span class=\"arabpsy-math-inline\" data-latex=\"V\">$V$<\/span> significativement sup\u00e9rieures \u00e0 1 signalent un \u00e9loignement marqu\u00e9 par rapport \u00e0 la normalit\u00e9.<\/li>\n<li><code>z<\/code> : Le score <span class=\"arabpsy-math-inline\" data-latex=\"Z\">$Z$<\/span> r\u00e9sultant de la transformation de la statistique <span class=\"arabpsy-math-inline\" data-latex=\"V\">$V$<\/span> en une d\u00e9viation normale centr\u00e9e r\u00e9duite standard. Cette valeur permet d&rsquo;ordonner l&rsquo;intensit\u00e9 des rejets entre plusieurs variables.<\/li>\n<li><code>Prob &gt; z<\/code> : La valeur critique (<span class=\"arabpsy-math-inline\" data-latex=\"p\">$p$<\/span>-value) associ\u00e9e au test d&rsquo;hypoth\u00e8se unilat\u00e9ral.<\/li>\n<\/ul>\n<p>Dans le cas d&rsquo;esp\u00e8ce de la variable <code>price<\/code>, la statistique <span class=\"arabpsy-math-inline\" data-latex=\"W\">$W$<\/span> s&rsquo;\u00e9tablit \u00e0 un niveau tr\u00e8s bas, g\u00e9n\u00e9rant un indice <span class=\"arabpsy-math-inline\" data-latex=\"V\">$V$<\/span> \u00e9lev\u00e9 et une <span class=\"arabpsy-math-inline\" data-latex=\"p\">$p$<\/span>-value associ\u00e9e largement inf\u00e9rieure au seuil alpha traditionnel de 0,05 (affich\u00e9e sous la forme <code>0.00000<\/code>). L&rsquo;hypoth\u00e8se nulle de normalit\u00e9 est donc rejet\u00e9e avec une marge d&rsquo;erreur infinit\u00e9simale. \u00c0 l&rsquo;inverse, l&rsquo;application de l&rsquo;instruction sur la variable <code>weight<\/code> :<\/p>\n<p><code>swilk weight<\/code><\/p>\n<p>g\u00e9n\u00e8re une valeur de <span class=\"arabpsy-math-inline\" data-latex=\"W = 0.941\">$W = 0.941$<\/span> avec une p-valeur de 0.0017, ce qui conduit \u00e9galement au rejet formel de la normalit\u00e9 au seuil de 5 %, bien que l&rsquo;amplitude de la d\u00e9viation soit nettement moins marqu\u00e9e que pour le prix. L&rsquo;analyste peut \u00e9galement ex\u00e9cuter le test simultan\u00e9ment sur une s\u00e9rie \u00e9tendue de variables continues :<\/p>\n<p><code>swilk price mpg weight displacement trunk<\/code><\/p>\n<figure class=\"wpt-preserved-image\"><img loading=\"lazy\" decoding=\"async\" src=\"http:\/\/french.arabpsychology.com\/wp-content\/uploads\/2026\/09\/normalityTestStata4.png\" alt=\"Multiple Shapiro-Wilk tests at once in Stata\" \/><figcaption>Multiple Shapiro-Wilk tests at once in Stata<\/figcaption><\/figure>\n<p>Cette approche group\u00e9e permet de comparer rapidement les degr\u00e9s relatifs de conformit\u00e9 distributionnelle au sein d&rsquo;une base de donn\u00e9es.<\/p>\n<h3>4.3 Limites inh\u00e9rentes \u00e0 la taille de l&rsquo;\u00e9chantillon avec swilk<\/h3>\n<p>Bien que le test de Shapiro-Wilk affiche d&rsquo;excellentes propri\u00e9t\u00e9s statistiques face \u00e0 des alternatives tr\u00e8s diverses, son interpr\u00e9tation pratique d\u00e9pend fortement de la taille de l&rsquo;\u00e9chantillon <span class=\"arabpsy-math-inline\" data-latex=\"N\">$N$<\/span>. Cette sensibilit\u00e9 refl\u00e8te le compromis classique entre la puissance statistique du test et l&rsquo;erreur d&rsquo;\u00e9chantillonnage.<\/p>\n<p>Dans les \u00e9chantillons de grande taille (<span class=\"arabpsy-math-inline\" data-latex=\"N &gt; 500\">$N &gt; 500$<\/span> ou $1,000$), la puissance statistique de <code>swilk<\/code> devient si \u00e9lev\u00e9e que le test d\u00e9tecte la moindre d\u00e9viation mineure par rapport \u00e0 une distribution normale th\u00e9orique id\u00e9ale. D\u00e8s lors, des micro-asym\u00e9tries sans impact pratique sur l&rsquo;inf\u00e9rence des estimateurs OLS g\u00e9n\u00e8rent syst\u00e9matiquement des valeurs de <span class=\"arabpsy-math-inline\" data-latex=\"p &lt; 0{,}001\">$p &lt; 0{,}001$<\/span>, conduisant \u00e0 rejeter l&#039;hypoth\u00e8se nulle \u00e0 tort d&#039;un point de vue fonctionnel. \u00c0 l&#039;oppos\u00e9, sur de tr\u00e8s petits \u00e9chantillons (<span class=\"arabpsy-math-inline\" data-latex=\"N &lt; 30\">$N &lt; 30$<\/span>), le test manque de puissance : de r\u00e9elles d\u00e9formations distributionnelles peuvent \u00e9chouer \u00e0 franchir le seuil de significativit\u00e9 statistique, exposant le praticien \u00e0 une erreur de seconde esp\u00e8ce (non-rejet d&#039;une hypoth\u00e8se nulle pourtant fausse).<\/p>\n<p>L&rsquo;application d&rsquo;un seuil critique universel fix\u00e9 \u00e0 5 % doit donc \u00eatre modul\u00e9e en fonction du volume d&rsquo;observations disponibles :<\/p>\n<ul>\n<li><strong>Petits \u00e9chantillons (<span class=\"arabpsy-math-inline\" data-latex=\"N &lt; 50\">$N &lt; 50$<\/span>) :<\/strong> Un seuil de tol\u00e9rance plus lib\u00e9ral (<span class=\"arabpsy-math-inline\" data-latex=\"\\alpha = 0{,}10\">$\\alpha = 0{,}10$<\/span>) est parfois envisag\u00e9 pour compenser le d\u00e9ficit structurel de puissance statistique.<\/li>\n<li><strong>Grands \u00e9chantillons (<span class=\"arabpsy-math-inline\" data-latex=\"N &gt; 500\">$N &gt; 500$<\/span>) :<\/strong> La valeur <span class=\"arabpsy-math-inline\" data-latex=\"p\">$p$<\/span> brute doit \u00eatre mise en perspective avec l&rsquo;amplitude du coefficient <span class=\"arabpsy-math-inline\" data-latex=\"W\">$W$<\/span> ; des valeurs de <span class=\"arabpsy-math-inline\" data-latex=\"W ge 0{,}98\">$W ge 0{,}98$<\/span> traduisent souvent une proximit\u00e9 empirique suffisante pour garantir la robustesse des mod\u00e8les param\u00e9triques, m\u00eame en pr\u00e9sence d&rsquo;une <span class=\"arabpsy-math-inline\" data-latex=\"p\">$p$<\/span>-value formellement significative.<\/li>\n<\/ul>\n<h2>5. Le test de Shapiro-Francia pour les \u00e9chantillons de grande taille<\/h2>\n<h3>5.1 Distinction algorithmique entre Shapiro-Wilk et Shapiro-Francia<\/h3>\n<p>Pour rem\u00e9dier aux contraintes computationnelles du calcul de la matrice inverse de variance-covariance des statistiques d&rsquo;ordre sur des \u00e9chantillons de taille mod\u00e9r\u00e9e \u00e0 grande, <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1080\/01621459.1972.10481232\">Shapiro et Francia ont propos\u00e9 en 1972<\/a> une variante algorithmique simplifi\u00e9e du test original. Alors que la statistique <span class=\"arabpsy-math-inline\" data-latex=\"W\">$W$<\/span> de Shapiro-Wilk incorpore l&rsquo;ensemble des covariances entre les statistiques d&rsquo;ordre th\u00e9oriques, le test de Shapiro-Francia fait l&rsquo;hypoth\u00e8se d&rsquo;une ind\u00e9pendance asymptotique relative de ces composantes ordonn\u00e9es, ce qui revient \u00e0 poser une matrice de covariance diagonale proportionnelle \u00e0 l&rsquo;identit\u00e9.<\/p>\n<p>La statistique r\u00e9sultante, formellement d\u00e9sign\u00e9e sous le symbole <span class=\"arabpsy-math-inline\" data-latex=\"W&#039;\">$W&rsquo;$<\/span>, repr\u00e9sente simplement le carr\u00e9 du coefficient de corr\u00e9lation lin\u00e9aire ordinaire calcul\u00e9 entre les statistiques d&rsquo;ordre observ\u00e9es dans l&rsquo;\u00e9chantillon et les esp\u00e9rances des statistiques d&rsquo;ordre correspondantes d&rsquo;une loi normale univari\u00e9e standard :<\/p>\n<p><code>W' = (sum(m_i * x_(i)))^2 \/ (sum(m_i^2) * sum((x_i - mean(x))^2))<\/code><\/p>\n<p>Cette simplification computationnelle, optimis\u00e9e par les travaux ult\u00e9rieurs de Royston en 1983, conf\u00e8re au test de Shapiro-Francia des propri\u00e9t\u00e9s asymptotiques remarquables. Il pr\u00e9sente une puissance statistique \u00e9quivalente \u2014 voire l\u00e9g\u00e8rement sup\u00e9rieure \u2014 \u00e0 celle du test <span class=\"arabpsy-math-inline\" data-latex=\"W\">$W$<\/span> original face aux distributions dot\u00e9es de queues lourdes ou platykurtiques, tout en supportant sans instabilit\u00e9 num\u00e9rique des \u00e9chantillons comprenant jusqu&rsquo;\u00e0 5 000 observations sous Stata.<\/p>\n<h3>5.2 Ex\u00e9cution et analyse sous Stata avec la commande sfrancia<\/h3>\n<p>L&rsquo;impl\u00e9mentation op\u00e9rationnelle sous Stata s&rsquo;effectue au moyen de la syntaxe <code>sfrancia<\/code>. La structure d&rsquo;appel demeure similaire \u00e0 celle de <code>swilk<\/code>, garantissant une coh\u00e9rence d&rsquo;utilisation pour l&rsquo;analyste :<\/p>\n<p><code>sfrancia price<\/code><\/p>\n<figure class=\"wpt-preserved-image\"><img loading=\"lazy\" decoding=\"async\" src=\"http:\/\/french.arabpsychology.com\/wp-content\/uploads\/2026\/09\/normalityTestStata5.png\" alt=\"Shapiro-Francia Test output in Stata\" \/><figcaption>Shapiro-Francia Test output in Stata<\/figcaption><\/figure>\n<p>Le tableau de r\u00e9sultats g\u00e9n\u00e9r\u00e9 par Stata fournit les indicateurs correspondants :<\/p>\n<ul>\n<li><code>W'<\/code> : Le coefficient d&rsquo;ajustement de Shapiro-Francia (<span class=\"arabpsy-math-inline\" data-latex=\"W&#039; = 0.751\">$W&rsquo; = 0.751$<\/span> pour la variable <code>price<\/code>).<\/li>\n<li><code>V'<\/code> : L&rsquo;indice d&rsquo;\u00e9chelle standardis\u00e9 analogue au facteur <span class=\"arabpsy-math-inline\" data-latex=\"V\">$V$<\/span> de Royston, calibr\u00e9 de sorte que sa valeur moyenne th\u00e9orique soit rigoureusement \u00e9gale \u00e0 1 sous l&rsquo;hypoth\u00e8se de normalit\u00e9.<\/li>\n<li><code>z<\/code> : La statistique de d\u00e9viation standardis\u00e9e utilis\u00e9e pour le positionnement sur la fonction de survie normale.<\/li>\n<li><code>Prob &gt; z<\/code> : La probabilit\u00e9 critique d&rsquo;erreur associ\u00e9e au rejet de l&rsquo;hypoth\u00e8se d&rsquo;ad\u00e9quation.<\/li>\n<\/ul>\n<p>Sur le plan des directives d&rsquo;application empirique, la r\u00e8gle de choix entre ces deux tests repose sur la taille de l&rsquo;\u00e9chantillon <span class=\"arabpsy-math-inline\" data-latex=\"N\">$N$<\/span> :<\/p>\n<ul>\n<li>Pour des cohortes comprises entre <span class=\"arabpsy-math-inline\" data-latex=\"N = 4\">$N = 4$<\/span> et <span class=\"arabpsy-math-inline\" data-latex=\"N = 200\">$N = 200$<\/span>, la commande <code>swilk<\/code> est g\u00e9n\u00e9ralement recommand\u00e9e en raison de l&rsquo;exactitude de ses coefficients de pond\u00e9ration sur les petites s\u00e9ries.<\/li>\n<li>Pour des \u00e9chantillons plus vastes, s&rsquo;\u00e9tendant de <span class=\"arabpsy-math-inline\" data-latex=\"N = 200\">$N = 200$<\/span> jusqu&rsquo;\u00e0 <span class=\"arabpsy-math-inline\" data-latex=\"N = 5,000\">$N = 5,000$<\/span>, la commande <code>sfrancia<\/code> est souvent privil\u00e9gi\u00e9e pour sa robustesse algorithmique et sa stabilit\u00e9 face aux distributions affichant des queues \u00e9paisses.<\/li>\n<\/ul>\n<h3>5.3 Applications empiriques par sous-groupes de variables<\/h3>\n<p>Il est fr\u00e9quent en recherche empirique que l&rsquo;hypoth\u00e8se de normalit\u00e9 doive \u00eatre examin\u00e9e distinctement selon les modalit\u00e9s d&rsquo;une variable de stratification ou de classification cat\u00e9gorielle. Par exemple, l&rsquo;analyste peut souhaiter \u00e9valuer la normalit\u00e9 de la consommation \u00e9nerg\u00e9tique (<code>mpg<\/code>) s\u00e9par\u00e9ment pour les v\u00e9hicules fabriqu\u00e9s sur le march\u00e9 int\u00e9rieur am\u00e9ricain et pour les v\u00e9hicules import\u00e9s, codifi\u00e9s par la variable binaire <code>foreign<\/code>.<\/p>\n<p>Stata permet d&rsquo;ex\u00e9cuter ces \u00e9valuations conditionnelles de fa\u00e7on \u00e9l\u00e9gante en associant les commandes de normalit\u00e9 au pr\u00e9fixe de partitionnement <code>bysort<\/code> :<\/p>\n<p><code>bysort foreign: sfrancia mpg<\/code><\/p>\n<p>L&rsquo;ex\u00e9cution de cette instruction partitionne les observations selon le statut d&rsquo;origine du v\u00e9hicule, puis applique de fa\u00e7on s\u00e9quentielle le test de Shapiro-Francia sur chaque sous-population. Cette approche permet de d\u00e9tecter si une non-normalit\u00e9 globale observ\u00e9e sur l&rsquo;\u00e9chantillon agr\u00e9g\u00e9 ne r\u00e9sulte pas simplement de la superposition de deux sous-populations homog\u00e8nes mais d\u00e9cal\u00e9es en moyenne (m\u00e9lange de distributions normales). Si le d\u00e9s\u00e9quilibre d&rsquo;effectif entre les groupes est important, l&rsquo;analyste devra interpr\u00e9ter les statistiques <span class=\"arabpsy-math-inline\" data-latex=\"W&#039;\">$W&rsquo;$<\/span> et <span class=\"arabpsy-math-inline\" data-latex=\"V&#039;\">$V&rsquo;$<\/span> en tenant compte de la perte relative de puissance dans le sous-groupe le moins repr\u00e9sent\u00e9.<\/p>\n<h2>6. L&rsquo;analyse formelle des coefficients d&rsquo;asym\u00e9trie et d&rsquo;aplatissement<\/h2>\n<h3>6.1 D\u00e9finition math\u00e9matique de l&rsquo;asym\u00e9trie et de l&rsquo;aplatissement<\/h3>\n<p>L&rsquo;architecture g\u00e9om\u00e9trique d&rsquo;une fonction de densit\u00e9 de probabilit\u00e9 peut \u00eatre d\u00e9compos\u00e9e au moyen de ses moments statistiques centr\u00e9s d&rsquo;ordre sup\u00e9rieur. Les moments du troisi\u00e8me et du quatri\u00e8me ordre d\u00e9crivent respectivement les d\u00e9formations d&rsquo;asym\u00e9trie (<a rel=\"noopener noreferrer\" href=\"https:\/\/mathworld.wolfram.com\/Skewness.html\">skewness<\/a>) et d&rsquo;aplatissement ou voussure (<a rel=\"noopener noreferrer\" href=\"https:\/\/mathworld.wolfram.com\/Kurtosis.html\">kurtosis<\/a>). Ces deux param\u00e8tres caract\u00e9risent la nature des \u00e9carts par rapport \u00e0 la distribution normale.<\/p>\n<p>Le coefficient th\u00e9orique d&rsquo;asym\u00e9trie de Fisher, not\u00e9 <span class=\"arabpsy-math-inline\" data-latex=\"\\gamma_1\">$\\gamma_1$<\/span> ou <span class=\"arabpsy-math-inline\" data-latex=\"\\sqrt{\\beta_1}\">$\\sqrt{\\beta_1}$<\/span>, correspond au troisi\u00e8me moment centr\u00e9 normalis\u00e9 par la troisi\u00e8me puissance de l&rsquo;\u00e9cart-type :<\/p>\n<p><code>gamma_1 = E[(X - mu)^3] \/ (sigma^2)^(3\/2)<\/code><\/p>\n<p>Une distribution normale th\u00e9orique est rigoureusement sym\u00e9trique par rapport \u00e0 sa moyenne, ce qui implique une asym\u00e9trie nulle (<span class=\"arabpsy-math-inline\" data-latex=\"\\gamma_1 = 0\">$\\gamma_1 = 0$<\/span>). Un coefficient positif (<span class=\"arabpsy-math-inline\" data-latex=\"\\gamma_1 &gt; 0\">$\\gamma_1 &gt; 0$<\/span>) signale une queue de distribution \u00e9tir\u00e9e vers les valeurs positives (asym\u00e9trie \u00e0 droite), tandis qu&rsquo;une valeur n\u00e9gative (<span class=\"arabpsy-math-inline\" data-latex=\"\\gamma_1 &lt; 0\">$\\gamma_1 &lt; 0$<\/span>) traduit un \u00e9tirement vers les valeurs n\u00e9gatives (asym\u00e9trie \u00e0 gauche).<\/p>\n<p>Le coefficient d&rsquo;aplatissement, not\u00e9 <span class=\"arabpsy-math-inline\" data-latex=\"\\beta_2\">$\\beta_2$<\/span>, d\u00e9coule quant \u00e0 lui du quatri\u00e8me moment centr\u00e9 standardis\u00e9 :<\/p>\n<p><code>beta_2 = E[(X - mu)^4] \/ (sigma^2)^2<\/code><\/p>\n<p>Pour une loi normale univari\u00e9e, ce ratio s&rsquo;\u00e9tablit exactement \u00e0 <span class=\"arabpsy-math-inline\" data-latex=\"\\beta_2 = 3\">$\\beta_2 = 3$<\/span>. En cons\u00e9quence, la litt\u00e9rature statistique emploie fr\u00e9quemment la notion de kurtosis excessif (ou aplatissement centr\u00e9), d\u00e9fini par la soustraction de cette constante :<\/p>\n<p><code>Kurtosis_excessif = beta_2 - 3<\/code><\/p>\n<p>Une distribution leptokurtique (<span class=\"arabpsy-math-inline\" data-latex=\"\\beta_2 &gt; 3\">$\\beta_2 &gt; 3$<\/span>) pr\u00e9sente un pic central plus aigu et des queues de distribution plus \u00e9paisses que la r\u00e9f\u00e9rence normale, ce qui accro\u00eet la probabilit\u00e9 de survenue d&rsquo;\u00e9v\u00e9nements extr\u00eames. Inversement, une distribution platykurtique (<span class=\"arabpsy-math-inline\" data-latex=\"\\beta_2 &lt; 3\">$\\beta_2 &lt; 3$<\/span>) affiche une morphologie plus aplatie et des queues plus courtes.<\/p>\n<h3>6.2 Extraction des moments empiriques avec la commande summarize, detail<\/h3>\n<p>Le moyen le plus direct d&rsquo;extraire les estimations d&rsquo;\u00e9chantillon des moments d&rsquo;asym\u00e9trie et d&rsquo;aplatissement sous Stata consiste \u00e0 utiliser l&rsquo;option d\u00e9taill\u00e9e de la proc\u00e9dure descriptive classique : <code>summarize, detail<\/code>.<\/p>\n<p>Consid\u00e9rons l&rsquo;application de cette instruction \u00e0 la variable <code>displacement<\/code> :<\/p>\n<p><code>summarize displacement, detail<\/code><\/p>\n<figure class=\"wpt-preserved-image\"><img loading=\"lazy\" decoding=\"async\" src=\"http:\/\/french.arabpsychology.com\/wp-content\/uploads\/2026\/09\/normalityTestStata6.png\" alt=\"Skewness and kurtosis for normality in Stata\" \/><figcaption>Skewness and kurtosis for normality in Stata<\/figcaption><\/figure>\n<p>Le tableau g\u00e9n\u00e9r\u00e9 par Stata structure les descripteurs statistiques en quatre colonnes principales :<\/p>\n<ul>\n<li>Les centiles caract\u00e9ristiques de la distribution (du 1er au 99e centile, incluant la m\u00e9diane au 50e rang).<\/li>\n<li>Les quatre plus petites et quatre plus grandes valeurs ordonn\u00e9es de la s\u00e9rie (utiles pour rep\u00e9rer d&rsquo;\u00e9ventuelles valeurs aberrantes).<\/li>\n<li>Les statistiques de tendance centrale et de dispersion \u00e9l\u00e9mentaires (le nombre d&rsquo;observations, la moyenne arithm\u00e9tique et l&rsquo;\u00e9cart-type).<\/li>\n<li>Les deux moments d&rsquo;ordre sup\u00e9rieur recherch\u00e9s : <strong>Skewness<\/strong> et <strong>Kurtosis<\/strong>.<\/li>\n<\/ul>\n<p>Pour la variable <code>displacement<\/code>, la sortie affiche un coefficient d&rsquo;asym\u00e9trie de <code>Skewness = 0.6726<\/code> et un coefficient d&rsquo;aplatissement de <code>Kurtosis = 2.2647<\/code>. L&rsquo;asym\u00e9trie est donc mod\u00e9r\u00e9ment positive, tandis que la valeur de kurtosis, inf\u00e9rieure au rep\u00e8re th\u00e9orique de 3, confirme une distribution platykurtique, d\u00e9pourvue de queues lourdes prononc\u00e9es.<\/p>\n<h3>6.3 Seuils critiques et r\u00e8gles de d\u00e9cision empiriques<\/h3>\n<p>L&rsquo;\u00e9valuation des moments empiriques permet d&rsquo;appr\u00e9cier la s\u00e9v\u00e9rit\u00e9 pratique des \u00e9carts \u00e0 la normalit\u00e9 sans d\u00e9pendre uniquement des tests d&rsquo;hypoth\u00e8se formels, souvent trop sensibles dans les grands \u00e9chantillons. La communaut\u00e9 statistique a \u00e9tabli des rep\u00e8res empiriques d&rsquo;acceptabilit\u00e9 pour l&rsquo;application des estimateurs lin\u00e9aires param\u00e9triques :<\/p>\n<ul>\n<li><strong>Zone de stricte normalit\u00e9 :<\/strong> L&rsquo;asym\u00e9trie empirique est comprise entre -0,5 et +0,5, et le kurtosis global se situe entre 2,5 et 3,5. Dans cet intervalle, les distorsions inf\u00e9rentielles sont consid\u00e9r\u00e9es comme n\u00e9gligeables.<\/li>\n<li><strong>Zone de d\u00e9viation mod\u00e9r\u00e9e :<\/strong> L&rsquo;asym\u00e9trie oscille entre l&rsquo;intervalle [-1, -0,5] ou [+0,5, +1], et le kurtosis se situe entre 2 et 4. Les m\u00e9thodes param\u00e9triques conservent une bonne robustesse sur des \u00e9chantillons de taille mod\u00e9r\u00e9e (<span class=\"arabpsy-math-inline\" data-latex=\"N &gt; 100\">$N &gt; 100$<\/span>).<\/li>\n<li><strong>Zone de d\u00e9viation s\u00e9v\u00e8re :<\/strong> L&rsquo;asym\u00e9trie absolue exc\u00e8de |1,5| ou |2|, et le kurtosis d\u00e9passe le seuil critique de 5 \u00e0 6 (ou tombe sous 1,5). Dans ce contexte, l&rsquo;inf\u00e9rence fond\u00e9e sur la loi de Student devient instable, imposant des corrections de type r\u00e9\u00e9chantillonnage ou transformations non lin\u00e9aires.<\/li>\n<\/ul>\n<p>Il est \u00e9galement possible de tester la nullit\u00e9 statistique des \u00e9carts par un ratio critique standardis\u00e9 simple :<\/p>\n<p><code>z_skewness = Skewness \/ sqrt(6 \/ N)<\/code><br \/>\n<code>z_kurtosis = (Kurtosis - 3) \/ sqrt(24 \/ N)<\/code><\/p>\n<p>Si la valeur absolue de ces statistiques d&rsquo;\u00e9cart standardis\u00e9 d\u00e9passe la valeur seuil de 1,96, l&rsquo;\u00e9cart par rapport \u00e0 l&rsquo;asym\u00e9trie ou \u00e0 l&rsquo;aplatissement normal est consid\u00e9r\u00e9 comme statistiquement significatif au seuil bilat\u00e9ral de 5 %.<\/p>\n<h2>7. Le test omnibus conjoint de D&rsquo;Agostino-Pearson dans Stata<\/h2>\n<h3>7.1 M\u00e9canisme inf\u00e9rentiel et formulation du test sktest<\/h3>\n<p>Les tests d&rsquo;hypoth\u00e8ses s\u00e9par\u00e9s portant sur l&rsquo;asym\u00e9trie ou l&rsquo;aplatissement pr\u00e9sentent l&rsquo;inconv\u00e9nient d&rsquo;ignorer la covariance distributionnelle conjointe de ces deux moments. Pour pallier cette limite, <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1093\/biomet\/60.3.613\">Ralph D&rsquo;Agostino et Egon Pearson ont formalis\u00e9 en 1973<\/a> un test omnibus d&rsquo;ad\u00e9quation. Ce test transforme simultan\u00e9ment les coefficients d&rsquo;asym\u00e9trie et de kurtosis empiriques en variables standardis\u00e9es approximativement gaussiennes, puis les combine au sein d&rsquo;une m\u00e9trique synth\u00e9tique globale.<\/p>\n<p>Dans sa version impl\u00e9ment\u00e9e sous Stata, la commande <code>sktest<\/code> int\u00e8gre les corrections analytiques d&rsquo;asym\u00e9trie d\u00e9velopp\u00e9es par D&rsquo;Agostino combin\u00e9es \u00e0 l&rsquo;algorithme d&rsquo;approximation polynomiale du kurtosis con\u00e7u par Patrick Royston en 1991. Soit <span class=\"arabpsy-math-inline\" data-latex=\"Z_1(s)\">$Z_1(s)$<\/span> le score gaussien issu de la transformation de l&rsquo;asym\u00e9trie empirique et <span class=\"arabpsy-math-inline\" data-latex=\"Z_2(k)\">$Z_2(k)$<\/span> le score correspondant d\u00e9riv\u00e9 de l&rsquo;aplatissement ; la statistique de test conjointe s&rsquo;exprime comme la somme quadratique de ces deux composantes ind\u00e9pendantes :<\/p>\n<p><code>K^2 = [Z_1(s)]^2 + [Z_2(k)]^2<\/code><\/p>\n<p>Sous l&rsquo;hypoth\u00e8se nulle <span class=\"arabpsy-math-inline\" data-latex=\"H_0\">$H_0$<\/span> stipulant que la population sous-jacente est r\u00e9gie par une loi normale, cette statistique quadratique <span class=\"arabpsy-math-inline\" data-latex=\"K^2\">$K^2$<\/span> converge asymptotiquement vers une loi du Chi-deux \u00e0 deux degr\u00e9s de libert\u00e9 (<span class=\"arabpsy-math-inline\" data-latex=\"\\chi^2(2)\">$\\chi^2(2)$<\/span>). Ce test omnibus fournit ainsi une \u00e9valuation globale : il suffit que la distribution s&rsquo;\u00e9carte significativement de la trajectoire normale soit par son asym\u00e9trie, soit par son aplatissement, soit par une combinaison des deux pour que la statistique conjointe conduise au rejet formel de l&rsquo;hypoth\u00e8se de normalit\u00e9.<\/p>\n<h3>7.2 Mise en \u0153uvre technique et interpr\u00e9tation des r\u00e9sultats de sktest<\/h3>\n<p>L&rsquo;ex\u00e9cution du test omnibus sous Stata repose sur une syntaxe directe : la commande <code>sktest<\/code> suivie de la liste des variables continues \u00e0 diagnostiquer. Appliquons cette analyse aux dimensions physiques des v\u00e9hicules :<\/p>\n<p><code>sktest weight length headroom<\/code><\/p>\n<p>Par d\u00e9faut, Stata applique la correction de Royston pour ajuster les moments sur les \u00e9chantillons de taille modeste. Le tableau r\u00e9capitulatif pr\u00e9sente les sorties suivantes pour chaque variable audit\u00e9e :<\/p>\n<ul>\n<li><code>Pr(Skewness)<\/code> : La p-valeur unilat\u00e9rale associ\u00e9e au test de conformit\u00e9 de l&rsquo;asym\u00e9trie seule.<\/li>\n<li><code>Pr(Kurtosis)<\/code> : La p-valeur unilat\u00e9rale testant l&rsquo;ad\u00e9quation de l&rsquo;aplatissement seul.<\/li>\n<li><code>adj chi2(2)<\/code> : La statistique globale combin\u00e9e du Chi-deux ajust\u00e9 \u00e0 deux degr\u00e9s de libert\u00e9.<\/li>\n<li><code>Prob &gt; chi2<\/code> : La valeur critique omnibus finale r\u00e9gissant la d\u00e9cision d&rsquo;inf\u00e9rence globale.<\/li>\n<\/ul>\n<p>Prenons l&rsquo;exemple de la variable <code>length<\/code> : si la p-valeur associ\u00e9e \u00e0 l&rsquo;asym\u00e9trie ressort \u00e0 <code>Pr(Skewness) = 0.3845<\/code> (indiquant une sym\u00e9trie satisfaisante), mais que le kurtosis affiche <code>Pr(Kurtosis) = 0.0210<\/code> (signalant un aplatissement non conforme), la statistique globale ajust\u00e9e <code>adj chi2(2)<\/code> affichera une probabilit\u00e9 globale de <code>Prob &gt; chi2 = 0.0450<\/code>. L&rsquo;hypoth\u00e8se nulle de normalit\u00e9 est alors rejet\u00e9e au seuil conventionnel de 5 %, le test identifiant avec pr\u00e9cision que l&rsquo;infraction distributionnelle provient exclusivement d&rsquo;une anomalie au niveau des queues de distribution, et non d&rsquo;un d\u00e9faut de sym\u00e9trie centrale.<\/p>\n<h3>7.3 Comparaison m\u00e9thodologique entre sktest et le test de Shapiro-Wilk<\/h3>\n<p>Le choix entre le test omnibus de D&rsquo;Agostino-Pearson (<code>sktest<\/code>) et le test d&rsquo;alignement de quantiles de Shapiro-Wilk (<code>swilk<\/code>) d\u00e9pend de la structure de l&rsquo;\u00e9chantillon et des objectifs du diagnostic m\u00e9thodologique. Ces deux tests reposent sur des prismes analytiques fondamentalement distincts :<\/p>\n<ul>\n<li><strong>Diagnostic des sources de d\u00e9viation :<\/strong> <code>sktest<\/code> pr\u00e9sente l&rsquo;avantage p\u00e9dagogique et pratique d&rsquo;isoler explicitement la part de non-normalit\u00e9 imputable \u00e0 l&rsquo;asym\u00e9trie de celle caus\u00e9e par l&rsquo;\u00e9paisseur des queues de distribution. \u00c0 l&rsquo;inverse, <code>swilk<\/code> condense l&rsquo;ensemble de l&rsquo;information dans un scalaire unique <span class=\"arabpsy-math-inline\" data-latex=\"W\">$W$<\/span>, sans expliciter la g\u00e9om\u00e9trie de la d\u00e9viation constat\u00e9e.<\/li>\n<li><strong>Puissance face aux distributions bimodales :<\/strong> Face \u00e0 des m\u00e9langes de distributions sym\u00e9triques mais nettement bimodales, <code>swilk<\/code> affiche g\u00e9n\u00e9ralement une puissance de d\u00e9tection sup\u00e9rieure \u00e0 celle de <code>sktest<\/code>, ce dernier pouvant \u00eatre tromp\u00e9 par une asym\u00e9trie apparente nulle compens\u00e9e par un kurtosis ambigu.<\/li>\n<li><strong>Comportement asymptotique :<\/strong> Sur des \u00e9chantillons interm\u00e9diaires (<span class=\"arabpsy-math-inline\" data-latex=\"100 le N le 1,000\">$100 le N le 1,000$<\/span>), <code>sktest<\/code> conserve d&rsquo;excellentes propri\u00e9t\u00e9s de calibrage et pr\u00e9sente une saturation moins rapide de sa p-valeur que <code>swilk<\/code>.<\/li>\n<\/ul>\n<p>Une bonne pratique m\u00e9thodologique consiste \u00e0 pratiquer une triangulation analytique : confronter les r\u00e9sultats de <code>sktest<\/code> \u00e0 ceux de <code>swilk<\/code>. Une convergence des conclusions vers le non-rejet conf\u00e8re un haut niveau de confiance quant \u00e0 la r\u00e9gularit\u00e9 distributionnelle des s\u00e9ries examin\u00e9es.<\/p>\n<h2>8. Les tests d&rsquo;ajustement globaux : Kolmogorov-Smirnov et Lilliefors<\/h2>\n<h3>8.1 Fondements du test de Kolmogorov-Smirnov \u00e0 un \u00e9chantillon<\/h3>\n<p>Le test d&rsquo;ajustement global de <a rel=\"noopener noreferrer\" href=\"https:\/\/mathworld.wolfram.com\/Kolmogorov-SmirnovTest.html\">Kolmogorov-Smirnov<\/a> repose sur une m\u00e9trique non param\u00e9trique \u00e9valuant la distance maximale observ\u00e9e entre la fonction de r\u00e9partition empirique d&rsquo;un \u00e9chantillon et la fonction de r\u00e9partition cumulative th\u00e9orique assign\u00e9e a priori sous l&rsquo;hypoth\u00e8se nulle.<\/p>\n<p>Soit <span class=\"arabpsy-math-inline\" data-latex=\"S_n(x)\">$S_n(x)$<\/span> la fonction d&rsquo;escalier repr\u00e9sentant la proportion cumul\u00e9e des observations empiriques inf\u00e9rieures ou \u00e9gales au r\u00e9el <span class=\"arabpsy-math-inline\" data-latex=\"x\">$x$<\/span>, et soit <span class=\"arabpsy-math-inline\" data-latex=\"F_0(x)\">$F_0(x)$<\/span> la fonction de r\u00e9partition th\u00e9orique normale continue postul\u00e9e. La statistique de Kolmogorov-Smirnov correspond \u00e0 la borne sup\u00e9rieure (distance supr\u00eame) des \u00e9carts absolus verticaux :<\/p>\n<p><code>D = sup_x |S_n(x) - F_0(x)|<\/code><\/p>\n<p>Dans l&rsquo;environnement Stata, la commande native correspondante est <code>ksmirnov<\/code>. Toutefois, une condition th\u00e9orique fondamentale restreint son application directe : le test de Kolmogorov-Smirnov classique impose que la distribution th\u00e9orique <span class=\"arabpsy-math-inline\" data-latex=\"F_0(x)\">$F_0(x)$<\/span> soit <em>enti\u00e8rement sp\u00e9cifi\u00e9e a priori<\/em>, c&rsquo;est-\u00e0-dire que ses param\u00e8tres de moyenne (<span class=\"arabpsy-math-inline\" data-latex=\"\\mu_0\">$\\mu_0$<\/span>) et de variance (<span class=\"arabpsy-math-inline\" data-latex=\"\\sigma^2_0\">$\\sigma^2_0$<\/span>) soient connus ind\u00e9pendamment des donn\u00e9es de l&rsquo;\u00e9chantillon. Pour tester une variable continue avec cette proc\u00e9dure, il est donc n\u00e9cessaire de calculer au pr\u00e9alable les valeurs standardis\u00e9es z-score sur la base d&rsquo;hypoth\u00e8ses fix\u00e9es :<\/p>\n<p><code>summarize weight<\/code><br \/>\n<code>generate double z_weight = (weight - r(mean)) \/ r(sd)<\/code><br \/>\n<code>ksmirnov z_weight = normal(z_weight)<\/code><\/p>\n<p>La statistique <span class=\"arabpsy-math-inline\" data-latex=\"D\">$D$<\/span> retourn\u00e9e correspond au plus grand \u00e9cart absolu constat\u00e9 sur l&rsquo;ensemble de la fonction de r\u00e9partition empirique par rapport \u00e0 la loi normale centr\u00e9e r\u00e9duite th\u00e9orique standardis\u00e9e.<\/p>\n<h3>8.2 Le probl\u00e8me des param\u00e8tres estim\u00e9s et l&rsquo;alternative de Lilliefors<\/h3>\n<p>L&rsquo;utilisation de la commande standard <code>ksmirnov<\/code> pour tester la normalit\u00e9 d&rsquo;une s\u00e9rie pose un probl\u00e8me th\u00e9orique majeur lorsque les param\u00e8tres de localisation (<span class=\"arabpsy-math-inline\" data-latex=\"\\mu\">$\\mu$<\/span>) et d&rsquo;\u00e9chelle (<span class=\"arabpsy-math-inline\" data-latex=\"\\sigma\">$\\sigma$<\/span>) ne sont pas connus d&rsquo;avance, mais estim\u00e9s \u00e0 partir de l&rsquo;\u00e9chantillon lui-m\u00eame via la moyenne <span class=\"arabpsy-math-inline\" data-latex=\"\\bar{x}\">$\\bar{x}$<\/span> et l&rsquo;\u00e9cart-type empirique <span class=\"arabpsy-math-inline\" data-latex=\"s\">$s$<\/span>. Cette pratique tr\u00e8s courante invalide les tables de probabilit\u00e9 critiques de Kolmogorov-Smirnov standard.<\/p>\n<p>En effet, le fait d&rsquo;ajuster la fonction th\u00e9orique normale sur les donn\u00e9es observ\u00e9es rapproche m\u00e9caniquement la courbe th\u00e9orique de la distribution empirique. La distance supr\u00eame calcul\u00e9e <span class=\"arabpsy-math-inline\" data-latex=\"D\">$D$<\/span> est alors syst\u00e9matiquement sous-estim\u00e9e par rapport \u00e0 ce qu&rsquo;elle serait face \u00e0 une fonction ind\u00e9pendante. En utilisant les valeurs critiques usuelles de la commande <code>ksmirnov<\/code>, le test devient excessivement conservateur : la probabilit\u00e9 r\u00e9elle de commettre une erreur de premi\u00e8re esp\u00e8ce s&rsquo;effondre bien en de\u00e7\u00e0 du seuil nominal de 5 %, ce qui conduit \u00e0 accepter \u00e0 tort l&rsquo;hypoth\u00e8se de normalit\u00e9 dans des situations de d\u00e9viation marqu\u00e9e (d\u00e9ficit prononc\u00e9 de puissance statistique).<\/p>\n<p>Pour corriger ce biais d&rsquo;estimation, Hubert Lilliefors a d\u00e9riv\u00e9 en 1967 une table de valeurs critiques calibr\u00e9es par simulations de Monte Carlo, adapt\u00e9e \u00e0 la situation o\u00f9 les param\u00e8tres sont estim\u00e9s sur l&rsquo;\u00e9chantillon. Bien que Stata ne propose pas de commande native interne d\u00e9di\u00e9e au test de Lilliefors, la communaut\u00e9 m\u00e9thodologique a d\u00e9velopp\u00e9 des modules externes ex\u00e9cutables, accessibles via le d\u00e9p\u00f4t SSC :<\/p>\n<p><code>ssc install lilliefors<\/code><br \/>\n<code>lilliefors weight<\/code><\/p>\n<p>L&rsquo;algorithme de Lilliefors recalibre correctement la valeur critique du test en fonction de la taille <span class=\"arabpsy-math-inline\" data-latex=\"N\">$N$<\/span>, r\u00e9tablissant une taille de test conforme au seuil <span class=\"arabpsy-math-inline\" data-latex=\"\\alpha\">$\\alpha$<\/span> nominal fix\u00e9.<\/p>\n<h3>8.3 \u00c9valuation comparative de l&rsquo;efficacit\u00e9 face \u00e0 swilk et sktest<\/h3>\n<p>Sur le plan de l&rsquo;efficacit\u00e9 statistique compar\u00e9e, la litt\u00e9rature m\u00e9thodologique sp\u00e9cialis\u00e9e d\u00e9conseille de mani\u00e8re unanime l&#8217;emploi du test de Kolmogorov-Smirnov (m\u00eame corrig\u00e9 par Lilliefors) pour diagnostiquer la normalit\u00e9 d&rsquo;une variable continue. Les \u00e9tudes comparatives de puissance par simulations de Monte Carlo mettent syst\u00e9matiquement en \u00e9vidence l&rsquo;inf\u00e9riorit\u00e9 de la statistique de Kolmogorov-Smirnov face aux tests de Shapiro-Wilk (<code>swilk<\/code>), de Shapiro-Francia (<code>sfrancia<\/code>) et de D&rsquo;Agostino-Pearson (<code>sktest<\/code>).<\/p>\n<p>Ce d\u00e9ficit structurel de puissance s&rsquo;explique par la nature g\u00e9om\u00e9trique de la statistique <span class=\"arabpsy-math-inline\" data-latex=\"D\">$D$<\/span> :<\/p>\n<ul>\n<li>Le test de Kolmogorov-Smirnov ne prend en compte que la distance verticale maximale unique observ\u00e9e sur l&rsquo;ensemble du domaine de distribution, ignorant les \u00e9carts cumulatifs r\u00e9guliers r\u00e9partis sur l&rsquo;ensemble du support.<\/li>\n<li>Les d\u00e9formations asym\u00e9triques ou les queues de distribution \u00e9paisses g\u00e9n\u00e8rent souvent des \u00e9carts verticaux faibles en valeur absolue sur la fonction de r\u00e9partition cumulative, mais fondamentaux quant \u00e0 la forme de la densit\u00e9. Ces \u00e9carts passent inaper\u00e7us sous la statistique <span class=\"arabpsy-math-inline\" data-latex=\"D\">$D$<\/span>, alors qu&rsquo;ils sont imm\u00e9diatement d\u00e9tect\u00e9s par la corr\u00e9lation des quantiles de <code>swilk<\/code> ou par les moments d&rsquo;ordre sup\u00e9rieur de <code>sktest<\/code>.<\/li>\n<\/ul>\n<p>Le recours \u00e0 la commande <code>ksmirnov<\/code> pour \u00e9valuer la normalit\u00e9 rel\u00e8ve donc d&rsquo;un h\u00e9ritage historique plut\u00f4t que d&rsquo;un choix m\u00e9thodologique moderne. Son usage dans les protocoles contemporains sous Stata doit \u00eatre r\u00e9serv\u00e9 \u00e0 l&rsquo;\u00e9valuation d&rsquo;autres lois th\u00e9oriques sp\u00e9cifiques enti\u00e8rement param\u00e9tr\u00e9es a priori, au profit exclusif de <code>swilk<\/code> et <code>sktest<\/code> pour la distribution gaussienne.<\/p>\n<h2>9. Diagnostics de la normalit\u00e9 des r\u00e9sidus dans les mod\u00e8les de r\u00e9gression lin\u00e9aire<\/h2>\n<h3>9.1 Le postulat de Gauss-Markov et la distribution des erreurs<\/h3>\n<p>L&rsquo;estimation des coefficients d&rsquo;un mod\u00e8le de r\u00e9gression lin\u00e9aire par les moindres carr\u00e9s ordinaires s&rsquo;appuie sur le th\u00e9or\u00e8me classique de <a rel=\"noopener noreferrer\" href=\"https:\/\/mathworld.wolfram.com\/Gauss-MarkovTheorem.html\">Gauss-Markov<\/a>. Selon ce th\u00e9or\u00e8me, si les perturbations stochastiques inobservables ont une esp\u00e9rance conditionnelle nulle (<span class=\"arabpsy-math-inline\" data-latex=\"E[epsilon|X] = 0\">$E[epsilon|X] = 0$<\/span>) et sont sph\u00e9riques (homosc\u00e9dastiques et non autocorr\u00e9l\u00e9es), les estimateurs OLS constituent les meilleurs estimateurs lin\u00e9aires sans biais (BLUE, <em>Best Linear Unbiased Estimators<\/em>). Il est fondamental de souligner que cette propri\u00e9t\u00e9 d&rsquo;efficacit\u00e9 lin\u00e9aire minimale <strong>ne requiert en aucun cas<\/strong> que les perturbations suivent une distribution normale.<\/p>\n<p>L&rsquo;hypoth\u00e8se de normalit\u00e9 des termes d&rsquo;erreur :<\/p>\n<p><code>epsilon | X ~ N(0, sigma^2 * I)<\/code><\/p>\n<p>n&rsquo;intervient qu&rsquo;au stade de l&rsquo;inf\u00e9rence statistique sur les param\u00e8tres estim\u00e9s. Elle est indispensable pour \u00e9tablir la distribution exacte des coefficients sous des \u00e9chantillons de dimension finie. Ce n&rsquo;est que sous cette hypoth\u00e8se que les ratios des coefficients estim\u00e9s divis\u00e9s par leurs erreurs-types respectives suivent exactement une distribution <span class=\"arabpsy-math-inline\" data-latex=\"t\">$t$<\/span> de Student \u00e0 <span class=\"arabpsy-math-inline\" data-latex=\"(N - k)\">$(N &#8211; k)$<\/span> degr\u00e9s de libert\u00e9, et que les tests d&rsquo;hypoth\u00e8ses lin\u00e9aires conjointes suivent une loi <span class=\"arabpsy-math-inline\" data-latex=\"F\">$F$<\/span> de Fisher-Snedecor.<\/p>\n<p>Une distinction m\u00e9thodologique essentielle s\u00e9pare la distribution marginale inconditionnelle de la variable expliqu\u00e9e (<span class=\"arabpsy-math-inline\" data-latex=\"Y\">$Y$<\/span>) de sa distribution conditionnelle sachant les r\u00e9gresseurs (<span class=\"arabpsy-math-inline\" data-latex=\"Y|X\">$Y|X$<\/span>). Exiger que la variable d\u00e9pendante brute soit normalement distribu\u00e9e constitue une erreur conceptuelle : une variable d\u00e9pendante peut pr\u00e9senter une distribution bimodale ou fortement \u00e9tir\u00e9e simplement parce qu&rsquo;elle d\u00e9pend d&rsquo;un pr\u00e9dicteur binaire ou d&rsquo;une covariable asym\u00e9trique. C&rsquo;est uniquement la composante r\u00e9siduelle inexpliqu\u00e9e par le mod\u00e8le qui doit satisfaire au postulat de normalit\u00e9 gaussienne.<\/p>\n<h3>9.2 Extraction, standardisation et stockage des r\u00e9sidus post-estimation<\/h3>\n<p>Dans la pratique analytique sous Stata, l&rsquo;audit de normalit\u00e9 intervient toujours comme une \u00e9tape diagnostique post-estimation, ex\u00e9cut\u00e9e cons\u00e9cutivement \u00e0 l&rsquo;ajustement du mod\u00e8le de r\u00e9gression via l&rsquo;instruction <code>regress<\/code>. Consid\u00e9rons l&rsquo;estimation du prix d&rsquo;un v\u00e9hicule en fonction de son kilom\u00e9trage, de son poids et de sa provenance :<\/p>\n<p><code>regress price mpg weight i.foreign<\/code><\/p>\n<p>Une fois le vecteur des coefficients obtenu, l&rsquo;analyste extrait les perturbations r\u00e9siduelles \u00e0 l&rsquo;aide de la commande <code>predict<\/code>. Stata propose plusieurs types de r\u00e9sidus adapt\u00e9s \u00e0 des finalit\u00e9s diagnostiques compl\u00e9mentaires :<\/p>\n<ul>\n<li><code>predict res_bruts, residuals<\/code> : Calcule les r\u00e9sidus ordinaires d\u00e9finis par la simple diff\u00e9rence alg\u00e9brique entre la valeur observ\u00e9e et la valeur ajust\u00e9e : <span class=\"arabpsy-math-inline\" data-latex=\"e_i = y_i - \\hat{y}_i\">$e_i = y_i &#8211; \\hat{y}_i$<\/span>.<\/li>\n<li><code>predict res_std, rstandard<\/code> : Calcule les r\u00e9sidus standardis\u00e9s (internement studentis\u00e9s), obtenus en divisant les r\u00e9sidus bruts par leur \u00e9cart-type d&rsquo;\u00e9chantillon corrig\u00e9 de l&rsquo;effet de levier individuel :\n<p><code>r_i = e_i \/ (s * sqrt(1 - h_ii))<\/code><\/p>\n<p>o\u00f9 <span class=\"arabpsy-math-inline\" data-latex=\"h_{ii}\">$h_{ii}$<\/span> correspond au <span class=\"arabpsy-math-inline\" data-latex=\"i\">$i$<\/span>-\u00e8me terme diagonal de la matrice de projection chapeau (hat matrix).<\/li>\n<li><code>predict res_stud, rstudent<\/code> : Calcule les r\u00e9sidus studentis\u00e9s (externement studentis\u00e9s ou jackknife), o\u00f9 la variance de l&rsquo;erreur est estim\u00e9e en omettant s\u00e9lectivement la <span class=\"arabpsy-math-inline\" data-latex=\"i\">$i$<\/span>-\u00e8me observation. Ce calcul \u00e9vite qu&rsquo;une valeur aberrante n&rsquo;att\u00e9nue elle-m\u00eame l&rsquo;ampleur de son r\u00e9sidu standardis\u00e9.<\/li>\n<\/ul>\n<p>Pour auditer la normalit\u00e9 de la composante d&rsquo;erreur r\u00e9siduelle, l&rsquo;utilisation des r\u00e9sidus standardis\u00e9s ou studentis\u00e9s est formellement recommand\u00e9e : elle neutralise les variations artificielles de variance induites par l&rsquo;\u00e9loignement spatial des observations au sein de l&rsquo;espace des covariables (h\u00e9t\u00e9rog\u00e9n\u00e9it\u00e9 des leviers <span class=\"arabpsy-math-inline\" data-latex=\"h_{ii}\">$h_{ii}$<\/span>).<\/p>\n<h3>9.3 Batterie compl\u00e8te d&rsquo;\u00e9valuation des r\u00e9sidus g\u00e9n\u00e9r\u00e9s<\/h3>\n<p>Une fois la variable r\u00e9siduelle standardis\u00e9e stock\u00e9e en m\u00e9moire, l&rsquo;analyste d\u00e9ploie l&rsquo;ensemble des outils graphiques et inf\u00e9rentiels pour valider la conformit\u00e9 distributionnelle de l&rsquo;erreur du mod\u00e8le. Cette d\u00e9marche diagnostique s&rsquo;articule autour de trois \u00e9tapes m\u00e9thodologiques compl\u00e9mentaires :<\/p>\n<p>1. <strong>Inspection de densit\u00e9 non param\u00e9trique :<\/strong> L&rsquo;utilisation combin\u00e9e d&rsquo;un estimateur \u00e0 noyau de densit\u00e9 de type Epanechnikov via <code>kdensity<\/code>, avec superposition de la r\u00e9f\u00e9rence gaussienne th\u00e9orique :<\/p>\n<p><code>kdensity res_stud, normal kernel(epanechnikov) title(\"Densit\u00e9 \u00e0 noyau des r\u00e9sidus studentis\u00e9s\")<\/code><\/p>\n<p>Ce graphique permet de rep\u00e9rer d&rsquo;\u00e9ventuelles asym\u00e9tries locales ou aplatissements excessifs que le lissage d&rsquo;un histogramme standard pourrait masquer.<\/p>\n<p>2. <strong>Tests inf\u00e9rentiels formels sur les r\u00e9sidus :<\/strong> Les tests de Shapiro-Wilk et de D&rsquo;Agostino-Pearson sont appliqu\u00e9s directement sur la s\u00e9rie des r\u00e9sidus standardis\u00e9s g\u00e9n\u00e9r\u00e9s :<\/p>\n<p><code>swilk res_stud<\/code><br \/>\n<code>sktest res_stud<\/code><\/p>\n<p>Si la statistique conjointe de <code>sktest<\/code> et le coefficient <span class=\"arabpsy-math-inline\" data-latex=\"W\">$W$<\/span> de <code>swilk<\/code> indiquent un rejet formel de l&rsquo;hypoth\u00e8se nulle (<span class=\"arabpsy-math-inline\" data-latex=\"p &lt; 0{,}05\">$p &lt; 0{,}05$<\/span>), la validit\u00e9 des intervalles de confiance param\u00e9triques des coefficients de r\u00e9gression se trouve compromise sur petit \u00e9chantillon.<\/p>\n<p>3. <strong>V\u00e9rification conjointe de l&rsquo;homosc\u00e9dasticit\u00e9 :<\/strong> L&rsquo;analyse de la distribution r\u00e9siduelle gagne \u00e0 \u00eatre confront\u00e9e au graphique de dispersion reliant les r\u00e9sidus aux valeurs pr\u00e9dites :<\/p>\n<p><code>rvfplot, yline(0) title(\"Dispersion des r\u00e9sidus vs valeurs ajust\u00e9es\")<\/code><\/p>\n<p>La pr\u00e9sence d&rsquo;h\u00e9t\u00e9rosc\u00e9dasticit\u00e9 (structure en entonnoir sur le graphique <code>rvfplot<\/code>) biaise fr\u00e9quemment les indicateurs de normalit\u00e9 globale des r\u00e9sidus. Dans ce cas, l&rsquo;anomalie distributionnelle apparente n&rsquo;est souvent que la cons\u00e9quence directe d&rsquo;une instabilit\u00e9 sous-jacente de la variance des erreurs conditionnelles.<\/p>\n<h2>10. \u00c9valuation de la normalit\u00e9 multivari\u00e9e sous Stata<\/h2>\n<h3>10.1 Concepts fondamentaux de la loi normale multivari\u00e9e<\/h3>\n<p>Dans les architectures statistiques fond\u00e9es sur l&rsquo;analyse simultan\u00e9e de plusieurs variables interd\u00e9pendantes \u2014 telles que les r\u00e9gressions multivari\u00e9es, les mod\u00e8les d&rsquo;\u00e9quations structurelles, l&rsquo;analyse multivari\u00e9e de la variance (<a rel=\"noopener noreferrer\" href=\"https:\/\/en.wikipedia.org\/wiki\/Multivariate_analysis_of_variance\">MANOVA<\/a>) ou l&rsquo;analyse en composantes principales \u2014, l&rsquo;hypoth\u00e8se de travail repose sur le respect d&rsquo;une loi normale multivari\u00e9e par le vecteur al\u00e9atoire <span class=\"arabpsy-math-inline\" data-latex=\"X = (X_1, X_2, dots, X_p)&#039;\">$X = (X_1, X_2, dots, X_p)&rsquo;$<\/span>. Une distribution normale multivari\u00e9e d&rsquo;ordre <span class=\"arabpsy-math-inline\" data-latex=\"p\">$p$<\/span>, not\u00e9e <span class=\"arabpsy-math-inline\" data-latex=\"N_p(\\mu, \\Sigma)\">$N_p(\\mu, \\Sigma)$<\/span>, est enti\u00e8rement caract\u00e9ris\u00e9e par son vecteur d&rsquo;esp\u00e9rances math\u00e9matiques <span class=\"arabpsy-math-inline\" data-latex=\"\\mu\">$\\mu$<\/span> de dimension <span class=\"arabpsy-math-inline\" data-latex=\"(p \\times 1)\">$(p \\times 1)$<\/span> et sa matrice sym\u00e9trique semi-d\u00e9finie positive de variance-covariance <span class=\"arabpsy-math-inline\" data-latex=\"\\Sigma\">$\\Sigma$<\/span> de dimension <span class=\"arabpsy-math-inline\" data-latex=\"(p \\times p)\">$(p \\times p)$<\/span>.<\/p>\n<p>Une propri\u00e9t\u00e9 g\u00e9om\u00e9trique fondamentale de cette distribution r\u00e9side dans le fait que ses contours d&rsquo;\u00e9quidensit\u00e9 \u00e9pousent des formes ellipso\u00efdales dans l&rsquo;espace euclidien \u00e0 <span class=\"arabpsy-math-inline\" data-latex=\"p\">$p$<\/span> dimensions. Pour chaque individu statistique <span class=\"arabpsy-math-inline\" data-latex=\"i\">$i$<\/span>, l&rsquo;\u00e9loignement spatial pond\u00e9r\u00e9 par rapport au centre de gravit\u00e9 du nuage de points est mesur\u00e9 par le carr\u00e9 de la distance statistique de Mahalanobis :<\/p>\n<p><code>D_i^2 = (x_i - mean(x))' * Sigma^(-1) * (x_i - mean(x))<\/code><\/p>\n<p>Sous l&rsquo;hypoth\u00e8se nulle d&rsquo;une stricte conformit\u00e9 \u00e0 la loi normale multivari\u00e9e, ces distances de Mahalanobis quadratiques <span class=\"arabpsy-math-inline\" data-latex=\"D_i^2\">$D_i^2$<\/span> suivent th\u00e9oriquement une distribution du Chi-deux \u00e0 <span class=\"arabpsy-math-inline\" data-latex=\"p\">$p$<\/span> degr\u00e9s de libert\u00e9 :<\/p>\n<p><code>D_i^2 ~ chi^2(p)<\/code><\/p>\n<p>La simple juxtaposition de normalit\u00e9s univari\u00e9es marginales ne suffit pas \u00e0 garantir la normalit\u00e9 conjointe du vecteur al\u00e9atoire. Bien que la normalit\u00e9 multivari\u00e9e implique obligatoirement que chaque composante prise individuellement suive une loi normale univari\u00e9e, la r\u00e9ciproque math\u00e9matique est fausse : il est tout \u00e0 fait possible de concevoir un syst\u00e8me o\u00f9 toutes les marges univari\u00e9es sont gaussiennes sans que leur distribution conjointe ne respecte la structure matricielle normale.<\/p>\n<h3>10.2 La commande mvtest normality et les tests de Mardia<\/h3>\n<p>Stata dispose d&rsquo;un module d&rsquo;analyse d\u00e9di\u00e9 \u00e0 l&rsquo;\u00e9valuation de la normalit\u00e9 multivari\u00e9e, accessible via l&rsquo;instruction native <code>mvtest normality<\/code>. Ce module int\u00e8gre les statistiques multivari\u00e9es d\u00e9velopp\u00e9es par <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1093\/biomet\/57.3.519\">Kanti Mardia en 1970<\/a>, qui d\u00e9composent l&rsquo;\u00e9valuation conjointe selon deux indices sp\u00e9cifiques : l&rsquo;asym\u00e9trie multivari\u00e9e de Mardia (<span class=\"arabpsy-math-inline\" data-latex=\"b_{1,p}\">$b_{1,p}$<\/span>) et l&rsquo;aplatissement multivari\u00e9 de Mardia (<span class=\"arabpsy-math-inline\" data-latex=\"b_{2,p}\">$b_{2,p}$<\/span>).<\/p>\n<p>Pour soumettre simultan\u00e9ment les variables <code>price<\/code>, <code>weight<\/code>, <code>mpg<\/code> et <code>length<\/code> \u00e0 cette proc\u00e9dure, la syntaxe s&rsquo;articule comme suit :<\/p>\n<p><code>mvtest normality price weight mpg length<\/code><\/p>\n<p>Par d\u00e9faut, l&rsquo;instruction ex\u00e9cute le test de Mardia et retourne une matrice de sortie structur\u00e9e autour de trois diagnostics compl\u00e9mentaires :<\/p>\n<ul>\n<li><strong>Asym\u00e9trie de Mardia (Mardia Skewness) :<\/strong> Ce test \u00e9value la somme des cubes des distances statistiques crois\u00e9es entre observations. Stata affiche la statistique calcul\u00e9e <span class=\"arabpsy-math-inline\" data-latex=\"b_{1,p}\">$b_{1,p}$<\/span>, sa conversion standardis\u00e9e sous la forme d&rsquo;un score <span class=\"arabpsy-math-inline\" data-latex=\"\\chi^2\">$\\chi^2$<\/span>, ainsi que la p-valeur associ\u00e9e. Une asym\u00e9trie multivari\u00e9e significative indique une d\u00e9formation d&rsquo;orientation dans l&rsquo;espace multidimensionnel.<\/li>\n<li><strong>Aplatissement de Mardia (Mardia Kurtosis) :<\/strong> Ce test examine la dispersion du carr\u00e9 des distances de Mahalanobis autour de leur esp\u00e9rance th\u00e9orique <span class=\"arabpsy-math-inline\" data-latex=\"p(p + 2)\">$p(p + 2)$<\/span>. Le logiciel produit la statistique d&rsquo;aplatissement <span class=\"arabpsy-math-inline\" data-latex=\"b_{2,p}\">$b_{2,p}$<\/span>, un score standardis\u00e9 <span class=\"arabpsy-math-inline\" data-latex=\"z\">$z$<\/span>, et la probabilit\u00e9 critique bilat\u00e9rale correspondante. Une d\u00e9viation significative \u00e0 ce niveau r\u00e9v\u00e8le la pr\u00e9sence de queues lourdes multidimensionnelles ou de valeurs aberrantes multivari\u00e9es.<\/li>\n<li><strong>Test conjointe d&rsquo;asym\u00e9trie et de kurtosis :<\/strong> Combine les deux composantes pr\u00e9c\u00e9dentes en un crit\u00e8re omnibus unique.<\/li>\n<\/ul>\n<p>Le rejet de l&rsquo;hypoth\u00e8se nulle sous la commande <code>mvtest normality<\/code> avertit l&rsquo;analyste que les estimateurs du maximum de vraisemblance standard risquent de produire des erreurs-types biais\u00e9es lors de l&rsquo;estimation de mod\u00e8les d&rsquo;\u00e9quations structurelles multivari\u00e9s.<\/p>\n<h3>10.3 Tests de Doornik-Hansen et de Henze-Zirkler<\/h3>\n<p>Outre la m\u00e9thodologie de Mardia, Stata propose des proc\u00e9dures alternatives r\u00e9put\u00e9es pour leurs propri\u00e9t\u00e9s de puissance statistique sup\u00e9rieures face \u00e0 des corr\u00e9lations inter-variables \u00e9lev\u00e9es, accessibles via les options de la commande <code>mvtest<\/code> :<\/p>\n<p><code>mvtest normality price weight mpg length, doornikhansen<\/code><br \/>\n<code>mvtest normality price weight mpg length, henzezirkler<\/code><\/p>\n<p>Le test de Doornik-Hansen (1994) s&rsquo;appuie sur une transformation des coefficients d&rsquo;asym\u00e9trie et de kurtosis univari\u00e9s apr\u00e8s orthogonalisation compl\u00e8te des variables par la d\u00e9composition de Cholesky de leur matrice de corr\u00e9lation. La statistique globale r\u00e9sultante converge vers une loi du Chi-deux \u00e0 <span class=\"arabpsy-math-inline\" data-latex=\"2p\">$2p$<\/span> degr\u00e9s de libert\u00e9. Cette approche affiche de meilleures propri\u00e9t\u00e9s d&rsquo;\u00e9chantillonnage sur les s\u00e9ries r\u00e9duites (<span class=\"arabpsy-math-inline\" data-latex=\"N &lt; 100\">$N &lt; 100$<\/span>) que les statistiques asymptotiques de Mardia.<\/p>\n<p>Le test de Henze-Zirkler (1990) mobilise quant \u00e0 lui une fonction de distance globale bas\u00e9e sur la transform\u00e9e de Fourier de la fonction caract\u00e9ristique multivari\u00e9e empirique. Une statistique <span class=\"arabpsy-math-inline\" data-latex=\"HZ\">$HZ$<\/span> \u00e9lev\u00e9e indique un \u00e9cart significatif par rapport \u00e0 l&rsquo;hypoth\u00e8se de normalit\u00e9 conjointe. En mod\u00e9lisation causale (SEM), la violation av\u00e9r\u00e9e de la normalit\u00e9 multivari\u00e9e sous ces tests impose d&rsquo;abandonner l&rsquo;estimateur classique du maximum de vraisemblance au profit d&rsquo;estimateurs quasi-maximum de vraisemblance robustes (QML) ou d&rsquo;algorithmes asymptotiquement ind\u00e9pendants de la distribution (ADF \/ Satorra-Bentler).<\/p>\n<h2>11. Strat\u00e9gies m\u00e9thodologiques face \u00e0 la violation de l&rsquo;hypoth\u00e8se de normalit\u00e9<\/h2>\n<h3>11.1 Transformations math\u00e9matiques univari\u00e9es classiques<\/h3>\n<p>Lorsqu&rsquo;une variable continue s&rsquo;\u00e9carte substantiellement de la loi normale par suite d&rsquo;une forte asym\u00e9trie unimodale, la solution m\u00e9thodologique la plus directe consiste \u00e0 appliquer une transformation math\u00e9matique non lin\u00e9aire monotone. L&rsquo;objectif est de recomprimer les \u00e9carts relatifs entre les observations de mani\u00e8re \u00e0 rapprocher la distribution d&rsquo;une forme sym\u00e9trique gaussienne.<\/p>\n<p>Pour les variables affichant une asym\u00e9trie positive \u00e9tir\u00e9e vers la droite (telles que le revenu, le co\u00fbt ou les concentrations biologiques), la transformation logarithmique n\u00e9p\u00e9rienne constitue le standard empirique :<\/p>\n<p><code>generate double ln_price = ln(price)<\/code><br \/>\n<code>histogram ln_price, normal<\/code><br \/>\n<code>swilk ln_price<\/code><\/p>\n<p>L&rsquo;application du logarithme comprime l&rsquo;\u00e9ventail des grandes valeurs tout en \u00e9tirant les \u00e9carts entre les valeurs proches de z\u00e9ro, att\u00e9nuant m\u00e9caniquement l&rsquo;asym\u00e9trie positive. D&rsquo;autres fonctions analytiques sont fr\u00e9quemment mobilis\u00e9es selon l&rsquo;intensit\u00e9 de la d\u00e9viation constat\u00e9e :<\/p>\n<ul>\n<li><strong>Racine carr\u00e9e (<code>sqrt(x)<\/code>) :<\/strong> Transformation mod\u00e9r\u00e9e particuli\u00e8rement adapt\u00e9e aux variables de comptage ou aux asym\u00e9tries positives l\u00e9g\u00e8res.<\/li>\n<li><strong>Inverse (<code>1 \/ x<\/code>) :<\/strong> Transformation plus radicale destin\u00e9e aux asym\u00e9tries positives extr\u00eames, inversant l&rsquo;ordre de la s\u00e9rie (sauf utilisation de <span class=\"arabpsy-math-inline\" data-latex=\"-1\/x\">$-1\/x$<\/span>).<\/li>\n<li><strong>Puissance carr\u00e9e (<code>x^2<\/code>) ou cubique :<\/strong> Indiqu\u00e9e pour traiter les distributions pr\u00e9sentant une asym\u00e9trie n\u00e9gative \u00e9tir\u00e9e vers la gauche.<\/li>\n<\/ul>\n<p>Pour identifier la transformation optimale sans proc\u00e9der par t\u00e2tonnements, Stata propose l&rsquo;algorithme de <a rel=\"noopener noreferrer\" href=\"https:\/\/www.jstor.org\/stable\/2984418\">George Box et David Cox (1964)<\/a> via la commande <code>boxcox<\/code>. Cette approche estime par le maximum de vraisemblance le param\u00e8tre <span class=\"arabpsy-math-inline\" data-latex=\"lambda\">$lambda$<\/span> r\u00e9gissant la famille de transformations param\u00e9triques :<\/p>\n<p><code>boxcox price<\/code><\/p>\n<p>La valeur estim\u00e9e du param\u00e8tre d&rsquo;ajustement oriente le chercheur vers la fonction la plus appropri\u00e9e : un <span class=\"arabpsy-math-inline\" data-latex=\"lambda\">$lambda$<\/span> proche de z\u00e9ro justifie le recours au logarithme n\u00e9p\u00e9rien naturel, une valeur avoisinant 0,5 pr\u00e9conise la racine carr\u00e9e, tandis qu&rsquo;une valeur proche de 1 signale qu&rsquo;aucune transformation math\u00e9matique d&rsquo;\u00e9chelle n&rsquo;est n\u00e9cessaire.<\/p>\n<h3>11.2 M\u00e9thodes d&rsquo;estimation robustes et alternatives non param\u00e9triques<\/h3>\n<p>La transformation d&rsquo;une variable pr\u00e9sente l&rsquo;inconv\u00e9nient d&rsquo;alt\u00e9rer l&rsquo;\u00e9chelle de mesure originale, ce qui complique l&rsquo;interpr\u00e9tation directe des coefficients estim\u00e9s (passage \u00e0 des semi-\u00e9lasticit\u00e9s ou \u00e9lasticit\u00e9s pures). Lorsque l&rsquo;analyste souhaite conserver l&rsquo;\u00e9chelle d&rsquo;origine en pr\u00e9sence de r\u00e9sidus non normaux ou h\u00e9t\u00e9rosc\u00e9dastiques, l&rsquo;\u00e9conom\u00e9trie moderne privil\u00e9gie le recours \u00e0 des m\u00e9thodes d&rsquo;estimation robustes.<\/p>\n<p>Dans le cadre de la mod\u00e9lisation lin\u00e9aire par les moindres carr\u00e9s ordinaires, l&rsquo;adjonction de l&rsquo;option <code>vce(robust)<\/code> \u00e0 l&rsquo;instruction <code>regress<\/code> permet de mobiliser les estimateurs de variance asymptotique de White-Huber :<\/p>\n<p><code>regress price mpg weight i.foreign, vce(robust)<\/code><\/p>\n<p>Cette sp\u00e9cification n&rsquo;alt\u00e8re en rien la valeur ponctuelle des coefficients estim\u00e9s <span class=\"arabpsy-math-inline\" data-latex=\"\\hat{\\beta}\">$\\hat{\\beta}$<\/span>, mais recalcule la matrice de variance-covariance pour corriger l&rsquo;h\u00e9t\u00e9rosc\u00e9dasticit\u00e9 et mitiger l&rsquo;impact des d\u00e9viations de distribution des erreurs sur les erreurs-types des coefficients. La validit\u00e9 de cette correction repose toutefois sur l&rsquo;hypoth\u00e8se d&rsquo;\u00e9chantillons de taille suffisante pour assurer la convergence asymptotique.<\/p>\n<p>Dans l&rsquo;\u00e9ventualit\u00e9 o\u00f9 l&rsquo;analyste proc\u00e8de \u00e0 des comparaisons de localisation sur de tr\u00e8s petits effectifs (<span class=\"arabpsy-math-inline\" data-latex=\"N &lt; 30\">$N &lt; 30$<\/span>) o\u00f9 la normalit\u00e9 ne peut \u00eatre ni \u00e9tablie ni compens\u00e9e par les grands nombres, les tests d&#039;hypoth\u00e8ses param\u00e9triques usuels doivent \u00eatre remplac\u00e9s par leurs \u00e9quivalents non param\u00e9triques bas\u00e9s sur les rangs :<\/p>\n<ul>\n<li>Le test de Wilcoxon-Mann-Whitney remplace le test <span class=\"arabpsy-math-inline\" data-latex=\"t\">$t$<\/span> de Student pour deux \u00e9chantillons ind\u00e9pendants : <code>ranksum mpg, by(foreign)<\/code>.<\/li>\n<li>Le test des rangs sign\u00e9s de Wilcoxon remplace le test <span class=\"arabpsy-math-inline\" data-latex=\"t\">$t$<\/span> sur s\u00e9ries appari\u00e9es : <code>signrank pre_mesure = post_mesure<\/code>.<\/li>\n<li>Le test de Kruskal-Wallis remplace l&rsquo;ANOVA \u00e0 un facteur pour les comparaisons \u00e0 <span class=\"arabpsy-math-inline\" data-latex=\"k\">$k$<\/span> groupes : <code>kwallis price, by(rep78)<\/code>.<\/li>\n<\/ul>\n<p>Ces m\u00e9thodes \u00e9valuent les sommes de rangs sans imposer de distribution param\u00e9trique th\u00e9orique \u00e0 la population sous-jacente.<\/p>\n<h3>11.3 Inf\u00e9rence par r\u00e9\u00e9chantillonnage bootstrap dans Stata<\/h3>\n<p>L&rsquo;inf\u00e9rence par r\u00e9\u00e9chantillonnage bootstrap repr\u00e9sente une alternative moderne puissante lorsque les perturbations r\u00e9siduelles s&rsquo;\u00e9cartent manifestement du profil gaussien et que l&rsquo;\u00e9chantillon est trop modeste pour invoquer le th\u00e9or\u00e8me central limite avec confiance. Formalis\u00e9e par Bradley Efron en 1979, cette approche s&rsquo;abstient de formuler une hypoth\u00e8se th\u00e9orique pr\u00e9alable sur la loi d&rsquo;erreur : elle utilise l&rsquo;\u00e9chantillon observ\u00e9 lui-m\u00eame comme une approximation empirique de la population parente inobservable.<\/p>\n<p>La commande pr\u00e9fixe <code>bootstrap<\/code> sous Stata tire avec remise un nombre \u00e9lev\u00e9 <span class=\"arabpsy-math-inline\" data-latex=\"B\">$B$<\/span> de pseudo-\u00e9chantillons de m\u00eame taille <span class=\"arabpsy-math-inline\" data-latex=\"N\">$N$<\/span> que la base originale, r\u00e9estime les coefficients \u00e0 chaque r\u00e9plication, puis reconstruit la distribution empirique des param\u00e8tres d&rsquo;int\u00e9r\u00eat :<\/p>\n<p><code>bootstrap, reps(2000) seed(12345): regress price mpg weight i.foreign<\/code><\/p>\n<p>Ce processus g\u00e9n\u00e8re des intervalles de confiance non param\u00e9triques robustes aux non-normalit\u00e9s distributionnelles :<\/p>\n<ul>\n<li><code>estat bootstrap, normal<\/code> : Calcule l&rsquo;intervalle fond\u00e9 sur l&rsquo;erreur-type bootstrap standardis\u00e9e.<\/li>\n<li><code>estat bootstrap, percentile<\/code> : Calcule l&rsquo;intervalle bas\u00e9 directement sur les centiles de la distribution d&rsquo;\u00e9chantillonnage empirique r\u00e9pliqu\u00e9e.<\/li>\n<li><code>estat bootstrap, bc<\/code> : \u00c9tablit l&rsquo;intervalle de confiance corrig\u00e9 des biais d&rsquo;estimation (<em>bias-corrected<\/em>).<\/li>\n<li><code>estat bootstrap, bca<\/code> : Fournit la correction conjointe des biais et de l&rsquo;acc\u00e9l\u00e9ration (<em>bias-corrected and accelerated<\/em>), recommand\u00e9e pour absorber les asym\u00e9tries r\u00e9siduelles complexes.<\/li>\n<\/ul>\n<p>L&rsquo;inf\u00e9rence par bootstrap dispense l&rsquo;analyste de l&rsquo;hypoth\u00e8se de normalit\u00e9 des r\u00e9sidus pour d\u00e9river les intervalles de confiance et tester la significativit\u00e9 des estimateurs, sans imposer de modifications artificielles de l&rsquo;\u00e9chelle des variables d\u00e9pendantes.<\/p>\n<h2>12. Protocole syst\u00e9matique de d\u00e9cision et restitution acad\u00e9mique des r\u00e9sultats<\/h2>\n<h3>12.1 Arbre de d\u00e9cision int\u00e9gr\u00e9 pour l&rsquo;analyste de donn\u00e9es<\/h3>\n<p>Pour \u00e9viter les d\u00e9marches ad hoc ou l&rsquo;application d\u00e9sordonn\u00e9e de tests d&rsquo;ad\u00e9quation multiples, le diagnostic de normalit\u00e9 sous Stata gagne \u00e0 \u00eatre structur\u00e9 sous la forme d&rsquo;un protocole d&rsquo;analyse s\u00e9quentiel et standardis\u00e9. La d\u00e9marche d\u00e9cisionnelle pr\u00e9conis\u00e9e s&rsquo;articule autour des \u00e9tapes m\u00e9thodologiques suivantes :<\/p>\n<ul>\n<li><strong>Phase 1 : Examen graphique exploratoire :<\/strong> G\u00e9n\u00e9rer syst\u00e9matiquement un histogramme avec courbe normale superpos\u00e9e (<code>histogram, normal<\/code>) compl\u00e9t\u00e9 d&rsquo;un diagramme quantile-quantile (<code>qnorm<\/code>). Cette \u00e9tape permet d&rsquo;identifier la morphologie globale de la variable, la pr\u00e9sence de valeurs aberrantes isol\u00e9es et l&rsquo;\u00e9ventuelle bimodalit\u00e9 de la distribution.<\/li>\n<li><strong>Phase 2 : Calcul des moments empiriques :<\/strong> Ex\u00e9cuter l&rsquo;instruction <code>summarize, detail<\/code> pour extraire les valeurs de Skewness et de Kurtosis. V\u00e9rifier si les moments s&rsquo;inscrivent dans l&rsquo;intervalle empirique acceptable de mod\u00e9lisation lin\u00e9aire (<span class=\"arabpsy-math-inline\" data-latex=\"|\\text{Skewness}| &lt; 1\">$|\\text{Skewness}| &lt; 1$<\/span> et <span class=\"arabpsy-math-inline\" data-latex=\"\\text{Kurtosis} in [2, 4]\">$\\text{Kurtosis} in [2, 4]$<\/span>).<\/li>\n<li><strong>Phase 3 : \u00c9valuation inf\u00e9rentielle adapt\u00e9e \u00e0 l&rsquo;effectif :<\/strong>\n<ul>\n<li>Si <span class=\"arabpsy-math-inline\" data-latex=\"N &lt; 200\">$N &lt; 200$<\/span> : Ex\u00e9cuter le test de Shapiro-Wilk (<code>swilk<\/code>) en fixant le seuil alpha nominal \u00e0 0,05.<\/li>\n<li>Si <span class=\"arabpsy-math-inline\" data-latex=\"200 le N le 5,000\">$200 le N le 5,000$<\/span> : Ex\u00e9cuter le test de Shapiro-Francia (<code>sfrancia<\/code>) conjointement au test omnibus de D&rsquo;Agostino-Pearson (<code>sktest<\/code>).<\/li>\n<li>Si <span class=\"arabpsy-math-inline\" data-latex=\"N &gt; 5,000\">$N &gt; 5,000$<\/span> : Relativiser la significativit\u00e9 des p-valeurs formelles et prioriser l&rsquo;examen de l&rsquo;amplitude du coefficient <span class=\"arabpsy-math-inline\" data-latex=\"W&#039;\">$W&rsquo;$<\/span> et des moments standardis\u00e9s.<\/li>\n<\/ul>\n<\/li>\n<li><strong>Phase 4 : Choix de la strat\u00e9gie de traitement :<\/strong>\n<ul>\n<li>En cas de non-rejet de l&rsquo;hypoth\u00e8se nulle : Conserver les m\u00e9thodes param\u00e9triques classiques sous hypoth\u00e8se gaussienne.<\/li>\n<li>En pr\u00e9sence d&rsquo;asym\u00e9trie s\u00e9v\u00e8re sans changement de signe : Appliquer la transformation logarithmique ou la proc\u00e9dure de Box-Cox (<code>boxcox<\/code>).<\/li>\n<li>En pr\u00e9sence de leptokurticit\u00e9 mod\u00e9r\u00e9e sur grand \u00e9chantillon : Invoquer le th\u00e9or\u00e8me central limite et adjoindre l&rsquo;estimateur de variance robuste <code>vce(robust)<\/code>.<\/li>\n<li>Sur petit \u00e9chantillon en infraction distributionnelle s\u00e9v\u00e8re : Basculer vers l&rsquo;inf\u00e9rence par r\u00e9\u00e9chantillonnage bootstrap (<code>bootstrap<\/code>) ou les tests non param\u00e9triques sur les rangs.<\/li>\n<\/ul>\n<\/li>\n<\/ul>\n<h3>12.2 Bonnes pratiques de communication des diagnostics selon les normes scientifiques<\/h3>\n<p>La restitution acad\u00e9mique des diagnostics de normalit\u00e9 dans les sections m\u00e9thodologiques des m\u00e9moires, th\u00e8ses ou publications scientifiques \u00e9valu\u00e9es par les pairs doit respecter des standards d&rsquo;exhaustivit\u00e9 et de transparence statistique. Les formulismes vagues du type \u00ab la normalit\u00e9 a \u00e9t\u00e9 v\u00e9rifi\u00e9e graphiquement \u00bb doivent \u00eatre proscrits au profit d&rsquo;\u00e9nonc\u00e9s pr\u00e9cis d\u00e9taillant les m\u00e9triques calcul\u00e9es, les statistiques de test mobilis\u00e9es et les niveaux de probabilit\u00e9 associ\u00e9s.<\/p>\n<p>Voici un mod\u00e8le r\u00e9dactionnel recommand\u00e9 pour pr\u00e9senter l&rsquo;audit de normalit\u00e9 d&rsquo;une variable d\u00e9pendante continue :<\/p>\n<blockquote>\n<p>\u00ab La conformit\u00e9 de la distribution de la variable de consommation \u00e9nerg\u00e9tique (mpg) \u00e0 l&rsquo;hypoth\u00e8se de normalit\u00e9 a fait l&rsquo;objet d&rsquo;un audit distributionnel approfondi. L&rsquo;\u00e9valuation de la morphologie empirique met en \u00e9vidence un coefficient d&rsquo;asym\u00e9trie mod\u00e9r\u00e9 de 0,39 (erreur-type = 0,28) et un coefficient d&rsquo;aplatissement de 3,17 (erreur-type = 0,55), s&rsquo;inscrivant dans l&rsquo;intervalle d&rsquo;ellipticit\u00e9 standard. Le test d&rsquo;ad\u00e9quation de Shapiro-Wilk confirme le non-rejet de l&rsquo;hypoth\u00e8se nulle de normalit\u00e9 univari\u00e9e au sein de notre \u00e9chantillon (W = 0,981, V = 1,18, z = 0,36, p = 0,359). Ces r\u00e9sultats sont corrobor\u00e9s par le trac\u00e9 du diagramme quantile-quantile (Q-Q plot), attestant d&rsquo;un alignement uniforme des quantiles observ\u00e9s le long du vecteur th\u00e9orique diagonal, justifiant l&rsquo;application ult\u00e9rieure des mod\u00e8les param\u00e9triques lin\u00e9aires ordinaires. \u00bb<\/p>\n<\/blockquote>\n<p>Dans le cas d&rsquo;une analyse multivari\u00e9e ou d&rsquo;une r\u00e9gression complexe ayant n\u00e9cessit\u00e9 le recours \u00e0 une correction robuste, la justification m\u00e9thodologique doit expliciter le choix des estimateurs retenus :<\/p>\n<blockquote>\n<p>\u00ab L&rsquo;examen post-estimation des perturbations r\u00e9siduelles standardis\u00e9es r\u00e9v\u00e8le une violation de l&rsquo;hypoth\u00e8se de normalit\u00e9 (test omnibus de D&rsquo;Agostino-Pearson : Chi-deux ajust\u00e9(2) = 14,82, p &lt; 0,001 ; Mardia Kurtosis = 28,45, p = 0,004). En cons\u00e9quence, et pour pr\u00e9munir notre inf\u00e9rence contre les risques de sur-rejet d&#039;hypoth\u00e8se de type I li\u00e9s aux queues \u00e9paisses r\u00e9siduelles, l&#039;ensemble des erreurs-types associ\u00e9es aux estimateurs a \u00e9t\u00e9 calcul\u00e9 au moyen d&#039;un r\u00e9\u00e9chantillonnage bootstrap empirique avec correction d&#039;acc\u00e9l\u00e9ration et de biais (BCa) fond\u00e9 sur 2 000 r\u00e9plications al\u00e9atoires avec remise. \u00bb<\/p>\n<\/blockquote>\n<h3>12.3 Automatisation et reproductibilit\u00e9 des diagnostics sous Stata<\/h3>\n<p>Dans un contexte de recherche moderne exigeant une stricte reproductibilit\u00e9 des calculs, le diagnostic de normalit\u00e9 ne doit pas \u00eatre conduit de mani\u00e8re interactive et morcel\u00e9e au sein de la ligne de commande, mais formalis\u00e9 \u00e0 travers un script de programmation complet (do-file) structur\u00e9 et document\u00e9. L&rsquo;automatisation via des structures it\u00e9ratives permet de tester rapidement un ensemble de variables tout en centralisant les r\u00e9sultats au sein de matrices d&rsquo;exportation r\u00e9utilisables.<\/p>\n<p>Le bloc de code ci-apr\u00e8s montre l&rsquo;utilisation d&rsquo;une boucle it\u00e9rative <code>foreach<\/code> pour \u00e9valuer s\u00e9quentiellement une liste de variables continues et compiler leurs m\u00e9triques d&rsquo;asym\u00e9trie, d&rsquo;aplatissement et de Shapiro-Wilk :<\/p>\n<p><code>capture log close<\/code><br \/>\n<code>log using \"audit_normalite.log\", text replace<\/code><br \/>\n<code>sysuse auto, clear<\/code><br \/>\n<code>local var_audit price mpg weight displacement trunk<\/code><br \/>\n<code>display \"----------------------------------------------------------------------\"<\/code><br \/>\n<code>display \" VARIABLE       | SKEWNESS | KURTOSIS | W SHAPIRO | P-VALUE (SWILK) \"<\/code><br \/>\n<code>display \"----------------------------------------------------------------------\"<\/code><br \/>\n<code>foreach v of local var_audit {<\/code><br \/>\n<code>&nbsp;&nbsp;&nbsp;&nbsp;quietly summarize `v', detail<\/code><br \/>\n<code>&nbsp;&nbsp;&nbsp;&nbsp;local sk = r(skewness)<\/code><br \/>\n<code>&nbsp;&nbsp;&nbsp;&nbsp;local ku = r(kurtosis)<\/code><br \/>\n<code>&nbsp;&nbsp;&nbsp;&nbsp;quietly swilk `v'<\/code><br \/>\n<code>&nbsp;&nbsp;&nbsp;&nbsp;local w_stat = r(W)<\/code><br \/>\n<code>&nbsp;&nbsp;&nbsp;&nbsp;local p_stat = r(p)<\/code><br \/>\n<code>&nbsp;&nbsp;&nbsp;&nbsp;display as text %-15s \"`v'\" \" | \" as result %8.4f `sk' \" | \" %8.4f `ku' \" | \" %9.5f `w_stat' \" | \" %15.5f `p_stat'<\/code><br \/>\n<code>}<\/code><br \/>\n<code>display \"----------------------------------------------------------------------\"<\/code><br \/>\n<code>log close<\/code><\/p>\n<p>Ce script extrait dynamiquement les scalaires calcul\u00e9s stock\u00e9s en m\u00e9moire dans les structures syst\u00e8me <code>r()<\/code>, puis formate les r\u00e9sultats sous une table consolid\u00e9e claire. Cette approche \u00e9vite les saisies redondantes, garantit l&rsquo;auditabilit\u00e9 int\u00e9grale de la cha\u00eene de traitement et assure la tra\u00e7abilit\u00e9 des diagnostics pr\u00e9alables \u00e0 la mod\u00e9lisation statistique d\u00e9finitive.<\/p>\n<h2>R\u00e9f\u00e9rences<\/h2>\n<ul>\n<li>Box, G. E. P., &amp; Cox, D. R. (1964). An analysis of transformations. <em>Journal of the Royal Statistical Society: Series B (Methodological)<\/em>, 26(2), 211\u2013243. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1111\/j.2517-6161.1964.tb00553.x\">https:\/\/doi.org\/10.1111\/j.2517-6161.1964.tb00553.x<\/a><\/li>\n<li>D&rsquo;Agostino, R. B., &amp; Pearson, E. S. (1973). Tests for departure from normality. Empirical results for the distributions of <span class=\"arabpsy-math-inline\" data-latex=\"b_2\">$b_2$<\/span> and <span class=\"arabpsy-math-inline\" data-latex=\"\\sqrt{b_1}\">$\\sqrt{b_1}$<\/span>. <em>Biometrika<\/em>, 60(3), 613\u2013622. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1093\/biomet\/60.3.613\">https:\/\/doi.org\/10.1093\/biomet\/60.3.613<\/a><\/li>\n<li>Doornik, J. A., &amp; Hansen, H. (2008). An omnibus test for univariate and multivariate normality. <em>Oxford Bulletin of Economics and Statistics<\/em>, 70(s1), 927\u2013939. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1111\/j.1468-0084.2008.00537.x\">https:\/\/doi.org\/10.1111\/j.1468-0084.2008.00537.x<\/a><\/li>\n<li>Efron, B. (1979). Bootstrap methods: Another look at the jackknife. <em>The Annals of Statistics<\/em>, 7(1), 1\u201326. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1214\/aos\/1176344552\">https:\/\/doi.org\/10.1214\/aos\/1176344552<\/a><\/li>\n<li>Henze, N., &amp; Zirkler, B. (1990). A class of invariant consistent tests for multivariate normality. <em>Communications in Statistics &#8211; Theory and Methods<\/em>, 19(10), 3595\u20133617. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1080\/03610929008830400\">https:\/\/doi.org\/10.1080\/03610929008830400<\/a><\/li>\n<li>Lilliefors, H. W. (1967). On the Kolmogorov-Smirnov test for normality with mean and variance unknown. <em>Journal of the American Statistical Association<\/em>, 62(318), 399\u2013402. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1080\/01621459.1967.10482916\">https:\/\/doi.org\/10.1080\/01621459.1967.10482916<\/a><\/li>\n<li>Mardia, K. V. (1970). Measures of multivariate skewness and kurtosis with applications. <em>Biometrika<\/em>, 57(3), 519\u2013530. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1093\/biomet\/57.3.519\">https:\/\/doi.org\/10.1093\/biomet\/57.3.519<\/a><\/li>\n<li>Royston, P. (1982). An extension of Shapiro and Wilk&rsquo;s W test for normality to large samples. <em>Journal of the Royal Statistical Society: Series C (Applied Statistics)<\/em>, 31(2), 115\u2013124. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.2307\/2347973\">https:\/\/doi.org\/10.2307\/2347973<\/a><\/li>\n<li>Royston, P. (1991). Approximating the Shapiro-Wilk W-test for non-normality. <em>Statistics and Computing<\/em>, 2(3), 117\u2013119. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1007\/BF01891203\">https:\/\/doi.org\/10.1007\/BF01891203<\/a><\/li>\n<li>Shapiro, S. S., &amp; Francia, R. S. (1972). An approximate analysis of variance test for normality. <em>Journal of the American Statistical Association<\/em>, 67(337), 215\u2013216. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1080\/01621459.1972.10481232\">https:\/\/doi.org\/10.1080\/01621459.1972.10481232<\/a><\/li>\n<li>Shapiro, S. S., &amp; Wilk, M. B. (1965). An analysis of variance test for normality (complete samples). <em>Biometrika<\/em>, 52(3\/4), 591\u2013611. <a rel=\"noopener noreferrer\" href=\"https:\/\/doi.org\/10.1093\/biomet\/52.3-4.591\">https:\/\/doi.org\/10.1093\/biomet\/52.3-4.591<\/a><\/li>\n<li>StataCorp. (2023). <em>Stata Base Reference Manual: Release 18<\/em>. Stata Press. <a rel=\"noopener noreferrer\" href=\"https:\/\/www.stata.com\/manuals\/r.pdf\">https:\/\/www.stata.com\/manuals\/r.pdf<\/a><\/li>\n<\/ul>\n","protected":false},"excerpt":{"rendered":"<p>Guide m\u00e9thodologique complet pour \u00e9valuer la normalit\u00e9 univari\u00e9e et multivari\u00e9e dans Stata : tests statistiques, approches graphiques et rem\u00e9diations.<\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"open","ping_status":"closed","template":"","categories":[2244,2269],"tags":[8077,4660,2216,3292,2779],"class_list":["post-6447","psych_stat","type-psych_stat","status-publish","hentry","category-logiciel-stata","category-methodes-quantitatives","tag-asymetrie-et-kurtosis","tag-shapiro-wilk","tag-stata","tag-statistiques-parametriques","tag-test-de-normalite"],"_links":{"self":[{"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=\/wp\/v2\/psych_stat\/6447","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=\/wp\/v2\/psych_stat"}],"about":[{"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=\/wp\/v2\/types\/psych_stat"}],"author":[{"embeddable":true,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Fcomments&post=6447"}],"version-history":[{"count":1,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=\/wp\/v2\/psych_stat\/6447\/revisions"}],"predecessor-version":[{"id":6454,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=\/wp\/v2\/psych_stat\/6447\/revisions\/6454"}],"wp:attachment":[{"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=6447"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=6447"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/french.arabpsychology.com\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=6447"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}