1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25
| 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(2,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")
par(mfrow=c(2,2))
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")
par(mfrow=c(1,2))
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() |
Partager