Mercurial > repos > nilesh > rseqc

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

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
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
comparison
equal deleted inserted replaced
48:2e6190c29c54 49:6b33e31bda10
1 pdf('output.saturation.pdf')
2 par(mfrow=c(2,2))
3 name=c(5,10,15,20,25,30,35,40,45,50,55,60,65,70,75,80,85,90,95)
4 S5=c()
5 S10=c()
6 S15=c()
7 S20=c()
8 S25=c()
9 S30=c()
10 S35=c()
11 S40=c()
12 S45=c()
13 S50=c()
14 S55=c()
15 S60=c()
16 S65=c()
17 S70=c()
18 S75=c()
19 S80=c()
20 S85=c()
21 S90=c()
22 S95=c()
23 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')
24 name=c(5,10,15,20,25,30,35,40,45,50,55,60,65,70,75,80,85,90,95)
25 S5=c(1.0)
26 S10=c(1.0)
27 S15=c(1.0)
28 S20=c(1.0)
29 S25=c(1.0)
30 S30=c(1.0)
31 S35=c(1.0)
32 S40=c(1.0)
33 S45=c(0.259259259259)
34 S50=c(0.333333333334)
35 S55=c(0.39393939394)
36 S60=c(0.444444444444)
37 S65=c(0.0256410256403)
38 S70=c(0.0476190476199)
39 S75=c(0.111111111112)
40 S80=c(0.166666666666)
41 S85=c(0.21568627451)
42 S90=c(0.111111111112)
43 S95=c(0.0526315789469)
44 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')
45 name=c(5,10,15,20,25,30,35,40,45,50,55,60,65,70,75,80,85,90,95)
46 S5=c()
47 S10=c()
48 S15=c()
49 S20=c()
50 S25=c()
51 S30=c()
52 S35=c()
53 S40=c()
54 S45=c()
55 S50=c()
56 S55=c()
57 S60=c()
58 S65=c()
59 S70=c()
60 S75=c()
61 S80=c()
62 S85=c()
63 S90=c()
64 S95=c()
65 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')
66 name=c(5,10,15,20,25,30,35,40,45,50,55,60,65,70,75,80,85,90,95)
67 S5=c(4.00000000001)
68 S10=c(1.5)
69 S15=c(0.666666666666)
70 S20=c(1.5)
71 S25=c(1.00000000001)
72 S30=c(0.666666666666)
73 S35=c(0.428571428571)
74 S40=c(0.25)
75 S45=c(0.111111111111)
76 S50=c(0.0)
77 S55=c(0.09090909091)
78 S60=c(0.25)
79 S65=c(0.153846153846)
80 S70=c(0.0714285714295)
81 S75=c(0.0)
82 S80=c(0.0624999999991)
83 S85=c(0.117647058824)
84 S90=c(0.111111111111)
85 S95=c(0.0526315789465)
86 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')
87 dev.off()