Mercurial > repos > rnateam > locarna_pairwise_p
annotate test-data/archaea_relplot.scr @ 2:535f2e41d4f7 draft
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 10ccc47885ce71e602d66e157bd475f1facbd042
author | bgruening |
---|---|
date | Mon, 05 Dec 2022 14:29:23 +0000 |
parents | c5f150d2686f |
children | a1c9c960cdd5 |
rev | line source |
---|---|
0
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
1 pdf("/tmp/tmp5yxSuG/files/000/dataset_2.dat",width=12,height=4,version="1.4") |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
2 |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
3 rel <- read.table("mlocarna_results/results/result.bmreliability"); |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
4 seqrel <- rel[[2]] |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
5 strrel <- rel[[3]] |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
6 |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
7 if ("" != "") { |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
8 seq <- "" |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
9 seq <- strsplit(seq,split="") |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
10 tab <- unlist(seq)!="-" |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
11 |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
12 seqrel<-seqrel[tab] |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
13 strrel<-strrel[tab] |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
14 } |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
15 |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
16 len<-length(seqrel) |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
17 |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
18 if (0) { |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
19 seqrel<-seqrel[len:1] |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
20 strrel<-strrel[len:1] |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
21 } |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
22 |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
23 |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
24 |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
25 if (1) { |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
26 seqrel <- seqrel/1.0 |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
27 } |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
28 |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
29 |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
30 totalrel <- seqrel+strrel; |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
31 |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
32 anno_space<-0.075 |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
33 |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
34 maxy <- max(c(1,totalrel))+anno_space*(0+1); |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
35 |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
36 |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
37 firstpos <- 1 |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
38 lastpos <- 1+len-1 |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
39 |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
40 if (0) { |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
41 the_xlim <- c(lastpos,firstpos) |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
42 } else { |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
43 the_xlim <- c(firstpos,lastpos) |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
44 } |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
45 |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
46 # set margin |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
47 # b, l, t, r |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
48 par(mar=c(6,2.5,1,1)) |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
49 |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
50 # open plot (and draw threshold) |
2
535f2e41d4f7
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 10ccc47885ce71e602d66e157bd475f1facbd042
bgruening
parents:
0
diff
changeset
|
51 plot(c(0),c(0),type="l", |
0
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
52 xlab="",ylab="", |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
53 xlim=the_xlim,ylim=c(0,maxy), |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
54 yaxp=c(0,1,2)) |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
55 |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
56 ## title inside of plot |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
57 legend("topleft","",bty="n") |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
58 |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
59 |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
60 # total reliability |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
61 polygon(c(firstpos,firstpos:lastpos,lastpos),c(0,totalrel,0),col=rgb(0.8,0.8,0.9,0.5),lwd=2,border=FALSE) |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
62 lines(firstpos:lastpos,totalrel,col="blue",lwd=2) |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
63 |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
64 # plot structure reliability |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
65 polygon(c(firstpos,firstpos:lastpos,lastpos),c(0,strrel,0),col=rgb(0.3,0.3,0.5,0.8),lwd=1,border=FALSE) |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
66 |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
67 |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
68 ## draw other signals |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
69 signals<-c(); |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
70 signal_sizes<-c(); |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
71 |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
72 signal_starts <- 1:0 |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
73 |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
74 signal_starts[1]<-1; |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
75 if (0>1) { |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
76 for (i in 2:0) { |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
77 signal_starts[i]<-signal_starts[i-1]+signal_sizes[i-1]*2+1; |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
78 } |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
79 } |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
80 |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
81 colors <- c( |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
82 rgb(0.6,0.1,0.1,0.9), |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
83 rgb(0.6,0.6,0.1,0.9), |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
84 rgb(0.1,0.6,0.6,0.9), |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
85 rgb(0.6,0.1,0.6,0.9) |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
86 ); |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
87 colors<-c(colors,colors); |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
88 |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
89 if (0>0) { |
2
535f2e41d4f7
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 10ccc47885ce71e602d66e157bd475f1facbd042
bgruening
parents:
0
diff
changeset
|
90 |
0
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
91 for (i in 1:0) { |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
92 orientation <- signals[signal_starts[i]+signal_sizes[i]*2]; |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
93 sig_y <- maxy-i*anno_space; |
2
535f2e41d4f7
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 10ccc47885ce71e602d66e157bd475f1facbd042
bgruening
parents:
0
diff
changeset
|
94 |
0
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
95 for (j in 0:(signal_sizes[i]-1)) { |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
96 |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
97 sig_x <- c(signals[signal_starts[i]+j*2],signals[signal_starts[i]+j*2+1]); |
2
535f2e41d4f7
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 10ccc47885ce71e602d66e157bd475f1facbd042
bgruening
parents:
0
diff
changeset
|
98 |
0
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
99 ## draw arrows |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
100 if (orientation!=0) { |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
101 the_code <- 1+(orientation+1)/2; |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
102 arrows(sig_x[1],sig_y,sig_x[2],sig_y,lwd=4,col=colors[i],code=the_code,angle=20,length=0.15); |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
103 } else { |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
104 lines(sig_x,c(sig_y,sig_y),lwd=4,col=colors[i]); |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
105 } |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
106 } |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
107 } |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
108 } |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
109 |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
110 #draw inferred on-signal |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
111 hit_color <- rgb(0.1,0.6,0.1,0.9) |
2
535f2e41d4f7
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 10ccc47885ce71e602d66e157bd475f1facbd042
bgruening
parents:
0
diff
changeset
|
112 |
0
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
113 if (0!=1) { |
2
535f2e41d4f7
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 10ccc47885ce71e602d66e157bd475f1facbd042
bgruening
parents:
0
diff
changeset
|
114 |
0
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
115 on <- c(0,4,33,44,49); |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
116 off <- c(3,20,43,48,51); |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
117 |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
118 if (length(on)>0) { |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
119 for (i in 1:length(on)) { |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
120 lines(c(1+on[i],1+off[i]-1),c(maxy,maxy),lwd=7,col=hit_color); |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
121 } |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
122 } |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
123 |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
124 ### draw on/off values |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
125 if (0) { |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
126 lines(c(1,1+len),c(0.0355791,0.0355791),lty=2,lwd=1) |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
127 lines(c(1,1+len),c(0.702675,0.702675),lty=2,lwd=1) |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
128 } |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
129 } |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
130 |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
131 |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
132 signal_names<-c(); |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
133 |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
134 if (length(signal_names)>0 || (0!=1)) { |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
135 legend("bottom",c("LocARNA",signal_names),lwd=7,col=c(hit_color,colors),horiz=TRUE,inset=-0.4); |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
136 # ,xpd=TRUE |
c5f150d2686f
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 287021573c592fdb70fdbbc88943aa16a8740fc0
rnateam
parents:
diff
changeset
|
137 } |