Histogramme

La première phase est la collecte des données en cours de fabrication. Cette collecte peut être réalisée soit de façon exceptionnelle à l’occasion de l’étude de la variable soit en utilisant un relevé automatique ou manuel fait lors d’un contrôle réalisé dans le cadre de la surveillance du procédé de fabrication.

Sans qu’il soit réellement possible de donner un nombre minimum, il faut que le nombre de valeurs relevées soit suffisant. Plus on dispose d’un nombre élevé de valeurs, plus l’interprétation sera aisée.

Le choix des classes, soit leur nombre et leurs largeurs, n'est pas univoque. Il convient pour les déterminer de prendre en compte à la fois la nature de la distribution et le nombre de points de données. Souvent, dans le cadre d’une analyse de ce type, on utilise des classes de largeur identique.

On pourra trouver dans la littérature de nombreuses suggestions de choix pour le nombre de classe. Citons par exemple :

• Celle de Herbert Sturges (1926)[1] qui, pour N points de données répartis avec une distribution approximativement normale, suggère un nombre de classes K obtenu avec la formule suivante [2] :

${\displaystyle K=1+\log _{2}N\approx 1+{\frac {10}{3}}\log _{10}N}$

On pourra consulter utilement l'article Règle de Sturges à ce propos.

• L'alternative à la règle précédente est la règle dite de Rule [3] où

${\displaystyle K=2\cdot N^{\frac {1}{3}}}$

• Le choix simple de la racine carrée :

${\displaystyle K={\sqrt {N}}}$

En tout état de cause, l’histogramme étant un outil visuel, il est possible de faire varier le nombre de classes. Ceci permet de voir l’histogramme avec un nombre différent de classes et ainsi de trouver le meilleur compromis qui facilitera l’interprétation. L’utilisation d’un logiciel dédié ou plus simplement d’un tableur facilite cette opération.

L’amplitude (minimale) w de l’histogramme est

${\displaystyle w={\hbox{valeur maximale}}-{\hbox{valeur minimale}}\,}$

Il peut cependant être intéressant pour obtenir un histogramme plus parlant de choisir une amplitude plus large que l'amplitude minimale.

L’amplitude h théorique de chaque classe est alors :

${\displaystyle h={\frac {w}{K}}}$

Il faut arrondir cette valeur à un multiple de résolution de l’instrument de mesure (arrondi à l'excès).

Exemple : Soit la masse d’une préparation culinaire avant conditionnement. Le calcul d'amplitude de classe donne h_th = 0,014 kg. La résolution de la balance utilisée est de 0,001 kg. On arrondit la valeur h à 0,015 kg.

Les classes peuvent être du type [limite _inférieure ; limite _supérieure[ ou ] limite _inférieure ; limite _supérieure].

La valeur minimale de la première classe est donnée par la valeur minimale de la série moins une demi-résolution.

Exemple : la valeur la plus petite relevée lors de la fabrication de la préparation culinaire est de 0,498 kg. La limite inférieure sera : 0,498 – (0,001 / 2) = 0,497 5 kg.

Pour plus de facilité, il est préférable de prendre une valeurs « ronde » par exemple 0,495 kg

• Pour un exemple développé in extenso, voir article: Statistiques élémentaires continues

Plusieurs choix sont possibles pour déterminer la hauteur des rectangles.

1. Les hauteurs correspondent aux fréquences absolues, soit pour chaque rectangle le nombre d'observations dans la classe correspondante.
2. Les hauteurs correspondent aux fréquences relatives, soit pour chaque rectangle la proportion, donnée par exemple en pourcentage, d'observations dans la classe correspondante.
3. Les hauteurs sont déterminées de manière que la surface du rectangle corresponde à la fréquence relative d'observations dans la classe correspondante.

La troisième méthode permet en outre d'accommoder des classes de largeurs variables ce qui est utile lorsque les données sont peu denses dans certaines régions comme dans les queues de distribution.

Soit la fabrication de rations alimentaires, la pesée des rations avant emballage donne la série de mesures suivantes en kg :

Les caractéristiques du relevé sont les suivantes :

• Le nombre d'échantillons : N=64
• L'étendue : w=0,098 kg
• Valeur minimale : 0,498 kg
• Valeur maximale : 0,596 kg

On en déduit les paramètres suivants pour l'histogramme :

• Le nombre de classes est de 7 (en utilisant la formule avec le logarithme)
• L'amplitude de classe est 0,098/7 = 0,014 kg que l'on arrondit à 0,015 kg (résolution de la balance : 0,001 kg)
• La valeur minimale de la première classe est de 0,498 – (0,001/2) = 0,4975. Par souci de facilité pour l'interprétation, on peut arrondir cette valeur à 0,495 kg.

On obtient l'histogramme suivant :

Comparaison d'un histogramme avec la courbe d'une loi normale.

La distribution de beaucoup de paramètres industriels correspond souvent à une loi normale. On compare souvent l'histogramme obtenu au profil « en cloche » de la loi normale. Cette comparaison est visuelle et même si elle peut être une première approche, elle ne constitue pas un test de « normalité ». Pour cela, il faut exécuter un test dont un des plus classiques est la droite de Henry.

La distribution suivant la loi normale, si elle est extrêmement fréquente, n'est pas systématique. On vérifiera que la distribution ne correspond pas à une distribution de défaut de forme (exemple : mesure de l'excentration dans un tube, position d'objets lancés dans la direction d'un mur dont certains rebondissent sur ce mur).

L'interprétation peut, par exemple, donner les résultats suivants :

Histogramme montrant un mélange de deux lots.	Histogramme montrant un mélange de deux lots mais avec une moyenne proche. On veillera dans ce cas à faire aussi varier le nombre de classes pour vérifier qu'il ne s'agit pas d'un problème de construction.	Histogramme montrant que le lot a subi un tri. Tous les éléments pour lesquels la valeur du paramètre mesuré était inférieure à A ont été supprimés.

Dans le cas d'histogramme montrant un mélange de deux lots ayant une moyenne différente, il existe des cas où la dispersion présente cet aspect sans pour autant incriminer un mélange. C'est par exemple le cas de la mesure d'une pièce cylindrique mais qui présente un défaut de type ovalisation. Les deux moyennes représentent alors le grand diamètre et le petit diamètre. C'est la connaissance du procédé et/ou du produit qui permet de réaliser ce type d'interprétation.

On recherche les probabilités ${\displaystyle p_{i}}$ qui caractérisent la distribution. On note cette distribution f par abus. Un estimateur naturel est:

${\displaystyle {\widehat {f}}(x)={\frac {n(x)}{n}}}$

où ${\displaystyle n(x)}$ est le nombre d'observations de l'échantillon qui sont égales à x. Une manière alternative de noter cet estimateur est:

${\displaystyle {\widehat {f}}(x)={\frac {1}{n}}\sum _{i}^{n}I(x_{i}=x)}$

où ${\displaystyle I(\cdot )}$ est la fonction indicatrice : elle vaut 1 lorsque son argument est vrai.

L'estimateur précédent n'est plus valable, car dans le cas continu, on ne peut plus compter le nombre d'observations exactement égales à x. Par contre, on peut considérer une boîte centrée en x, et de largeur h, paramètre positif. On peut compter le nombre d'observations approximativement (et non plus exactement) égales à x, en comptant les observations tombant dans ladite boîte. L'estimateur[4] devient:

${\displaystyle {\widehat {f}}(x)={\frac {1}{nh}}\sum _{i}^{n}I\left(-{\frac {1}{2}}\leq {\frac {x_{i}-x}{h}}\leq {\frac {1}{2}}\right)}$

ou encore, en posant ${\displaystyle y_{i}=(x_{i}-x)/h}$ :

${\displaystyle {\widehat {f}}(x)={\frac {1}{nh}}\sum _{i}^{n}I\left(|y_{i}|\leq {\frac {1}{2}}\right)}$.

Le paramètre h contrôle le niveau de lissage de l'estimation et doit être recherché avec soin. L'estimateur précédent présente de bonnes propriétés comparables à celles d'une densité continue:

• il est positif ;
• il s'intègre à l'unité.

Toutefois, il présente un gros défaut pour pouvoir estimer une densité: il n'est pas continu. Pour gagner la continuité, on utilisera l'estimateur de Parzen (ou à noyau). Le principe est simple: il suffit de remplacer la fonction indicatrice par une fonction réelle, qui attribue un poids d'autant plus important que les observations sont situées à proximité de x. À titre d'exemple voici deux estimations de densité par la méthode des noyaux avec des paramètres différents :

• Maurice Lethielleux, Statistique descriptive, éditions Dunod, Paris, 1999 (ISBN 2 10 003513 4), 124 pages.
• Maurice Pillet, Appliquer la maîtrise statistique des procédés MSP/SPC, Les Éditions d'Organisation, 1995 (ISBN 2-7081-1774-2), 336 pages.
• Pierre Souvay, Statistique et qualité, AFNOR, Paris, collection « À savoir », 1994, 40 pages.
• Pierre Souvay, Savoir utiliser la statistique, outil à la décision et à l'amélioration de la qualité, AFNOR, Saint-Denis-la-Plaine, 2002 (ISBN 2-12-475821-7), 434 pages

• Représentations graphiques de données statistiques
• Centile
• Décile
• Quantile
• Quartile
• Histogramme conjoint

• A Method for Selecting the Bin Size of a Histogram
• cours avec exemples de classes de largeurs distinctes