view test-data/output.saturation.r @ 49:6b33e31bda10 draft

Uploaded tar based on https://github.com/lparsons/galaxy_tools/tree/master/tools/rseqc 1a3c419bc0ded7c40cb2bc3e7c87bfb01ddfeba2

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()