Deviation nombres aleatoires

[Myriam Croze]

Bonjour,

Je travail sur R et je n'arrive pas a fair un graphe.
Je veux tirer au sort aleatoirement des valeurs d'une table de donnees.
Et je veux regarder a partir de quand la variance se stabilise.
C'est a dire, si je tire 10, 20, 30, ......valeurs, a partir de quand j'aurais la meme variation.
Pour cela, je voudrais faire plusieurs boxplot pour 10, 20, 30, ...., 100 valeurs tirees au sort aleatoirement.
Voila le script que j'ai fait jusqu'a present:

Code : Tout sélectionner

 # read in the data
S_table<-read.table("D:/PhD/HKA/Fen_2000/Input_MLHKA_FR_Chr3L_neut.csv", sep=",", header=TRUE)
# store S data in a vector
S<-S_table$S

boxplot (sample(S,10), sample(S,20), sample(S,30), sample(S,40), sample(S,50), sample(S,60), sample(S,80), sample(S,90), sample(S,100))


Je suppose que je dois faire une boucle pour repeter plusieurs fois le tirage au sort et avoir une la variation.
En effet je m'attend a une variation beaucoup plus grande pour 10 valeurs tirees aleatoirement et qui diminue avec l'augmentation du nombre de valeurs dans mon jeu de donnees tire au sort.

Helas, je ne vois pas comment faire et quoi ajouter dans le script.
Si quelqu'un pouvez m'aider, je lui en serai vraiment reconnaissante.
Merci.

Myriam

[Perceval Sondag]

Cela dépend de la variation que vous voulez observer. Variation sur les moyennes ? Les étendues ?

Si vous voulez observer une variation sur le boxplot 'complet' vous auriez plus intérêt à observer des graphiques comme ceci :

Code : Tout sélectionner

par(mfrow= c(3,4))
for (i in seq(10,100,by=10)) {
boxplot (sample(S,i), sample(S,i), sample(S,i), sample(S,i), sample(S,i), sample(S,i), sample(S,i), sample(S,i), sample(S,i), sample(S,i))
}


En espérant que ça aide,

Perceval

[Myriam Croze]

Merci pour ta reponse perceval. Mais ce n'est pas exactement ce que je veux.
En fait ce que je voudrais, c'est sur un seul graphes, avoir plusieurs boxplot pour differentes valeurs de tirage au sort.
Et je voudrais voir la variation par rapport a la moyenne.
En gros je voudrais savoir a partir de quel moment, on a une variation qui se stabilise et reste la meme.

Myriam

[Perceval Sondag]

je voudrais voir la variation par rapport a la moyenne

Dans ce cas :

Code : Tout sélectionner

Smean <- matrix (0,1000,10)
size <- seq(10,100,by=10)
colnames(Smean)<-size

for (j in  1:10){
   for (i in 1:1000){
      Smean[i,j]<-mean(sample(S,size[j]))

   }
}

boxplot(Smean)


[Perceval Sondag]

Ou encore :

Code : Tout sélectionner


Smean2 <- matrix (0,1000,100)
Svar <- numeric(100)

for (j in  1:100){
   for (i in 1:1000){
      Smean2[i,j]<-mean(sample(S,j))
   
   }
   Svar[j]<-var(Smean2[,j])
}

plot(Svar, type = "l")


[Myriam Croze]

Merci pour ton aide.
J'ai reussi a faire mon graphe.