# HG changeset patch
# User pravs
# Date 1529545356 14400
# Node ID e407b1a7a8de4438a349357d22884cae2f4935f2
# Parent 80892c607b1e33e6ae48dd643aae85dbdd188987
planemo upload
diff -r 80892c607b1e -r e407b1a7a8de protein_rna_correlation.r
--- a/protein_rna_correlation.r Wed Jun 20 21:33:37 2018 -0400
+++ b/protein_rna_correlation.r Wed Jun 20 21:42:36 2018 -0400
@@ -638,29 +638,29 @@
file = htmloutfile, append = TRUE);
outplot = paste(outdir,"/PE_GE_lm_1.png",sep="",collapse="");
- png(outplot);
- #bitmap(outplot,"png16m");
+ #png(outplot);
+ bitmap(outplot,"png16m");
par(mfrow=c(1,1));
plot(regmodel, 1, cex.lab=1.5);
dev.off();
outplot = paste(outdir,"/PE_GE_lm_2.png",sep="",collapse="");
- png(outplot);
- #bitmap(outplot,"png16m");
+ #png(outplot);
+ bitmap(outplot,"png16m");
par(mfrow=c(1,1));
plot(regmodel, 2, cex.lab=1.5);
dev.off();
outplot = paste(outdir,"/PE_GE_lm_3.png",sep="",collapse="");
- png(outplot);
- #bitmap(outplot,"png16m");
+ #png(outplot);
+ bitmap(outplot,"png16m");
par(mfrow=c(1,1));
plot(regmodel, 3, cex.lab=1.5);
dev.off();
outplot = paste(outdir,"/PE_GE_lm_5.png",sep="",collapse="");
- png(outplot);
- #bitmap(outplot,"png16m");
+ #png(outplot);
+ bitmap(outplot,"png16m");
par(mfrow=c(1,1));
plot(regmodel, 5, cex.lab=1.5);
dev.off();
@@ -711,8 +711,8 @@
cat("Generating cooksd plot\n",file=logfile, append=T);
outplot = paste(outdir,"/PE_GE_lm_cooksd.png",sep="",collapse="");
- png(outplot);
- #bitmap(outplot,"png16m");
+ #png(outplot);
+ bitmap(outplot,"png16m");
par(mfrow=c(1,1));
plot(cooksd, pch="*", cex=2, cex.lab=1.5,main="Influential Obs. by Cooks distance", ylab="Cook\'s distance", xlab="Observations") # plot cooks distance
abline(h = 4*mean(cooksd, na.rm=T), col="red") # add cutoff line
@@ -772,8 +772,8 @@
# Scatter plot
#=============================================================================================================
outplot = paste(outdir,"/AbundancePlot_scatter_without_outliers.png",sep="",collapse="");
- png(outplot);
- #bitmap(outplot,"png16m");
+ #png(outplot);
+ bitmap(outplot,"png16m");
par(mfrow=c(1,1));
scatter.smooth(PE_GE_data_no_outlier[,"GE_abundance"], PE_GE_data_no_outlier[,"PE_abundance"], xlab="Transcript Abundance", ylab="Protein Abundance", cex.lab=1.5);
@@ -815,8 +815,8 @@
cat("Generating heatmap plot\n",file=logfile, append=T);
outplot = paste(outdir,"/PE_GE_heatmap.png",sep="",collapse="");
- png(outplot);
- #bitmap(outplot,"png16m");
+ #png(outplot);
+ bitmap(outplot,"png16m");
par(mfrow=c(1,1));
#heatmap.2(as.matrix(PE_GE_data[,c("PE_abundance","GE_abundance")]), trace="none", cexCol=1, col=greenred(100),Colv=F, labCol=c("PE","GE"), scale="col");
my_palette <- colorRampPalette(c("green", "white", "red"))(299);
@@ -837,8 +837,8 @@
k1 = kmeans(PE_GE_data_kdata[,c("PE_abundance","GE_abundance")], 5);
cat("Generating kmeans plot\n",file=logfile, append=T);
outplot = paste(outdir,"/PE_GE_kmeans.png",sep="",collapse="");
- png(outplot);
- #bitmap(outplot,"png16m");
+ #png(outplot);
+ bitmap(outplot,"png16m");
par(mfrow=c(1,1));
scatter.smooth(PE_GE_data_kdata[,"GE_abundance"], PE_GE_data_kdata[,"PE_abundance"], xlab="Transcript Abundance", ylab="Protein Abundance", cex.lab=1.5);
diff -r 80892c607b1e -r e407b1a7a8de protein_rna_correlation.xml
--- a/protein_rna_correlation.xml Wed Jun 20 21:33:37 2018 -0400
+++ b/protein_rna_correlation.xml Wed Jun 20 21:42:36 2018 -0400
@@ -3,8 +3,8 @@
R
r-protrnacorr
-
- packagelibx
+ ghostscript
+
diff -r 80892c607b1e -r e407b1a7a8de tool_dependencies.xml
--- a/tool_dependencies.xml Wed Jun 20 21:33:37 2018 -0400
+++ b/tool_dependencies.xml Wed Jun 20 21:42:36 2018 -0400
@@ -38,10 +38,10 @@
Contains a tool dependency definition that downloads and installs protein rna correlation tool package for R library.
-
+
+
+
-