# HG changeset patch # User davidecangelosi # Date 1558020022 14400 # Node ID 3168db2e0ff5a5751f7ecf34c4311e79c0234cb6 # Parent 969d0bbbf5a9cab009c926e169a69e3f9e2df3af planemo upload for repository https://github.com/igg-molecular-biology-lab/pipe-t.git commit 4e1ae25e30206be75710403c7500980ddb474beb diff -r 969d0bbbf5a9 -r 3168db2e0ff5 pipe-t.R --- a/pipe-t.R Thu May 16 11:12:34 2019 -0400 +++ b/pipe-t.R Thu May 16 11:20:22 2019 -0400 @@ -795,8 +795,8 @@ PY = ecdf(gm) plot_colors <- c(rgb(r=0.0,g=0.0,b=0.9), "red", "forestgreen",rgb(r=0.0,g=0.0,b=0.0),rgb(r=0.5,g=0.0,b=0.3),rgb(r=0.0,g=0.4,b=0.4)) - -plot(P,cex.lab=3, cex.axis = 2,col=plot_colors[1],xlim=c(0.0,600), mar = c(8,8,8,8), ylim=c(0.0,1),xaxp = c(0.0, 600, 6),yaxp = c(0.0, 1, 10), cex=1.3, lwd=5, main=paste("p-value=", formatC(ks.test(vec,gm)$p.value, format = "e", digits = 2)),xlab="CV(%)",ylab="Empirical Cumulative Distribution") +par(mar = c(8,8,8,8)) +plot(P,cex.lab=3, cex.axis = 2,col=plot_colors[1],xlim=c(0.0,600), ylim=c(0.0,1),xaxp = c(0.0, 600, 6),yaxp = c(0.0, 1, 10), cex=1.3, lwd=5, main=paste("p-value=", formatC(ks.test(vec,gm)$p.value, format = "e", digits = 2)),xlab="CV(%)",ylab="Empirical Cumulative Distribution") lines(PY, lwd=5, col=plot_colors[6],cex=1.3) legend("bottomright", c("not normalized", "normalized"), cex=1.3, col=c(plot_colors[1],plot_colors[6]), lwd=c(5,5)); dev.off()