Discussion :

[Ekhron]

Bonjour


Je me suis mis à R il y a très peu de temps et j'ai pas desx difficultés pour faire des graphiques.

Mon problèmes est que je veux représenter l'approximation de pi, de la même manière que l'on représente la loi des grands nombres. C'est à dire sur un graphique, représenter plusieurs courbes qui prennent différentes valeurs et qui s'approche de la valeur Pi, pour N suffisement grand

J'ai déjà représenté une simulation pour approximer la valeur de pi avec un carré et un quart de cercle.

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
n=1000
x<-seq(0,1,0.0001)
XY=matrix(runif(2*n),ncol=2)
p=mean(XY[,1]^2+XY[,2]^2<1)
pi=4*p
pi
plot(XY[,1],XY[,2],axes=FALSE,pch=1,main="Estimation de pi",xlab="X",ylab="Y");axis(1);axis(2)
polygon(c(x,0),c(sqrt(1-x^2),0),density=10,col="red")

Je voudrais obtenir cinq courbes différentes (ou trois), de couleurs différentes, pour représenter l'approximation de pi, pour un nombre N assez grand, mais j'y arrive pas.
Je sais simuler plusieurs approximation à la suite

i<-1
x<-0
while(i<=5){
n=1000
XY=matrix(runif(2*n),ncol=2)
p=mean(XY[,1]^2+XY[,2]^2<1)
pi=4*p
x[i]<-pi
print(pi)
i<-i+1
}

mais pour faire un graphe, là j'y arrive pas.


Merci d'avance pour votre aide.

[tomaprice]

c'est pas cela que tu cherches à faire ?

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
pi=vector()
for(n in 1:1000){
	XY=matrix(runif(2*n),ncol=2)
	p=mean(XY[,1]^2+XY[,2]^2<1)	
	pi[n]=4*p
}
plot(pi,type="l")

en tout cas il s'agit bien d'une représentation de la loi des grands nombres

En statistiques, la loi des grands nombres exprime le fait que les caractéristiques d'un échantillon aléatoire se rapprochent d'autant plus des caractéristiques statistiques de la population que la taille de l'échantillon augmente @wiki