[R] Simuler une loi des grands nombres

Question

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.

Hxyp · Answer

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))

[R] Simuler une loi des grands nombres

1 réponse

Votre réponse

Discussions similaires

Newsletters