Mercurial > repos > artbio > small_rna_signatures

annotate signature.r @ 0:3e0ea204d09e draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
"planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
author artbio
date Tue, 07 Jan 2020 11:59:34 +0000
parents
children ca2b04cdbf4d
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
1 ## Setup R error handling to go to stderr
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
2 #options(show.error.messages=F,
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
3 #error = function () { cat( geterrmessage(), file=stderr() ); q( "no", 1, F ) } )
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
4 warnings()
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
5
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
6 library(RColorBrewer)
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
7 library(lattice)
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
8 library(latticeExtra)
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
9 library(grid)
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
10 library(gridExtra)
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
11 library(optparse)
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
12
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
13 option_list <- list(
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
14 make_option("--h_dataframe", type="character", help="path to h-signature dataframe"),
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
15 make_option("--z_dataframe", type="character", help="path to z-signature dataframe"),
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
16 make_option("--plot_method", type = "character", help="How data should be plotted (global or lattice)"),
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
17 make_option("--pdf", type = "character", help="path to the pdf file with plots"),
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
18 make_option("--title", type = "character", help="Graph Title")
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
19 )
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
20
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
21 parser <- OptionParser(usage = "%prog [options] file", option_list = option_list)
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
22 args = parse_args(parser)
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
23
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
24
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
25
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
26
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
27 # data frames implementation
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
28 h_dataframe = read.delim(args$h_dataframe, header=F)
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
29 colnames(h_dataframe) = c("chrom", "overlap", "sig", "z-score")
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
30 h_dataframe$sig = h_dataframe$sig * 100 # to get probs in %
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
31 z_dataframe = read.delim(args$z_dataframe, header=F)
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
32 colnames(z_dataframe) = c("chrom", "overlap", "sig", "z-score")
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
33
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
34 # functions
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
35 globalgraph = function () {
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
36 pdf(args$pdf)
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
37 par(mfrow=c(2,2),oma = c(0, 0, 3, 0))
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
38
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
39 plot(z_dataframe[z_dataframe$chrom == "all_chromosomes", c(2,3)],
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
40 type = "h", main="Numbers of pairs", cex.main=1, xlab="overlap (nt)",
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
41 ylab="Numbers of pairs", col="darkslateblue", lwd=4)
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
42
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
43 plot(z_dataframe[z_dataframe$chrom == "all_chromosomes", c(2,4)],
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
44 type = "l", main="Number of pairs Z-scores", cex.main=1, xlab="overlap (nt)",
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
45 ylab="z-score", pch=19, cex=0.2, col="darkslateblue", lwd=2)
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
46
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
47 plot(h_dataframe[h_dataframe$chrom == "all_chromosomes", c(2,3)],
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
48 type = "l", main="Overlap probabilities", cex.main=1, xlab="overlap (nt)",
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
49 ylab="Probability [%]", ylim=c(0,50), pch=19, col="darkslateblue", lwd=2)
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
50
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
51 plot(h_dataframe[h_dataframe$chrom == "all_chromosomes", c(2,4)],
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
52 type = "l", main="Overlap Probability Z-scores", cex.main=1,
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
53 xlab="overlap (nt)", ylab="z-score", pch=19, cex=0.2,
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
54 col="darkslateblue", lwd=2)
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
55
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
56 mtext(args$title, outer = TRUE, cex=1)
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
57 dev.off()
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
58 }
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
59
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
60 treillisgraph = function (df, ...) {
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
61 pdf(args$pdf, paper="special", height=11.69, width=6 ) # 8.2677
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
62 p = xyplot(sig ~ overlap|factor(method, levels=unique(method))+chrom, data = df,
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
63 type = "l",
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
64 col='darkblue',
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
65 cex=0.5,
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
66 scales=list(y=list(tick.number=4, relation="free", cex=0.6, rot=0), x=list(cex=0.6, alternating=FALSE)),
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
67 xlab = "Overlap",
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
68 ylab = "signature (Nbr of pairs / Overlap prob.)",
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
69 main = args$title,
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
70 par.strip.text=list(cex=.5),
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
71 pch=19, lwd =2,
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
72 as.table=TRUE,
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
73 layout=c(2,12),
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
74 newpage = T,
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
75 ...)
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
76 plot(p)
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
77 dev.off()
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
78 }
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
79
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
80 # main
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
81
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
82 if (args$plot_method=="global") {
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
83 globalgraph()
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
84 }
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
85
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
86 if(args$plot_method=="lattice") {
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
87 # rearrange dataframes
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
88 h_sig = h_dataframe[,c(1,2,3)]
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
89 h_sig = cbind(rep("Overlap Prob (%)", length(h_sig[,1])), h_sig)
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
90 colnames(h_sig) = c("method", "chrom", "overlap", "sig")
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
91 z_pairs = z_dataframe[,c(1,2,3)]
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
92 z_pairs = cbind(rep("Nbr of pairs", length(z_pairs[,1])), z_pairs)
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
93 colnames(z_pairs) = c("method", "chrom", "overlap", "sig")
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
94 lattice_df = rbind(z_pairs, h_sig)
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
95 par.settings.treillis=list(strip.background = list(
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
96 col = c("lightblue", "lightgreen")))
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
97
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
98 treillisgraph(lattice_df, par.settings=par.settings.treillis)
3e0ea204d09e "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit 4e84386b619a7941f23d175d7fc86aba7990ac36"
artbio
parents:
diff changeset
99 }