Mercurial > repos > youngkim > ezbamqc
diff ezBAMQC/test-data/output/data/smp1.mapq_profile.r @ 0:dfa3745e5fd8
Uploaded
author | youngkim |
---|---|
date | Thu, 24 Mar 2016 17:12:52 -0400 |
parents | |
children |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/ezBAMQC/test-data/output/data/smp1.mapq_profile.r Thu Mar 24 17:12:52 2016 -0400 @@ -0,0 +1,12 @@ +png("/sonas-hs/bsr/hpc/data/yjin/test_BAMqc/exp/test1/figs/smp1.mapq_profile.png",width=500,height=500,units="px") +mapq_val=c(0,1,3,255) +mapq_count=c(589,3042,8955,66672) +xname=c("<3","<10","<20","<30","30-255") +freq = rep(0,5) +freq[1] = sum(mapq_count[which(mapq_val<3)])/79258*100 +freq[2] = sum(mapq_count[which(mapq_val<10)])/79258*100 +freq[3] = sum(mapq_count[which(mapq_val<20)])/79258*100 +freq[4] = sum(mapq_count[which(mapq_val<30)])/79258*100 +freq[5] = 100 +barplot(freq,beside=T,xlab="Mapping Quality",border="NA",space=1.5,main="Mapping Quality",ylim=c(0,100),ylab="Cumulative proportion (%)",col="blue",names.arg=xname) +dev.state=dev.off()