Skip to content

Commit a61c8fc

Browse files
committed
Update length plot to improve poster readibility
1 parent 3fcfec4 commit a61c8fc

File tree

1 file changed

+6
-4
lines changed

1 file changed

+6
-4
lines changed

length_plot.r

Lines changed: 6 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -133,8 +133,10 @@ dens.mat <- sapply(fileNames, function(x){
133133
plain.hist(data,
134134
main=sprintf("Read Count Distribution Plot (%s)",subName));
135135
dummy <- dev.off();
136+
par(mar=c(5.5,6.5,2.5,1.5));
136137
sequence.hist(data,
137138
main=sprintf("Read Length Distribution Plot (%s)",subName));
139+
par(mar=c(5.5,1.5,2.5,6.5));
138140
plain.hist(data,
139141
main=sprintf("Read Count Distribution Plot (%s)",subName));
140142
res.out;
@@ -300,10 +302,10 @@ bpdens.mat <- dens.mat * 10^as.numeric(rownames(dens.mat));
300302
ylab = "");
301303
for(col in seq_len(ncol(bpdens.mat))){
302304
points(as.numeric(rownames(bpdens.mat)),bpdens.mat[,col] /
303-
sum(bpdens.mat[,col]), type="l", lwd=3,
305+
sum(bpdens.mat[,col]), type="l", lwd=5,
304306
col=hcl(h=col/ncol(bpdens.mat)*360, l=70, c=80));
305307
}
306-
legend("topleft", cex=0.71, legend=colnames(bpdens.mat), ncol=ceiling(ncol(bpdens.mat)/16),
308+
legend("topleft", cex=1.5, legend=colnames(bpdens.mat), ncol=ceiling(ncol(bpdens.mat)/16),
307309
fill=hcl(h=(1:ncol(bpdens.mat))/ncol(bpdens.mat)*360, l=70, c=80),
308310
inset=0.05);
309311
mtext("Base Density (arbitrary scale)",2,3);
@@ -382,10 +384,10 @@ bpdens.mat <- dens.mat * 10^as.numeric(rownames(dens.mat));
382384
ylab = "");
383385
for(col in seq_len(ncol(bpdens.mat))){
384386
points(as.numeric(rownames(bpdens.mat)),bpdens.mat[,col] /
385-
sum(bpdens.mat[,col]), type="l", lwd=3,
387+
sum(bpdens.mat[,col]), type="l", lwd=4,
386388
col=hcl(h=col/ncol(bpdens.mat)*360, l=70, c=80));
387389
}
388-
legend("topleft", cex=0.71, legend=colnames(bpdens.mat), ncol=ceiling(ncol(bpdens.mat)/16),
390+
legend("topleft", cex=1.5, legend=colnames(bpdens.mat), ncol=ceiling(ncol(bpdens.mat)/16),
389391
fill=hcl(h=(1:ncol(bpdens.mat))/ncol(bpdens.mat)*360, l=70, c=80),
390392
inset=0.05);
391393
mtext("Base Density (arbitrary scale)",2,3);

0 commit comments

Comments
 (0)