Les posteurs les plus actifs de la semaine
Aucun utilisateur |
Sujets les plus vus
Code.
2 participants
Page 1 sur 1
Code.
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.
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
FMarwen- Nombre de messages : 226
Date d'inscription : 17/04/2013
Re: Code.
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
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- Nombre de messages : 1156
Date d'inscription : 04/09/2009
Re: Code.
Merci beaucoup pour l'explication c'est très claire.
FMarwen- Nombre de messages : 226
Date d'inscription : 17/04/2013
Page 1 sur 1
Permission de ce forum:
Vous ne pouvez pas répondre aux sujets dans ce forum
|
|