Discussion :

[Rawley]

Bonjour,

J’ai un jeu de données dont j’ai représenté l’histogramme. Désormais je veux ajouter à mon histogramme une distribution de loi normale pour voir si mes données suivent ce type de loi.

Je sais faire les deux graphiques indépendamment avec les fonctions hist et histfit mais je n’arrive pas à faire un seul graphique. Il faudrait que ma courbe rouge du graphique 2 soit sur mon histogramme du graphique 1 tout en conservant l’axe des abscisses du graphique 1…

Est-ce que quelqu’un a déjà confronté à ce genre de chose ?

Raw

[ol9245]

entre tes deux commandes de graphique, insérer la commande suivante :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

hold on

[shaiHulud]

bonjour,
N'oublie pas de renormaliser l'histogramme, de sorte que son aire soit 1

[Rawley]

Bonjour,

"Hold on" ça va pas de problème par contre il doit sans doute falloir que je normalise ma courbe de Gauss comment je fais ça ?

Merci pour votre réponse !

Raw

[shaiHulud]

Ce n'est pas la courbe de Gauss qu'il faut normaliser, mais bien l'histogramme.
La fonction hist() est censée avoir une options 'normalize', mais à la main ca marche aussi: par défaut, l'aire sous un histogramme régulier de pas dx vaut N*dx.

[Rawley]

Re-bonjour,

J'ai tenté une autre option : j'ai fait un histogramme de mes données, j'ai calculé la moyenne et l'écart-type et le vecteur x qui représente l'équation d'une courbe de Gauss :

x = (1./(ecart_type.*sqrt(2.*pi())))*(exp((-1/2).*((data-moy_data)./(ecart_typea)).^2));

Je fais un plot avec l'histogramme et la courbe x par dessus, j'obtiens la bonne forme mais les amplitudes de la courbe sont fausses. J'obtiens environ 7 en effectif pour le sommet de ma cloche de Gauss (voir graph ci-joint) alors que lorsque je fais histfit j'obtiens un sommet vers 8 en effectif (voir graph quelque posts plus haut)...

Est-ce que quelqu'un saurait où est mon ereur ?

En vour remerciant,

Raw

[shaiHulud]

Le différence est surement due à l'arrondie de l'histogramme (puisqu'il calcul une approximation en escalier).

Change les bins (nombre ou centres) de l'histogramme pour le vérifier