[R] Simuler une loi des grands nombres

Fermé
Keldras Messages postés 1 Date d'inscription dimanche 29 janvier 2012 Statut Membre Dernière intervention 29 janvier 2012 - 29 janv. 2012 à 12:58
Hxyp Messages postés 401 Date d'inscription vendredi 28 janvier 2011 Statut Membre Dernière intervention 27 avril 2014 - 3 févr. 2012 à 14:02
Bonjour,

Je voudrais simuler une lois des grands nombres autour de la valeur pi.

Je voudrais en simuler plusieurs et les mettre sur le même graphique.

Le problème c'est que je n'arrive même pas à en simuler une seule.

J'ai fait ce code:

n=0
while (n<1000)
{
XY=matrix(runif(2*n),ncol=2)
p=mean(XY[,1]^2+XY[,2]^2<1)
pi=4*p
plot(n,pi,ylim=c(2,4))
n<-n+1
}


Dans le fond ça marche mais ça fait juste un oscilloscope, il affiche la valeur, il passe à l'autre, mais sans l'imprimer sur le graphique et permettre de voir la courbe finale.


Merci d'avance pour votre aide.

1 réponse

Hxyp Messages postés 401 Date d'inscription vendredi 28 janvier 2011 Statut Membre Dernière intervention 27 avril 2014 54
3 févr. 2012 à 14:02
Bonjour,
En mettant les valeurs dans un tableau et l'affichant à la fin :
n=0
Pi<-array(0,c(1000))
while (n<1000)
{
XY=matrix(runif(2*n),ncol=2)
p=mean(XY[,1]^2+XY[,2]^2<1)
Pi[n]=4*p
n<-n+1
}
plot(Pi,ylim=c(2,4))
0