Mercurial > repos > pjbriggs > ceas
annotate test-data/ceas_out2.log @ 11:d3bd9caee1b7 draft default tip
Uploaded version 1.0.2-4.
author | pjbriggs |
---|---|
date | Wed, 26 Sep 2018 06:37:20 -0400 |
parents | 4e2883bb058d |
children |
rev | line source |
---|---|
4
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
1 ceas -- 0.9.9.7 (package version 1.0.2) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
2 INFO @ Tue, 23 Jun 2015 11:24:09: |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
3 # ARGUMENTS: |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
4 # name = ceas |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
5 # gene annotation table = galGal3.refGene |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
6 # BED file = ceas_in.bed |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
7 # WIG file = ceas_in_stp1000.wig |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
8 # extra BED file = None |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
9 # ChIP annotation = On |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
10 # gene-centered annotation = On |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
11 # average profiling = On |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
12 # dump profiles = Off |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
13 # re-annotation for genome background (ChIP region annotation) = False |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
14 # promoter sizes (ChIP region annotation) = 1000,2000,3000 bp |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
15 # downstream sizes (ChIP region annotation) = 1000,2000,3000 bp |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
16 # bidrectional promoter sizes (ChIP region annotation) = 2500,5000 bp |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
17 # span size (gene-centered annotation) = 3000 bp |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
18 # profiling resolution (average profiling) = 50 bp |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
19 # relative distance wrt TSS and TTS (average profiling) = 3000 bp |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
20 INFO @ Tue, 23 Jun 2015 11:24:09: #1 read the gene table... |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
21 INFO @ Tue, 23 Jun 2015 11:24:09: #2 read the bed file of ChIP regions... |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
22 INFO @ Tue, 23 Jun 2015 11:24:09: #3 perform gene-centered annotation... |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
23 INFO @ Tue, 23 Jun 2015 11:24:09: #4 See ceas.xls for gene-centered annotation! |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
24 INFO @ Tue, 23 Jun 2015 11:24:09: #5 read the pre-computed genome bg annotation... |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
25 INFO @ Tue, 23 Jun 2015 11:24:09: #6 perform ChIP region annotation... |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
26 INFO @ Tue, 23 Jun 2015 11:24:09: #7 write a R script of ChIP region annotation... |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
27 INFO @ Tue, 23 Jun 2015 11:24:09: #8-1 run wig profiling of chr26... |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
28 INFO @ Tue, 23 Jun 2015 11:24:09: #9 append an R script of wig profiling... |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
29 |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
30 R version 3.1.2 (2014-10-31) -- "Pumpkin Helmet" |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
31 Copyright (C) 2014 The R Foundation for Statistical Computing |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
32 Platform: x86_64-redhat-linux-gnu (64-bit) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
33 |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
34 R is free software and comes with ABSOLUTELY NO WARRANTY. |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
35 You are welcome to redistribute it under certain conditions. |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
36 Type 'license()' or 'licence()' for distribution details. |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
37 |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
38 Natural language support but running in an English locale |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
39 |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
40 R is a collaborative project with many contributors. |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
41 Type 'contributors()' for more information and |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
42 'citation()' on how to cite R or R packages in publications. |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
43 |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
44 Type 'demo()' for some demos, 'help()' for on-line help, or |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
45 'help.start()' for an HTML browser interface to help. |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
46 Type 'q()' to quit R. |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
47 |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
48 > # ARGUMENTS: |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
49 > # name = ceas |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
50 > # gene annotation table = galGal3.refGene |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
51 > # BED file = ceas_in.bed |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
52 > # WIG file = ceas_in_stp1000.wig |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
53 > # extra BED file = None |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
54 > # ChIP annotation = On |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
55 > # gene-centered annotation = On |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
56 > # average profiling = On |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
57 > # dump profiles = Off |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
58 > # re-annotation for genome background (ChIP region annotation) = False |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
59 > # promoter sizes (ChIP region annotation) = 1000,2000,3000 bp |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
60 > # downstream sizes (ChIP region annotation) = 1000,2000,3000 bp |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
61 > # bidrectional promoter sizes (ChIP region annotation) = 2500,5000 bp |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
62 > # span size (gene-centered annotation) = 3000 bp |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
63 > # profiling resolution (average profiling) = 50 bp |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
64 > # relative distance wrt TSS and TTS (average profiling) = 3000 bp |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
65 > pdf("ceas.pdf",height=11.5,width=8.5) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
66 > |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
67 > # 11:24:09 Tue, 23 Jun 2015 |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
68 > # |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
69 > # ChIP annotation |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
70 > # |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
71 > |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
72 > |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
73 > # |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
74 > # Chromosomal Distribution |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
75 > # |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
76 > |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
77 > par(mar=c(4, 4, 5, 3.8),oma=c(4, 2, 4, 2)) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
78 > r0<-c(100.0) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
79 > r1<-c(100.0) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
80 > height<-rbind(r0,r1) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
81 > names=c("26") |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
82 > mp<-barplot(height=height,names=names,beside=TRUE,horiz=TRUE,col=c("#5FA1C1","#EB9D86"),main="Chromosomal Distribution of ChIP Regions",xlab="Percentage %",ylab="Chromosome",border=FALSE,xlim=c(0.000000,183.333333),cex.names=1) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
83 > text(x=c(100.0),y=mp[1,],label=c("100.0 %"),pos=4,offset=0.2,cex=0.9) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
84 > text(x=c(100.0),y=mp[2,],label=c("100.0 % (<=4.9e-324)"),pos=4,offset=0.2,cex=0.9) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
85 > legend("right",legend=c("Genome","ChIP (p-value)"),col=c("#5FA1C1","#EB9D86"),pch=15,bty="n") |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
86 > |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
87 > # |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
88 > # Promoter,Bipromoter,Downstream, Gene and Regions of interest |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
89 > # |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
90 > |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
91 > par(mfrow=c(4, 1),mar=c(4, 4, 5, 3.8),oma=c(4, 2, 4, 2)) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
92 > r0<-c(1.8532425688606797, 3.616851183410451, 5.322318854623416) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
93 > r1<-c(0.0, 0.0, 0.0) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
94 > height<-rbind(r0,r1) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
95 > names=c("<=1000 bp","<=2000 bp","<=3000 bp") |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
96 > mp<-barplot(height=height,names=names,beside=TRUE,horiz=FALSE,col=c("#5FA1C1","#EB9D86"),main="Promoter",ylab="Percentage %",border=FALSE,ylim=c(0.000000,9.757585),cex.names=1) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
97 > text(x=mp[1,],y=c(1.8532425688606797, 3.616851183410451, 5.322318854623416),label=c("1.9 %","3.6 %","5.3 %"),pos=3,offset=0.2) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
98 > text(x=mp[2,],y=c(0.0, 0.0, 0.0),label=c("0.000 % |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
99 + (0.981)","0.000 % |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
100 + (0.964)","0.000 % |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
101 + (0.947)"),pos=3,offset=0.2) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
102 > legend("topleft",legend=c("Genome","ChIP (p-value)"),col=c("#5FA1C1","#EB9D86"),pch=15,bty="n") |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
103 > r0<-c(0.03876062889120376, 0.03876062889120376) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
104 > r1<-c(0.0, 0.0) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
105 > height<-rbind(r0,r1) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
106 > names=c("<=2500 bp","<=5000 bp") |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
107 > mp<-barplot(height=height,names=names,beside=TRUE,horiz=FALSE,col=c("#5FA1C1","#EB9D86"),main="Bidirectional Promoter",ylab="Percentage %",border=FALSE,ylim=c(0.000000,0.071061),cex.names=1) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
108 > text(x=mp[1,],y=c(0.03876062889120376, 0.03876062889120376),label=c("0.04 %","0.04 %"),pos=3,offset=0.2) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
109 > text(x=mp[2,],y=c(0.0, 0.0),label=c("0.000 % |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
110 + (1.000)","0.000 % |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
111 + (1.000)"),pos=3,offset=0.2) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
112 > legend("topleft",legend=c("Genome","ChIP (p-value)"),col=c("#5FA1C1","#EB9D86"),pch=15,bty="n") |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
113 > r0<-c(1.8290171758036773, 3.4690762857627364, 4.980740812519683) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
114 > r1<-c(0.0, 0.0, 0.0) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
115 > height<-rbind(r0,r1) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
116 > names=c("<=1000 bp","<=2000 bp","<=3000 bp") |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
117 > mp<-barplot(height=height,names=names,beside=TRUE,horiz=FALSE,col=c("#5FA1C1","#EB9D86"),main="Downstream",ylab="Percentage %",border=FALSE,ylim=c(0.000000,9.131358),cex.names=1) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
118 > text(x=mp[1,],y=c(1.8290171758036773, 3.4690762857627364, 4.980740812519683),label=c("1.8 %","3.5 %","5.0 %"),pos=3,offset=0.2) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
119 > text(x=mp[2,],y=c(0.0, 0.0, 0.0),label=c("0.000 % |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
120 + (0.982)","0.000 % |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
121 + (0.965)","0.000 % |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
122 + (0.950)"),pos=3,offset=0.2) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
123 > legend("topleft",legend=c("Genome","ChIP (p-value)"),col=c("#5FA1C1","#EB9D86"),pch=15,bty="n") |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
124 > r0<-c(0.2034933016788197, 1.3978051793890356, 2.359553283752029, 19.734005184234114, 23.694856949054) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
125 > r1<-c(0.0, 0.0, 0.0, 0.0, 0.0) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
126 > height<-rbind(r0,r1) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
127 > names=c("5'UTR","3'UTR","Coding Exon","Intron","All") |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
128 > mp<-barplot(height=height,names=names,beside=TRUE,horiz=FALSE,col=c("#5FA1C1","#EB9D86"),main="Gene",ylab="Percentage %",border=FALSE,ylim=c(0.000000,43.440571),cex.names=1) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
129 > text(x=mp[1,],y=c(0.2034933016788197, 1.3978051793890356, 2.359553283752029, 19.734005184234114, 23.694856949054),label=c("0.2 %","1.4 %","2.4 %","19.7 %","23.7 %"),pos=3,offset=0.2) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
130 > text(x=mp[2,],y=c(0.0, 0.0, 0.0, 0.0, 0.0),label=c("0.000 % |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
131 + (0.998)","0.000 % |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
132 + (0.986)","0.000 % |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
133 + (0.976)","0.000 % |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
134 + (0.803)","0.000 % |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
135 + (0.763)"),pos=3,offset=0.2) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
136 > legend("topleft",legend=c("Genome","ChIP (p-value)"),col=c("#5FA1C1","#EB9D86"),pch=15,bty="n") |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
137 > |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
138 > # |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
139 > # Distribution of Genome and ChIP regions over cis-regulatory element |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
140 > # Note that the x may be modified for better graphics in case a value is too small |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
141 > # Thus, look at the labels of the pie chart to get the real percentage values |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
142 > # |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
143 > |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
144 > par(mfcol=c(2, 2),mar=c(3, 3, 4, 2.8),oma=c(4, 2, 4, 2)) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
145 > x<-c(0.018532,0.017055,0.016037,0.017830,0.015092,0.014051,0.010000,0.013833,0.023014,0.192592,0.670292) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
146 > pie(x=x,labels=c("1.9 %","1.7 %","1.6 %","1.8 %","1.5 %","1.4 %","0.2 %","1.4 %","2.3 %","19.3 %","67.0 %"),main="Genome",col=c("#445FA2","#EB9D86","#799F7A","#6C527F","#5FA1C1","#E8BB77","#A8C5EF","#FDCDB9","#C6E6B5","#F1D5EE","#B4E1F6"),clockwise=TRUE,border=FALSE,radius=0.9,cex=0.8,init.angle=90,density=100) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
147 > x<-c(0.000000,1.000000) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
148 > y<-c(0.000000,1.000000) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
149 > plot(x, y,type="n",main="",xlab="",ylab="",frame=FALSE,axes=FALSE,xaxt="s",yaxt="s") |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
150 > legend("top",legend=c("Promoter (<=1000 bp): 1.9 %","Promoter (1000-2000 bp): 1.7 %","Promoter (2000-3000 bp): 1.6 %","Downstream (<=1000 bp): 1.8 %","Downstream (1000-2000 bp): 1.5 %","Downstream (2000-3000 bp): 1.4 %","5'UTR: 0.2 %","3'UTR: 1.4 %","Coding exon: 2.3 %","Intron: 19.3 %","Distal intergenic: 67.0 %"),col=c("#445FA2","#EB9D86","#799F7A","#6C527F","#5FA1C1","#E8BB77","#A8C5EF","#FDCDB9","#C6E6B5","#F1D5EE","#B4E1F6"),pch=15,bty="n") |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
151 > x<-c(0.010000,0.010000,0.010000,0.010000,0.010000,0.010000,0.010000,0.010000,0.010000,0.010000,1.000000) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
152 > pie(x=x,labels=c("0.000 %","0.000 %","0.000 %","0.000 %","0.000 %","0.000 %","0.000 %","0.000 %","0.000 %","0.000 %","100.0 %"),main="ChIP",col=c("#445FA2","#EB9D86","#799F7A","#6C527F","#5FA1C1","#E8BB77","#A8C5EF","#FDCDB9","#C6E6B5","#F1D5EE","#B4E1F6"),clockwise=TRUE,border=FALSE,radius=0.9,cex=0.8,init.angle=90,density=100) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
153 > x<-c(0.000000,1.000000) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
154 > y<-c(0.000000,1.000000) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
155 > plot(x, y,type="n",main="",xlab="",ylab="",frame=FALSE,axes=FALSE,xaxt="s",yaxt="s") |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
156 > legend("top",legend=c("Promoter (<=1000 bp): 0.000 %","Promoter (1000-2000 bp): 0.000 %","Promoter (2000-3000 bp): 0.000 %","Downstream (<=1000 bp): 0.000 %","Downstream (1000-2000 bp): 0.000 %","Downstream (2000-3000 bp): 0.000 %","5'UTR: 0.000 %","3'UTR: 0.000 %","Coding exon: 0.000 %","Intron: 0.000 %","Distal intergenic: 100.0 %"),col=c("#445FA2","#EB9D86","#799F7A","#6C527F","#5FA1C1","#E8BB77","#A8C5EF","#FDCDB9","#C6E6B5","#F1D5EE","#B4E1F6"),pch=15,bty="n") |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
157 > |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
158 > # |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
159 > # ChIP regions over the genome |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
160 > # |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
161 > |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
162 > par(mar=c(4, 4, 5, 3.8),oma=c(4, 2, 4, 2)) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
163 > layout(matrix(c(1, 0, 2, 2), 2, 2, byrow = TRUE),widths=c(1, 1),heights=c(1, 5)) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
164 > x<-c(0.000000,0.000000) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
165 > y<-c(0.000000,1.000000) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
166 > plot(x, y,type="n",main="Distribution of Peak Heights",xlab="",ylab="",xlim=c(0.000000,0.000000),ylim=c(0.000000,1.000000),frame=FALSE,xaxt="s",yaxt="n",cex=0.9) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
167 > x<-c(0.000000,0.000000,0.000000,0.000000) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
168 > y<-c(0.000000,0.000000,1.000000,1.000000) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
169 > polygon(x,y,col=c("black")) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
170 > x <- c(0.000000) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
171 > y<-c(0.800000) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
172 > lines(x, y,xlim=c(0, 0.0),ylim=c(0, 1),type="l",col=c("cyan"),lwd=2) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
173 > x<-c(0.000000,4127518.000000) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
174 > y<-c(0.855556,1.144444) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
175 > plot(x, y,type="n",main="ChIP Regions (Peaks) over Chromosomes",xlab="Chromosome Size (bp)",ylab="Chromosome",xlim=c(0.000000,4127518.000000),ylim=c(0.855556,1.144444),frame=FALSE,xaxt="s",yaxt="n") |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
176 > start <- c(4119129) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
177 > end <- c(4119130) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
178 > vals <- c(0.0) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
179 > vals[vals > 0.0] <- 0.0 |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
180 > vals[vals < 0] <- 0 |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
181 > heights <- 0.288889 * ((vals - 0)/(0.0 - 0)) + 0.855555555556 |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
182 > for (i in 1:length(heights)) { |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
183 + polygon(x=c(start[i], end[i], end[i], start[i]), y=c(0.855555555556, 0.855555555556, heights[i], heights[i]), col=c("#CC0000"), border=c("#CC0000")) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
184 + } |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
185 > mtext("26",side=2,line=0,outer=FALSE,at=1.0) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
186 > par(mar=c(4, 4, 5, 3.8),oma=c(4, 2, 4, 2)) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
187 > layout(matrix(c(1, 2, 3, 3, 4, 5), 3, 2, byrow = TRUE),widths=c(1, 1),heights=c(1, 1, 1)) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
188 > x<-c(-3000.000000,-2950.000000,-2900.000000,-2850.000000,-2800.000000,-2750.000000,-2700.000000,-2650.000000,-2600.000000,-2550.000000,-2500.000000,-2450.000000,-2400.000000,-2350.000000,-2300.000000,-2250.000000,-2200.000000,-2150.000000,-2100.000000,-2050.000000,-2000.000000,-1950.000000,-1900.000000,-1850.000000,-1800.000000,-1750.000000,-1700.000000,-1650.000000,-1600.000000,-1550.000000,-1500.000000,-1450.000000,-1400.000000,-1350.000000,-1300.000000,-1250.000000,-1200.000000,-1150.000000,-1100.000000,-1050.000000,-1000.000000,-950.000000,-900.000000,-850.000000,-800.000000,-750.000000,-700.000000,-650.000000,-600.000000,-550.000000,-500.000000,-450.000000,-400.000000,-350.000000,-300.000000,-250.000000,-200.000000,-150.000000,-100.000000,-50.000000,0.000000,50.000000,100.000000,150.000000,200.000000,250.000000,300.000000,350.000000,400.000000,450.000000,500.000000,550.000000,600.000000,650.000000,700.000000,750.000000,800.000000,850.000000,900.000000,950.000000,1000.000000,1050.000000,1100.000000,1150.000000,1200.000000,1250.000000,1300.000000,1350.000000,1400.000000,1450.000000,1500.000000,1550.000000,1600.000000,1650.000000,1700.000000,1750.000000,1800.000000,1850.000000,1900.000000,1950.000000,2000.000000,2050.000000,2100.000000,2150.000000,2200.000000,2250.000000,2300.000000,2350.000000,2400.000000,2450.000000,2500.000000,2550.000000,2600.000000,2650.000000,2700.000000,2750.000000,2800.000000,2850.000000,2900.000000,2950.000000,3000.000000) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
189 > y<-c(0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,217.391304,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
190 > plot(x, y,type="l",main="Average Profile near TSS",xlab="Relative Distance to TSS (bp)",ylab="Average Profile",col=c("#C8524D"),xaxt="s",yaxt="s",lwd=2) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
191 > abline(v=0.000000,lty=2,col=c("black")) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
192 > x<-c(-3000.000000,-2950.000000,-2900.000000,-2850.000000,-2800.000000,-2750.000000,-2700.000000,-2650.000000,-2600.000000,-2550.000000,-2500.000000,-2450.000000,-2400.000000,-2350.000000,-2300.000000,-2250.000000,-2200.000000,-2150.000000,-2100.000000,-2050.000000,-2000.000000,-1950.000000,-1900.000000,-1850.000000,-1800.000000,-1750.000000,-1700.000000,-1650.000000,-1600.000000,-1550.000000,-1500.000000,-1450.000000,-1400.000000,-1350.000000,-1300.000000,-1250.000000,-1200.000000,-1150.000000,-1100.000000,-1050.000000,-1000.000000,-950.000000,-900.000000,-850.000000,-800.000000,-750.000000,-700.000000,-650.000000,-600.000000,-550.000000,-500.000000,-450.000000,-400.000000,-350.000000,-300.000000,-250.000000,-200.000000,-150.000000,-100.000000,-50.000000,0.000000,50.000000,100.000000,150.000000,200.000000,250.000000,300.000000,350.000000,400.000000,450.000000,500.000000,550.000000,600.000000,650.000000,700.000000,750.000000,800.000000,850.000000,900.000000,950.000000,1000.000000,1050.000000,1100.000000,1150.000000,1200.000000,1250.000000,1300.000000,1350.000000,1400.000000,1450.000000,1500.000000,1550.000000,1600.000000,1650.000000,1700.000000,1750.000000,1800.000000,1850.000000,1900.000000,1950.000000,2000.000000,2050.000000,2100.000000,2150.000000,2200.000000,2250.000000,2300.000000,2350.000000,2400.000000,2450.000000,2500.000000,2550.000000,2600.000000,2650.000000,2700.000000,2750.000000,2800.000000,2850.000000,2900.000000,2950.000000,3000.000000) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
193 > y<-c(0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,217.391304,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,217.391304,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,217.391304,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,217.391304,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,217.391304,0.000000,0.000000,0.000000,0.000000,217.391304,0.000000,0.000000,0.000000,0.000000,0.000000,217.391304,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,217.391304,0.000000,0.000000) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
194 > plot(x, y,type="l",main="Average Profile near TTS",xlab="Relative Distance to TTS (bp)",ylab="Average Profile",col=c("#C8524D"),xaxt="s",yaxt="s",lwd=2) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
195 > abline(v=0.000000,lty=2,col=c("black")) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
196 > x<-c(-1000.000000,-950.000000,-900.000000,-850.000000,-800.000000,-750.000000,-700.000000,-650.000000,-600.000000,-550.000000,-500.000000,-450.000000,-400.000000,-350.000000,-300.000000,-250.000000,-200.000000,-150.000000,-100.000000,-50.000000,0.000000,50.000000,100.000000,150.000000,200.000000,250.000000,300.000000,350.000000,400.000000,450.000000,500.000000,550.000000,600.000000,650.000000,700.000000,750.000000,800.000000,850.000000,900.000000,950.000000,1000.000000,1050.000000,1100.000000,1150.000000,1200.000000,1250.000000,1300.000000,1350.000000,1400.000000,1450.000000,1500.000000,1550.000000,1600.000000,1650.000000,1700.000000,1750.000000,1800.000000,1850.000000,1900.000000,1950.000000,2000.000000,2050.000000,2100.000000,2150.000000,2200.000000,2250.000000,2300.000000,2350.000000,2400.000000,2450.000000,2500.000000,2550.000000,2600.000000,2650.000000,2700.000000,2750.000000,2800.000000,2850.000000,2900.000000,2950.000000,3000.000000,3050.000000,3100.000000,3150.000000,3200.000000,3250.000000,3300.000000,3350.000000,3400.000000,3450.000000,3500.000000,3550.000000,3600.000000,3650.000000,3700.000000,3750.000000,3800.000000,3850.000000,3900.000000,3950.000000,4000.000000) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
197 > y<-c(0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,238.095238,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,217.391304,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,217.391304,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
198 > plot(x, y,type="l",main="Average Gene Profile",xlab="Upstream (bp), 3000 bp of Meta-gene, Downstream (bp)",ylab="Average Profile",col=c("#C8524D"),xaxt="s",yaxt="s",lwd=2) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
199 > abline(v=0.000000,lty=2,col=c("black")) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
200 > abline(v=3000.000000,lty=2,col=c("black")) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
201 > x<-c(0.000000,3.333333,6.666667,10.000000,13.333333,16.666667,20.000000,23.333333,26.666667,30.000000,33.333333,36.666667,40.000000,43.333333,46.666667,50.000000,53.333333,56.666667,60.000000,63.333333,66.666667,70.000000,73.333333,76.666667,80.000000,83.333333,86.666667,90.000000,93.333333,96.666667,100.000000) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
202 > y<-c(0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,250.000000,0.000000,0.000000,0.000000,0.000000) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
203 > plot(x, y,type="l",main="Average Concatenated Exon Profile",xlab="Relative Location (%)",ylab="Average Profile",col=c("#C8524D"),ylim=c(0.000000,250.000000),xaxt="s",yaxt="s",lwd=2) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
204 > x<-c(0.000000,3.333333,6.666667,10.000000,13.333333,16.666667,20.000000,23.333333,26.666667,30.000000,33.333333,36.666667,40.000000,43.333333,46.666667,50.000000,53.333333,56.666667,60.000000,63.333333,66.666667,70.000000,73.333333,76.666667,80.000000,83.333333,86.666667,90.000000,93.333333,96.666667,100.000000) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
205 > y<-c(0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
206 > plot(x, y,type="l",main="Average Concatenated Intron Profile",xlab="Relative Location (%)",ylab="Average Profile",col=c("#C8524D"),ylim=c(0.000000,250.000000),xaxt="s",yaxt="s",lwd=2) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
207 > par(mfrow=c(3, 2),mar=c(4, 4, 5, 3.8),oma=c(4, 2, 4, 2)) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
208 > x<-c(0.000000,50.000000,100.000000) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
209 > y<-c(0.000000,0.000000,0.000000) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
210 > plot(x, y,type="l",main="Average Exon Profile |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
211 + (56 <= length < 109 bp)",xlab="Relative Location (%)",ylab="Average Profile",col=c("#C8524D"),ylim=c(0.000000,132.596685),xaxt="s",yaxt="s",lwd=2) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
212 > x<-c(0.000000,25.000000,50.000000,75.000000,100.000000) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
213 > y<-c(0.000000,0.000000,0.000000,0.000000,0.000000) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
214 > plot(x, y,type="l",main="Average Intron Profile |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
215 + (110 <= length < 345 bp)",xlab="Relative Location (%)",ylab="Average Profile",col=c("#C8524D"),ylim=c(0.000000,132.596685),xaxt="s",yaxt="s",lwd=2) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
216 > x<-c(0.000000,50.000000,100.000000) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
217 > y<-c(0.000000,0.000000,0.000000) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
218 > plot(x, y,type="l",main="Average Exon Profile |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
219 + (109 <= length < 160 bp)",xlab="Relative Location (%)",ylab="Average Profile",col=c("#C8524D"),ylim=c(0.000000,132.596685),xaxt="s",yaxt="s",lwd=2) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
220 > x<-c(0.000000,11.111111,22.222222,33.333333,44.444444,55.555556,66.666667,77.777778,88.888889,100.000000) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
221 > y<-c(0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
222 > plot(x, y,type="l",main="Average Intron Profile |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
223 + (344 <= length < 686 bp)",xlab="Relative Location (%)",ylab="Average Profile",col=c("#C8524D"),ylim=c(0.000000,132.596685),xaxt="s",yaxt="s",lwd=2) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
224 > x<-c(0.000000,33.333333,66.666667,100.000000) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
225 > y<-c(0.000000,0.000000,0.000000,0.000000) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
226 > plot(x, y,type="l",main="Average Exon Profile |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
227 + (160 <= length < 375 bp)",xlab="Relative Location (%)",ylab="Average Profile",col=c("#C8524D"),ylim=c(0.000000,132.596685),xaxt="s",yaxt="s",lwd=2) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
228 > x<-c(0.000000,4.761905,9.523810,14.285714,19.047619,23.809524,28.571429,33.333333,38.095238,42.857143,47.619048,52.380952,57.142857,61.904762,66.666667,71.428571,76.190476,80.952381,85.714286,90.476190,95.238095,100.000000) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
229 > y<-c(0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,110.497238,0.000000) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
230 > plot(x, y,type="l",main="Average Intron Profile |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
231 + (685 <= length < 2653 bp)",xlab="Relative Location (%)",ylab="Average Profile",col=c("#C8524D"),ylim=c(0.000000,132.596685),xaxt="s",yaxt="s",lwd=2) |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
232 > dev.off() |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
233 null device |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
234 1 |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
235 > |
4e2883bb058d
Synchronise version with one on main toolshed (version of 30/06/2015)
pjbriggs
parents:
diff
changeset
|
236 INFO @ Tue, 23 Jun 2015 11:24:09: #... cong! See ceas.pdf for the graphical results of CEAS! |