# HG changeset patch # User pravs # Date 1529544817 14400 # Node ID 80892c607b1e33e6ae48dd643aae85dbdd188987 # Parent fdd5707c61278cf37c9826e4013a0d2622b1b3b6 planemo upload diff -r fdd5707c6127 -r 80892c607b1e protein_rna_correlation.r --- a/protein_rna_correlation.r Wed Jun 20 21:23:00 2018 -0400 +++ b/protein_rna_correlation.r Wed Jun 20 21:33:37 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 fdd5707c6127 -r 80892c607b1e protein_rna_correlation.xml --- a/protein_rna_correlation.xml Wed Jun 20 21:23:00 2018 -0400 +++ b/protein_rna_correlation.xml Wed Jun 20 21:33:37 2018 -0400 @@ -3,8 +3,8 @@ R r-protrnacorr - ghostscript - + + packagelibx diff -r fdd5707c6127 -r 80892c607b1e tool_dependencies.xml --- a/tool_dependencies.xml Wed Jun 20 21:23:00 2018 -0400 +++ b/tool_dependencies.xml Wed Jun 20 21:33:37 2018 -0400 @@ -38,10 +38,10 @@ Contains a tool dependency definition that downloads and installs protein rna correlation tool package for R library. - - - -