LPNHE - UMR 7585 - calcul d’histogramme

Aller au menu
Aller à la recherche
Aller au contenu
Aller à la navigation

Partenaires

Rechercher

Accueil > Le LPNHE > Pages du personnel permanent > Beau Tristan > Divers > Vielleries > calcul d’histogramme

calcul d’histogramme

23 novembre 2005

par Tristan Beau - 7 juin 2017

Voici quelques scripts permettant, à partir d’une liste de nombre, d’effectuer un calcul d’histogramme.

script de base : `hist`

Voici le script hist qui effectue le calcul de base :

#!/bin/sh

awk 'BEGIN 
  {min=0; max=100; N=15;'$1';delta=(max-min)*1./(N); for(i=0;i<N;i++) a[i]=0;} 
  {idx=int(($1-min)/delta); if($1>=min && idx<N) a[idx]++; } END {for(i=0;i<N;i++) 
  { val=a[i]; print (min+i*delta+delta/2)" "val;} }' 
  < /dev/stdin

attention, la commande awk a besoin d’être écrite en une unique ligne, le copier-coller peut ne pas fonctionner correctement...

On lancera le script par exemple :

hist < liste_de_nombre.txt > histo.txt

Le fichier d’exemple liste_de_nombre.txt est donné ci-contre. On pourra ensuite créer l’histogramme graphique par exemple avec gnuplot avec la commande :

echo 'set terminal png; plot "histo.txt" with boxes' | gnuplot > histo.png

qui créera donc l’image au format png nommée histo.png, image que l’on visualisera avec l’outil que l’on préfère...

Par défaut, le binning s’effectuera entre 0 et 100, avec 15 bins. Pour changer ce fonctionnement, on pourra utiliser une option à hist.
Ainsi, pour binner entre 1 et 2 par exemple, sur 25 bins, on fera :

hist "min=10;max=20;N=25" < liste_de_nombre.txt > histo2.txt
echo 'set terminal png; plot "histo2.txt" with boxes' | gnuplot > histo2.png

histogramme avec bornes automatiques : `histauto`

Pour que le choix des bornes soit automatiques, il faut trouver les valeurs extrémales des données. On utilisera une petite fonction awk pour cela. Soit :

#!/bin/sh

TMP=`mktemp -p /tmp histauto.XXXXXXXX`

hist `cat /dev/stdin | tee $TMP | awk 'BEGIN {init=0} {if(init==0) { min=$1; max=$1;init=1;} else { if ($1>max) max=$1; if ($1<min) min=$1;} } END {print "min="min";max="max""}' ` < $TMP

rm -f $TMP

et l’on fera attention à ce que la commande awk se trouve également en une ligne unique.

On appelera ce script avec :

histauto < liste_de_nombre.txt > histo3.txt
echo 'set terminal png; plot "histo3.txt" with boxes' | gnuplot > histo3.png

Dans la même rubrique :

Partenaires

Menu

Rechercher

Navigation

calcul d’histogramme

script de base : `hist`

histogramme avec bornes automatiques : `histauto`

Partenaires

Menu

Rechercher

Navigation

calcul d’histogramme

script de base : hist

histogramme avec bornes automatiques : histauto

script de base : `hist`

histogramme avec bornes automatiques : `histauto`