annotate test-data/ceas_out1.log @ 4:cd6a48ffd093 draft default tip

Uploaded version 1.0.2-4.
author pjbriggs
date Wed, 26 Sep 2018 08:03:10 -0400
parents f411ce97a351
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
1 ceas -- 0.9.9.7 (package version 1.0.2)
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
2 INFO @ Tue, 23 Jun 2015 09:12:22:
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
3 # ARGUMENTS:
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
4 # name = ceas
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
5 # gene annotation table = galGal3.refGene
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
6 # BED file = ceas_in.bed
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
7 # WIG file = None
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
8 # extra BED file = None
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
9 # ChIP annotation = On
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
10 # gene-centered annotation = On
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
11 # average profiling = Off
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
12 # dump profiles = Off
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
13 # re-annotation for genome background (ChIP region annotation) = False
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
14 # promoter sizes (ChIP region annotation) = 1000,2000,3000 bp
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
15 # downstream sizes (ChIP region annotation) = 1000,2000,3000 bp
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
16 # bidrectional promoter sizes (ChIP region annotation) = 2500,5000 bp
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
17 # span size (gene-centered annotation) = 3000 bp
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
18 INFO @ Tue, 23 Jun 2015 09:12:22: #1 read the gene table...
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
19 INFO @ Tue, 23 Jun 2015 09:12:22: #2 read the bed file of ChIP regions...
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
20 INFO @ Tue, 23 Jun 2015 09:12:22: #3 perform gene-centered annotation...
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
21 INFO @ Tue, 23 Jun 2015 09:12:22: #4 See ceas.xls for gene-centered annotation!
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
22 INFO @ Tue, 23 Jun 2015 09:12:22: #5 read the pre-computed genome bg annotation...
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
23 INFO @ Tue, 23 Jun 2015 09:12:22: #6 perform ChIP region annotation...
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
24 INFO @ Tue, 23 Jun 2015 09:12:22: #7 write a R script of ChIP region annotation...
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
25
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
26 R version 3.1.2 (2014-10-31) -- "Pumpkin Helmet"
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
27 Copyright (C) 2014 The R Foundation for Statistical Computing
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
28 Platform: x86_64-redhat-linux-gnu (64-bit)
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
29
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
30 R is free software and comes with ABSOLUTELY NO WARRANTY.
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
31 You are welcome to redistribute it under certain conditions.
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
32 Type 'license()' or 'licence()' for distribution details.
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
33
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
34 Natural language support but running in an English locale
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
35
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
36 R is a collaborative project with many contributors.
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
37 Type 'contributors()' for more information and
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
38 'citation()' on how to cite R or R packages in publications.
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
39
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
40 Type 'demo()' for some demos, 'help()' for on-line help, or
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
41 'help.start()' for an HTML browser interface to help.
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
42 Type 'q()' to quit R.
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
43
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
44 > # ARGUMENTS:
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
45 > # name = ceas
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
46 > # gene annotation table = galGal3.refGene
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
47 > # BED file = ceas_in.bed
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
48 > # WIG file = None
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
49 > # extra BED file = None
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
50 > # ChIP annotation = On
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
51 > # gene-centered annotation = On
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
52 > # average profiling = Off
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
53 > # dump profiles = Off
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
54 > # re-annotation for genome background (ChIP region annotation) = False
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
55 > # promoter sizes (ChIP region annotation) = 1000,2000,3000 bp
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
56 > # downstream sizes (ChIP region annotation) = 1000,2000,3000 bp
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
57 > # bidrectional promoter sizes (ChIP region annotation) = 2500,5000 bp
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
58 > # span size (gene-centered annotation) = 3000 bp
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
59 > pdf("ceas.pdf",height=11.5,width=8.5)
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
60 >
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
61 > # 09:12:22 Tue, 23 Jun 2015
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
62 > #
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
63 > # ChIP annotation
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
64 > #
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
65 >
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
66 >
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
67 > #
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
68 > # Chromosomal Distribution
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
69 > #
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
70 >
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
71 > par(mar=c(4, 4, 5, 3.8),oma=c(4, 2, 4, 2))
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
72 > r0<-c(100.0)
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
73 > r1<-c(100.0)
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
74 > height<-rbind(r0,r1)
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
75 > names=c("26")
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
76 > 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)
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
77 > text(x=c(100.0),y=mp[1,],label=c("100.0 %"),pos=4,offset=0.2,cex=0.9)
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
78 > text(x=c(100.0),y=mp[2,],label=c("100.0 % (<=4.9e-324)"),pos=4,offset=0.2,cex=0.9)
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
79 > legend("right",legend=c("Genome","ChIP (p-value)"),col=c("#5FA1C1","#EB9D86"),pch=15,bty="n")
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
80 >
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
81 > #
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
82 > # Promoter,Bipromoter,Downstream, Gene and Regions of interest
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
83 > #
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
84 >
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
85 > par(mfrow=c(4, 1),mar=c(4, 4, 5, 3.8),oma=c(4, 2, 4, 2))
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
86 > r0<-c(1.8532425688606797, 3.616851183410451, 5.322318854623416)
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
87 > r1<-c(0.0, 0.0, 0.0)
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
88 > height<-rbind(r0,r1)
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
89 > names=c("<=1000 bp","<=2000 bp","<=3000 bp")
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
90 > 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)
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
91 > 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)
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
92 > text(x=mp[2,],y=c(0.0, 0.0, 0.0),label=c("0.000 %
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
93 + (0.981)","0.000 %
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
94 + (0.964)","0.000 %
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
95 + (0.947)"),pos=3,offset=0.2)
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
96 > legend("topleft",legend=c("Genome","ChIP (p-value)"),col=c("#5FA1C1","#EB9D86"),pch=15,bty="n")
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
97 > r0<-c(0.03876062889120376, 0.03876062889120376)
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
98 > r1<-c(0.0, 0.0)
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
99 > height<-rbind(r0,r1)
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
100 > names=c("<=2500 bp","<=5000 bp")
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
101 > 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)
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
102 > text(x=mp[1,],y=c(0.03876062889120376, 0.03876062889120376),label=c("0.04 %","0.04 %"),pos=3,offset=0.2)
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
103 > text(x=mp[2,],y=c(0.0, 0.0),label=c("0.000 %
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
104 + (1.000)","0.000 %
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
105 + (1.000)"),pos=3,offset=0.2)
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
106 > legend("topleft",legend=c("Genome","ChIP (p-value)"),col=c("#5FA1C1","#EB9D86"),pch=15,bty="n")
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
107 > r0<-c(1.8290171758036773, 3.4690762857627364, 4.980740812519683)
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
108 > r1<-c(0.0, 0.0, 0.0)
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
109 > height<-rbind(r0,r1)
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
110 > names=c("<=1000 bp","<=2000 bp","<=3000 bp")
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
111 > 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)
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
112 > 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)
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
113 > text(x=mp[2,],y=c(0.0, 0.0, 0.0),label=c("0.000 %
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
114 + (0.982)","0.000 %
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
115 + (0.965)","0.000 %
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
116 + (0.950)"),pos=3,offset=0.2)
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
117 > legend("topleft",legend=c("Genome","ChIP (p-value)"),col=c("#5FA1C1","#EB9D86"),pch=15,bty="n")
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
118 > r0<-c(0.2034933016788197, 1.3978051793890356, 2.359553283752029, 19.734005184234114, 23.694856949054)
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
119 > r1<-c(0.0, 0.0, 0.0, 0.0, 0.0)
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
120 > height<-rbind(r0,r1)
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
121 > names=c("5'UTR","3'UTR","Coding Exon","Intron","All")
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
122 > 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)
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
123 > 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)
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
124 > text(x=mp[2,],y=c(0.0, 0.0, 0.0, 0.0, 0.0),label=c("0.000 %
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
125 + (0.998)","0.000 %
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
126 + (0.986)","0.000 %
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
127 + (0.976)","0.000 %
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
128 + (0.803)","0.000 %
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
129 + (0.763)"),pos=3,offset=0.2)
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
130 > legend("topleft",legend=c("Genome","ChIP (p-value)"),col=c("#5FA1C1","#EB9D86"),pch=15,bty="n")
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
131 >
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
132 > #
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
133 > # Distribution of Genome and ChIP regions over cis-regulatory element
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
134 > # Note that the x may be modified for better graphics in case a value is too small
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
135 > # Thus, look at the labels of the pie chart to get the real percentage values
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
136 > #
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
137 >
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
138 > par(mfcol=c(2, 2),mar=c(3, 3, 4, 2.8),oma=c(4, 2, 4, 2))
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
139 > x<-c(0.018532,0.017055,0.016037,0.017830,0.015092,0.014051,0.010000,0.013833,0.023014,0.192592,0.670292)
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
140 > 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)
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
141 > x<-c(0.000000,1.000000)
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
142 > y<-c(0.000000,1.000000)
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
143 > plot(x, y,type="n",main="",xlab="",ylab="",frame=FALSE,axes=FALSE,xaxt="s",yaxt="s")
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
144 > 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")
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
145 > x<-c(0.010000,0.010000,0.010000,0.010000,0.010000,0.010000,0.010000,0.010000,0.010000,0.010000,1.000000)
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
146 > 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)
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
147 > x<-c(0.000000,1.000000)
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
148 > y<-c(0.000000,1.000000)
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
149 > plot(x, y,type="n",main="",xlab="",ylab="",frame=FALSE,axes=FALSE,xaxt="s",yaxt="s")
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
150 > 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")
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
151 >
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
152 > #
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
153 > # ChIP regions over the genome
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
154 > #
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
155 >
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
156 > par(mar=c(4, 4, 5, 3.8),oma=c(4, 2, 4, 2))
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
157 > layout(matrix(c(1, 0, 2, 2), 2, 2, byrow = TRUE),widths=c(1, 1),heights=c(1, 5))
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
158 > x<-c(0.000000,2.515610)
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
159 > y<-c(0.000000,1.000000)
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
160 > plot(x, y,type="n",main="Distribution of Peak Heights",xlab="",ylab="",xlim=c(0.000000,2.515610),ylim=c(0.000000,1.000000),frame=FALSE,xaxt="s",yaxt="n",cex=0.9)
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
161 > x<-c(0.000000,2.515610,2.515610,0.000000)
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
162 > y<-c(0.000000,0.000000,1.000000,1.000000)
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
163 > polygon(x,y,col=c("black"))
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
164 > x <- c(0.000000,0.169726,0.339451,0.509177,0.678903,0.848628,1.018354,1.188079,1.357805,1.527531,1.697256,1.866982,2.036708,2.206433,2.376159)
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
165 > 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.800000)
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
166 > lines(x, y,xlim=c(0, 2.51561),ylim=c(0, 1),type="l",col=c("cyan"),lwd=2)
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
167 > x<-c(4119129.000000,4119130.000000)
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
168 > y<-c(0.855556,1.144444)
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
169 > plot(x, y,type="n",main="ChIP Regions (Peaks) over Chromosomes",xlab="Chromosome Size (bp)",ylab="Chromosome",xlim=c(4119129.000000,4119130.000000),ylim=c(0.855556,1.144444),frame=FALSE,xaxt="s",yaxt="n")
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
170 > start <- c(4119129)
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
171 > end <- c(4119130)
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
172 > vals <- c(2.51561)
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
173 > vals[vals > 2.51561] <- 2.51561
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
174 > vals[vals < 0] <- 0
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
175 > heights <- 0.288889 * ((vals - 0)/(2.51561 - 0)) + 0.855555555556
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
176 > for (i in 1:length(heights)) {
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
177 + 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"))
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
178 + }
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
179 > mtext("26",side=2,line=0,outer=FALSE,at=1.0)
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
180 > dev.off()
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
181 null device
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
182 1
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
183 >
f411ce97a351 Uploaded initial version 1.0.2-2
pjbriggs
parents:
diff changeset
184 INFO @ Tue, 23 Jun 2015 09:12:22: #... cong! See ceas.pdf for the graphical results of CEAS!