comparison test-data/archaea_relplot.scr @ 3:e84d1fc87829 draft

planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 10ccc47885ce71e602d66e157bd475f1facbd042

2:dd991ca3d8ce

# set margin
# b, l, t, r
par(mar=c(6,2.5,1,1))

# open plot (and draw threshold)
plot(c(0),c(0),type="l",                                        
     xlab="",ylab="",
     xlim=the_xlim,ylim=c(0,maxy),
     yaxp=c(0,1,2))

## title inside of plot

    rgb(0.6,0.1,0.6,0.9)
);
colors<-c(colors,colors);

if (0>0) {
  
  for (i in 1:0) {
    orientation <- signals[signal_starts[i]+signal_sizes[i]*2];
    sig_y  <- maxy-i*anno_space;
    
    for (j in 0:(signal_sizes[i]-1)) {

      sig_x <- c(signals[signal_starts[i]+j*2],signals[signal_starts[i]+j*2+1]);
          
      ## draw arrows
      if (orientation!=0) {
        the_code <- 1+(orientation+1)/2;
        arrows(sig_x[1],sig_y,sig_x[2],sig_y,lwd=4,col=colors[i],code=the_code,angle=20,length=0.15);
      } else {

  }
}

#draw inferred on-signal
hit_color <- rgb(0.1,0.6,0.1,0.9)
  
if (0!=1) {
  
  on  <- c(0,4,33,44,49);
  off <- c(3,20,43,48,51);

  if (length(on)>0) {
    for (i in 1:length(on)) {

3:e84d1fc87829

# set margin
# b, l, t, r
par(mar=c(6,2.5,1,1))

# open plot (and draw threshold)
plot(c(0),c(0),type="l",
     xlab="",ylab="",
     xlim=the_xlim,ylim=c(0,maxy),
     yaxp=c(0,1,2))

## title inside of plot

    rgb(0.6,0.1,0.6,0.9)
);
colors<-c(colors,colors);

if (0>0) {

  for (i in 1:0) {
    orientation <- signals[signal_starts[i]+signal_sizes[i]*2];
    sig_y  <- maxy-i*anno_space;

    for (j in 0:(signal_sizes[i]-1)) {

      sig_x <- c(signals[signal_starts[i]+j*2],signals[signal_starts[i]+j*2+1]);

      ## draw arrows
      if (orientation!=0) {
        the_code <- 1+(orientation+1)/2;
        arrows(sig_x[1],sig_y,sig_x[2],sig_y,lwd=4,col=colors[i],code=the_code,angle=20,length=0.15);
      } else {

  }
}

#draw inferred on-signal
hit_color <- rgb(0.1,0.6,0.1,0.9)

if (0!=1) {

  on  <- c(0,4,33,44,49);
  off <- c(3,20,43,48,51);

  if (length(on)>0) {
    for (i in 1:length(on)) {