Mercurial > repos > nilesh > rseqc
view test-data/output.saturation.r @ 49:6b33e31bda10 draft
Uploaded tar based on https://github.com/lparsons/galaxy_tools/tree/master/tools/rseqc 1a3c419bc0ded7c40cb2bc3e7c87bfb01ddfeba2
author | lparsons |
---|---|
date | Thu, 16 Jul 2015 17:43:43 -0400 |
parents | |
children |
line wrap: on
line source
pdf('output.saturation.pdf') par(mfrow=c(2,2)) name=c(5,10,15,20,25,30,35,40,45,50,55,60,65,70,75,80,85,90,95) S5=c() S10=c() S15=c() S20=c() S25=c() S30=c() S35=c() S40=c() S45=c() S50=c() S55=c() S60=c() S65=c() S70=c() S75=c() S80=c() S85=c() S90=c() S95=c() boxplot(100*S5,100*S10,100*S15,100*S20,100*S25,100*S30,100*S35,100*S40,100*S45,100*S50,100*S55,100*S60,100*S65,100*S70,100*S75,100*S80,100*S85,100*S90,100*S95,names=name,outline=F,ylab='Percent Relative Error',main='Q1',xlab='Resampling percentage') name=c(5,10,15,20,25,30,35,40,45,50,55,60,65,70,75,80,85,90,95) S5=c(1.0) S10=c(1.0) S15=c(1.0) S20=c(1.0) S25=c(1.0) S30=c(1.0) S35=c(1.0) S40=c(1.0) S45=c(0.259259259259) S50=c(0.333333333334) S55=c(0.39393939394) S60=c(0.444444444444) S65=c(0.0256410256403) S70=c(0.0476190476199) S75=c(0.111111111112) S80=c(0.166666666666) S85=c(0.21568627451) S90=c(0.111111111112) S95=c(0.0526315789469) boxplot(100*S5,100*S10,100*S15,100*S20,100*S25,100*S30,100*S35,100*S40,100*S45,100*S50,100*S55,100*S60,100*S65,100*S70,100*S75,100*S80,100*S85,100*S90,100*S95,names=name,outline=F,ylab='Percent Relative Error',main='Q2',xlab='Resampling percentage') name=c(5,10,15,20,25,30,35,40,45,50,55,60,65,70,75,80,85,90,95) S5=c() S10=c() S15=c() S20=c() S25=c() S30=c() S35=c() S40=c() S45=c() S50=c() S55=c() S60=c() S65=c() S70=c() S75=c() S80=c() S85=c() S90=c() S95=c() boxplot(100*S5,100*S10,100*S15,100*S20,100*S25,100*S30,100*S35,100*S40,100*S45,100*S50,100*S55,100*S60,100*S65,100*S70,100*S75,100*S80,100*S85,100*S90,100*S95,names=name,outline=F,ylab='Percent Relative Error',main='Q3',xlab='Resampling percentage') name=c(5,10,15,20,25,30,35,40,45,50,55,60,65,70,75,80,85,90,95) S5=c(4.00000000001) S10=c(1.5) S15=c(0.666666666666) S20=c(1.5) S25=c(1.00000000001) S30=c(0.666666666666) S35=c(0.428571428571) S40=c(0.25) S45=c(0.111111111111) S50=c(0.0) S55=c(0.09090909091) S60=c(0.25) S65=c(0.153846153846) S70=c(0.0714285714295) S75=c(0.0) S80=c(0.0624999999991) S85=c(0.117647058824) S90=c(0.111111111111) S95=c(0.0526315789465) boxplot(100*S5,100*S10,100*S15,100*S20,100*S25,100*S30,100*S35,100*S40,100*S45,100*S50,100*S55,100*S60,100*S65,100*S70,100*S75,100*S80,100*S85,100*S90,100*S95,names=name,outline=F,ylab='Percent Relative Error',main='Q4',xlab='Resampling percentage') dev.off()