Bonjour,
J'aimerais savoir comment conserver plusieurs graphiques dans un meme pdf.
pour le moment voici mon code :
Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
 
mydata <- read.table("set01_annotes_R.txt")
mydata2 <- read.table("set01_non_annotes.txt")
pdf (file=paste ("set01_stat.pdf"),paper="a4")
par(mfrow=c(5,2)) 
hist(log10(mydata$V4),xlab="log10 exons length",ylab="frequency",col="blue",main="mean size of exons in transcripts")
hist(log10(mydata2$V4),xlab="log10 exons length",ylab="frequency",col="red",main="mean size of exons in transcripts")
hist(log10(mydata$V2),xlab="log10 transcripts length",ylab="frequency",col="blue",main="length of transcripts")
hist(log10(mydata2$V2),xlab="log10 transcripts length",ylab="frequency",col="red",main="length of transcripts")
hist(mydata$V3,xlab="exons number",ylab="frequency",col="blue",main="number of exons in transcripts")
hist(mydata2$V3,xlab="exons number",ylab="frequency",col="red",main="number of exons in transcripts")
hist((mydata$V5*100)/mydata$V2,xlab="N number",ylab="frequency",col="blue",main="percent N")
hist((mydata2$V5*100)/mydata2$V2,xlab="N number",ylab="frequency",col="red",main="percent N")
hist(mydata$V2,xlab="transcripts length",ylab="frequency",col="blue",main="length of mono-exons transcripts")
hist(mydata2$V2,xlab="transcripts length",ylab="frequency",col="red",main="length of mono-exons transcripts")
dev.off()
Mais par soucis de lecture et d'esthétique, cela ne me convient pas. J'aimerais donc savoir comment peut on conserver ses graphiques sur un pdf de 3 pages. Quelqu'un saurait il comment faire? Merci d'avance.