Les posteurs les plus actifs de la semaine
Aytan
Code. Vote_lcapCode. Voting_barCode. Vote_rcap 
c@ssoulet
Code. Vote_lcapCode. Voting_barCode. Vote_rcap 
Adddd
Code. Vote_lcapCode. Voting_barCode. Vote_rcap 
Eric Wajnberg
Code. Vote_lcapCode. Voting_barCode. Vote_rcap 
gg
Code. Vote_lcapCode. Voting_barCode. Vote_rcap 


Code.

Aller en bas

Code. Empty Code.

Message par FMarwen le 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

Voir le profil de l'utilisateur

Revenir en haut Aller en bas

Code. Empty Re: Code.

Message par droopy le 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 : 1126
Date d'inscription : 04/09/2009

Voir le profil de l'utilisateur

Revenir en haut Aller en bas

Code. Empty Re: Code.

Message par FMarwen le 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

Voir le profil de l'utilisateur

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