# HG changeset patch # User jeremyjliu # Date 1432658648 14400 # Node ID 4803f5186f1a3e576fcc6a5d6c519b9009920ff4 # Parent 53e45130e6f9a0cfbb5154595fc00f8967f5b72f Uploaded diff -r 53e45130e6f9 -r 4803f5186f1a region_motif_compare.r --- a/region_motif_compare.r Sat May 16 22:38:25 2015 -0400 +++ b/region_motif_compare.r Tue May 26 12:44:08 2015 -0400 @@ -28,8 +28,8 @@ # Set input and reference files inTab1 = args[9] inTab2 = args[10] -enrichTab = args[11] -depleteTab = args[12] +enrichTab1 = args[11] +enrichTab2 = args[12] plotsPng = args[13] # Load dependencies @@ -174,9 +174,9 @@ outDF$enrichment = ifelse(outDF$enrichment>1, round(outDF$enrichment,3), paste("1/",round(1/outDF$enrichment,3))) - write.table(outDF[indicesEnrich,], file=enrichTab, quote=FALSE, + write.table(outDF[indicesEnrich,], file=enrichTab1, quote=FALSE, sep="\t", append=FALSE, row.names=FALSE, col.names=TRUE) - write.table(outDF[indicesDeplete,], file=depleteTab, quote=FALSE, + write.table(outDF[indicesDeplete,], file=enrichTab2, quote=FALSE, sep="\t", append=FALSE, row.names=FALSE, col.names=TRUE) } diff -r 53e45130e6f9 -r 4803f5186f1a region_motif_compare.xml --- a/region_motif_compare.xml Sat May 16 22:38:25 2015 -0400 +++ b/region_motif_compare.xml Tue May 26 12:44:08 2015 -0400 @@ -4,14 +4,14 @@ RMOTIF_PATH - region_motif_compare.r --args \$RMOTIF_PATH $pwm_file.fields.path $in_tab_1 $in_tab_2 - $out_enriched - $out_depleted + region_motif_compare.r --args \$RMOTIF_PATH $pwm_file.fields.path $in_tab_a $in_tab_b + $out_enriched_a + $out_enriched_b $out_plots - - + + @@ -19,11 +19,12 @@ - - + + + This tools reads in two counts file and determines enriched and depleted