graphe.desimple<-function(n,translog=FALSE){ # simulation de n tirages successifs de dès # pour i compris entre 1 et 6, sont tracées six courbes, # chacune d'entre elles représentent les proportions d'apparition du numéro i, au cours des lancers. # si translog=TRUE, on représente log10(proportion) en fonction de log10(numéro de lancers) # ********************************************************* # graphe.desimple(n,translog=FALSE) : # * Variables d'entrées : # * n : nombre de tirages # * translog : booléen égal à FALSE par défaut # * Variables de sortie : # * pour 1 <=i<=6, les graphes (j,proportion(apparition numéro i)) où 1<=j<=n # 2009 by Jérôme BASTIEN # Université Claude Bernard Lyon I, UFRSTAPS, Laboratoire CRIS, Villeurbanne # E-Mail : jerome.bastien@univ-lyon1.fr echan<-sample(1:6,size=n,replace=T) tata<-matrix(nrow =6,ncol=n) for(i in 1:6) {auxi<-cumsum(echan==i) tata[i,]<-auxi/(1:n) if (translog) tata[i,]<-log10(abs(tata[i,]-1/6))} x<-if (translog) log10(1:n) else 1:n return({if (translog) plot(x,tata[1,],col=2,type='l',xlab="log10(n)",ylab="log10(|proportion-1/6|)") else plot(x,tata[1,],col=2,type='l',xlab="n",ylab="proportion",ylim=c(0,1/3)) for(i in 2:6) points(x,tata[i,],col=i+1,type='l') if (!translog) points(c(1,n),c(1/6,1/6),type='l',col=1) }) }