Discussion :

[wanaga3166]

Bonjour à tous les expeRts,

j'ai une liste de fichiers, et je souhaite faire un histogramme pour chacun de ces fichiers. J'ai utilisé la fonction list.files, mais elle ne me fait qu'un seul histogramme au lieu de tous.

J'ai commencé par faire le list.files, puis j'essaye de dire que tel fichier tu vas me faire l'histogramme.

Voici le code:

Code :

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

1
2
3
4
5
6
7
8
9
10
11
 
require(mclust)
load(".../Functions.RData")
 
# names of data series
names.spec <- list.files("../Trimmed", all.files=TRUE)
 
# load species i
i <- [1;83]
obs <- as.numeric(unlist(read.table(paste("../Trimmed/", names.spec[i], sep = ""))))
hist(obs, breaks=100, main=names.spec[i],  xlim=c(0,3))

Je pense que le problème vient d'ici :

Code :

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

1
2
3
 
# load species i
i <- [1;83]

Qu'en pensez-vous ? J'ai essayé d'autres types d'écritures, mais aucune n'a marché.

Merci pour votre aide.

[A. D.]

Bonjour,

Quelques pistes pour vous décoincer :

• Pour stocker tous les entiers de 1 à 83 dans une variable "x" par exemple :

  Code :

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

  x<-1:83

• Pour accéder à tous les éléments (83) contenus dans votre "names.spec", une des solutions est de passer par une boucle :

  Code :

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

  for(i in 1:83) { ... names.spec[i] ... }

• A chaque itération, par défaut, l'histogramme précédent sera écrasé par le nouveau. Pour y remédier, plusieurs solutions, par exemple ouvrir une nouvelle fenêtre à chaque itération :

  Code :

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

  x11()

  ouvrir un nouveau fichier d'image à chaque itération :

  Code :

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

  jpeg(paste("myfile",i,"".jpg",sep="")) ... dev.off()

  ou bien séparer l'image de sortie en plusieurs zones d'affichage (après à faire intelligemment car afficher 83 histogrammes sur une seule image, ce n'est peut-être pas terrible) :

  Code :

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

  par(mfrow=c(2,3)) # zone graphique séparer en 6 "cases", réparties en 2 lignes et 3 colonnes

Je vous laisse regarder un peu tout ça et n'hésitez pas à revenir poser des questions !

Bonne continuation


Cordialement,


A.D.