Mercurial > repos > greg > ideas

changeset 165:bb5544d1c85e draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author greg
date Thu, 18 Jan 2018 13:47:35 -0500
parents 68d909206a15
children 7d33a8157c6b
files ideas.xml
diffstat 1 files changed, 87 insertions(+), 127 deletions(-) [+]
line wrap: on
line diff
--- a/ideas.xml	Fri Jan 12 14:56:24 2018 -0500
+++ b/ideas.xml	Thu Jan 18 13:47:35 2018 -0500
@@ -17,6 +17,7 @@
 #set prep_input_config = "prep_input_config.txt"
 #set prep_output_config = "prep_output_config.txt"
 #set specify_genomic_window = $specify_genomic_window_cond.specify_genomic_window
+#set windows_positions_by_chroms_config = "windows_positions_by_chroms_config.txt"
 #set perform_training = $perform_training_cond.perform_training
 ##############################################
 ## Create the config file and prepare the data
@@ -25,9 +26,12 @@
     mkdir '$output_pdf_dir' &&
 #end if
 #if str($perform_training) == "yes":
+    #set output_dir = $output_training_dir
     mkdir '$output_training_dir' &&
+#else:
+    #set output_dir = $output_txt_dir
+    mkdir '$output_txt_dir' &&
 #end if
-mkdir '$output_txt_dir' &&
 cp '$gen_prep_input_config' $prep_input_config &&
 sort $prep_input_config -o $prep_input_config &&
 prepMat
@@ -55,7 +59,7 @@
 $norm
 &>prepmat_log.txt;
 if [[ $? -ne 0 ]]; then
-    cp prepmat_log.txt '$output_txt_dir';
+    cp prepmat_log.txt '$output_dir';
     exit 1;
 fi
 ##############################################
@@ -71,133 +75,68 @@
     ## the window positions by chromosome to enable
     ## the IDEAS -inv option.
     ##############################################
-    #import collections
-    #set window_positions_by_chroms = collections.OrderedDict()
-    #for count, line in enumerate(open($specify_genomic_window_cond.bed_input.file_name, 'r')):
-        #set $line = $line.strip()
-        #if not $line or $line.startswith('#'):
-            #continue
-        #end if
-        #set items = $line.split('\t')
-        #if $items[0] in $window_positions_by_chroms:
-            #set tup = $window_positions_by_chroms[$items[0]]
-            #set $tup[1] += 1
-            #set $window_positions_by_chroms[$items[0]] = $tup
-        #else:
-            #set $window_positions_by_chroms[$items[0]] = [$count, $count+1]
-        #end if
-    #end for
-    #for chrom, tup in $window_positions_by_chroms.items():
-        && Rscript '$__tool_directory__/ideas.R'
-        #if str($perform_training) == "yes":
-            --training_iterations $perform_training_cond.training_iterations
-            --training_windows $perform_training_cond.training_windows
-        #end if
-        --prep_output_config '$prep_output_config'
-        --windows_bed '$specify_genomic_window_cond.bed_input'
-        --hp $hp
-        --window_start $tup[0]
-        --window_end $tup[1]
-        #if str($log2) != "0.0":
-            --log2 $log2
-        #end if
-        #if str($max_states) != "0.0":
-            --max_states $max_states
-        #end if
-        #if str($initial_states) != "0":
-            --initial_states $initial_states
-        #end if
-        #if str($max_position_classes) != "0":
-            --max_position_classes $max_position_classes
-        #end if
-        #if str($max_cell_type_clusters) != "0":
-            --max_cell_type_clusters $max_cell_type_clusters
-        #end if
-        #if str($prior_concentration) != "0.0":
-            --prior_concentration $prior_concentration
-        #end if
-        --burnin_num $burnin_num
-        --mcmc_num $mcmc_num
-        #if str($minerr) != "0.0":
-            --minerr $minerr
-        #end if
-        #if str($maxerr) != "0.0":
-            --maxerr $maxerr
-        #end if
-        --rseed $rseed
-        --thread \${GALAXY_SLOTS:-4}
-        --project_name '$project_name.$chrom'
-        #if str($save_ideas_log) == "yes":
-            --save_ideas_log $save_ideas_log
-            --output_log '$output_log'
-        #end if
-    #end for
-#else:
-    ##############################################
-    ## Not using a genomic window bed file.
-    ##############################################
-    && Rscript '$__tool_directory__/ideas.R'
-    #if str($perform_training) == "yes":
-        --training_iterations $perform_training_cond.training_iterations
-        --training_windows $perform_training_cond.training_windows
-    #end if
-    --prep_output_config '$prep_output_config'
-    --windows_bed $tmp_dir/*.bed
-    --hp $hp
-    #if str($log2) != "0.0":
-        --log2 $log2
-    #end if
-    #if str($max_states) != "0.0":
-        --max_states $max_states
-    #end if
-    #if str($initial_states) != "0":
-        --initial_states $initial_states
-    #end if
-    #if str($max_position_classes) != "0":
-        --max_position_classes $max_position_classes
-    #end if
-    #if str($max_cell_type_clusters) != "0":
-        --max_cell_type_clusters $max_cell_type_clusters
-    #end if
-    #if str($prior_concentration) != "0.0":
-        --prior_concentration $prior_concentration
-    #end if
-    --burnin_num $burnin_num
-    --mcmc_num $mcmc_num
-    #if str($minerr) != "0.0":
-        --minerr $minerr
-    #end if
-    #if str($maxerr) != "0.0":
-        --maxerr $maxerr
-    #end if
-    --rseed $rseed
-    --thread \${GALAXY_SLOTS:-4}
-    --project_name '$project_name'
-    #if str($save_ideas_log) == "yes":
-        --save_ideas_log $save_ideas_log
-        --output_log '$output_log'
-    #end if
+    && cp '$gen_windows_positions_by_chroms_config' $windows_positions_by_chroms_config
+#end if
+&& Rscript '$__tool_directory__/ideas.R'
+#if str($specify_genomic_window) == "yes":
+    --windows_bed '$specify_genomic_window_cond.bed_input'
+    --windows_config $windows_positions_by_chroms_config
 #end if
 #if str($perform_training) == "yes":
-    && mv ./*.para0 '$output_training_dir'
-    && mv ./*.profile0 '$output_training_dir'
-    #if str($output_heatmaps) == "yes":
-        && Rscript '$__tool_directory__/create_heatmaps.R'
-        --input_dir 'output_training_dir'
-        --output_dir '$output_pdf_dir'
-        --script_dir '$__tool_directory__'
-        --in_training_mode true
-    #end if
+    --training_iterations $perform_training_cond.training_iterations
+    --training_windows $perform_training_cond.training_windows
+#end if
+--prep_output_config '$prep_output_config'
+--hp $hp
+#if str($log2) != "0.0":
+    --log2 $log2
+#end if
+#if str($max_states) != "0.0":
+    --max_states $max_states
+#end if
+#if str($initial_states) != "0":
+    --initial_states $initial_states
+#end if
+#if str($max_position_classes) != "0":
+    --max_position_classes $max_position_classes
+#end if
+#if str($max_cell_type_clusters) != "0":
+    --max_cell_type_clusters $max_cell_type_clusters
+#end if
+#if str($prior_concentration) != "0.0":
+    --prior_concentration $prior_concentration
+#end if
+--burnin_num $burnin_num
+--mcmc_num $mcmc_num
+#if str($minerr) != "0.0":
+    --minerr $minerr
+#end if
+#if str($maxerr) != "0.0":
+    --maxerr $maxerr
+#end if
+--rseed $rseed
+--thread \${GALAXY_SLOTS:-4}
+--project_name '$project_name'
+#if str($save_ideas_log) == "yes":
+    --save_ideas_log $save_ideas_log
+    --output_log '$output_log'
+#end if
+#if str($perform_training) == "yes":
+    && mv ./*.para0 '$output_dir'
+    && mv ./*.profile0 '$output_dir'
 #else:
-    && mv ./*.cluster '$output_txt_dir'
-    && mv ./*.para '$output_txt_dir'
-    && mv ./*.profile '$output_txt_dir'
-    && mv ./*.state '$output_txt_dir'
-    #if str($output_heatmaps) == "yes":
-        && Rscript '$__tool_directory__/create_heatmaps.R'
-        --input_dir '$output_txt_dir'
-        --output_dir '$output_pdf_dir'
-        --script_dir '$__tool_directory__'
+    && mv ./*.para '$output_dir'
+    && mv ./*.profile '$output_dir'
+#end if
+&& mv ./*.cluster '$output_dir'
+&& mv ./*.state '$output_dir'
+#if str($output_heatmaps) == "yes":
+    && Rscript '$__tool_directory__/create_heatmaps.R'
+    --input_dir '$output_dir'
+    --output_dir '$output_pdf_dir'
+    --script_dir '$__tool_directory__'
+    #if str($perform_training) == "yes":
+        --in_training_mode true
     #end if
 #end if
     ]]></command>
@@ -225,6 +164,27 @@
 ${input_items.cell_type_name} ${input_items.epigenetic_factor_name} ${input_items.input}
     #end for
 #end if]]></configfile>
+        <configfile name="gen_windows_positions_by_chroms_config"><![CDATA[#if str($specify_genomic_window_cond.specify_genomic_window) == "yes":
+    #import collections
+    #set window_positions_by_chroms_odict = $collections.OrderedDict()
+    #for count, line in enumerate(open($specify_genomic_window_cond.bed_input.file_name, 'r')):
+        #set $line = $line.strip()
+        #if not $line or $line.startswith('#'):
+            #continue
+        #end if
+        #set items = $line.split('\t')
+        #if $items[0] in $window_positions_by_chroms_odict:
+            #set tup = $window_positions_by_chroms_odict[$items[0]]
+            #set $tup[1] += 1
+            #set $window_positions_by_chroms_odict[$items[0]] = $tup
+        #else:
+            #set $window_positions_by_chroms_odict[$items[0]] = [$count, $count+1]
+        #end if
+    #end for
+    #for $chrom, $tup in $window_positions_by_chroms_odict.items():
+${chrom} ${tup[0]} ${tup[1]}
+    #end for
+#end if]]></configfile>
     </configfiles>
     <inputs>
         <conditional name="perform_training_cond">
@@ -233,7 +193,7 @@
                 <option value="no">No</option>
             </param>
             <when value="yes">
-                <param name="training_iterations" type="integer" value="20" min="2" label="Number of training iterations"/>
+                <param name="training_iterations" type="integer" value="20" min="3" label="Number of training iterations"/>
                 <param name="training_windows" type="integer" value="10000" min="2" label="Number of randomly selected windows for training"/>
             </when>
             <when value="no"/>