Forum de Statistiques
Vous souhaitez réagir à ce message ? Créez un compte en quelques clics ou connectez-vous pour continuer.
Les posteurs les plus actifs de la semaine
Aucun utilisateur

Le Deal du moment : -15%
(Adhérents) LEGO® Icons 10318 Le Concorde
Voir le deal
169.99 €

Code.

2 participants

Aller en bas

Code. Empty Code.

Message par FMarwen Ven 31 Jan 2014 - 14:54

Bonjour,

Je voudrais gérer 100 échantillons de taille 100 de loi normale et tracer l'histogramme des moyennes des échantillons et la densité correspondante, j'ai arrivé à faire le code suivant sauf que l'axe d'ordonné donne des valeurs >1 et moi je veux des probabilité ! Quelqu'un parmi vous peux me aider ? Merci d'avance.
Code:

Hist_Moy=function(m,s2){
data<-matrix(rnorm(10000,mean=m,sd=s2), nrow=100)
vect<-apply(data,1,mean)
par(mfrow=c(1,2))
{hist(vect,prob=T, breaks="FD", main="Histogramme_Moy de 100 tirages N(m,s2)")}
den <-density(vect)
plot(den, col = "red",add=T)
abline(v = m, col = "green", lty ="dotted")
}
Hist_Moy(0,1)

Fichiers joints
Code. Attachment
a.jpeg Vous n'avez pas la permission de télécharger les fichiers joints.(25 Ko) Téléchargé 0 fois

FMarwen

Nombre de messages : 226
Date d'inscription : 17/04/2013

Revenir en haut Aller en bas

Code. Empty Re: Code.

Message par droopy Mer 5 Fév 2014 - 8:34

Bonjour,

le code est bon tu pourrais juste remplacer la ligne vect <- apply(data,1,mean) par vect <- colMeans(data).
Ici il est question d'une fonction de densité, donc il est possible que les valeurs en ordonnées dépassent 1, puisque c'est l'aire sous la courbe qui doit faire 1. Ici tu as des valeurs entre -0.4 et 0.4 en gros, donc même si tu avais une répartition uniforme des valeurs, tu aurais une fonction de densité qui serait un rectangle compris en -0.4 et 0.4 soit une longueur de 0.8 en abscisse. Pour que l'aire de ton rectangle fasse 1 il te faudrait que l'autre côté du rectangle fasse 1/0.8 = 1.25 donc > 1. Ici tu n'as pas de rectangle puisque tes valeurs suivent une loi normale donc, il est d'autant plus normal que les valeurs sur les ordonnées dépassent 1. Les valeurs en ordonnées ne sont pas des probabilités.

Cdlt
droopy
droopy

Nombre de messages : 1156
Date d'inscription : 04/09/2009

Revenir en haut Aller en bas

Code. Empty Re: Code.

Message par FMarwen Ven 7 Fév 2014 - 17:53

Merci beaucoup pour l'explication c'est très claire.

FMarwen

Nombre de messages : 226
Date d'inscription : 17/04/2013

Revenir en haut Aller en bas

Code. Empty Re: Code.

Message par Contenu sponsorisé


Contenu sponsorisé


Revenir en haut Aller en bas

Revenir en haut

- Sujets similaires

 
Permission de ce forum:
Vous ne pouvez pas répondre aux sujets dans ce forum