Mercurial > repos > greg > multigps
comparison multigps.xml @ 25:32a59270224b draft
Uploaded
author | greg |
---|---|
date | Wed, 14 Dec 2016 07:54:08 -0500 |
parents | d5c1ab05259b |
children | fcc1cc3f3e3c |
comparison
equal
deleted
inserted
replaced
24:d5c1ab05259b | 25:32a59270224b |
---|---|
28 #set gmsc = $aoc.gauss_model_smoothing_cond | 28 #set gmsc = $aoc.gauss_model_smoothing_cond |
29 #set rbec = $aoc.report_binding_events_cond | 29 #set rbec = $aoc.report_binding_events_cond |
30 #set rloc = $aoc.reads_limits_options_cond | 30 #set rloc = $aoc.reads_limits_options_cond |
31 #set sdc = $aoc.scale_data_cond | 31 #set sdc = $aoc.scale_data_cond |
32 #set umc = $aoc.use_motif_cond | 32 #set umc = $aoc.use_motif_cond |
33 #set mpc = $umc.multigps_priors_cond | |
34 #set bmc = $mpc.both_motifs_cond | |
35 #set mfoc = $bmc.nomotifprior_cond | |
36 #set rgc = $umc.reference_genome_cond | |
37 #if str($aoc.advanced_options) == 'display': | 33 #if str($aoc.advanced_options) == 'display': |
38 ############################ | 34 ############################ |
39 ## Limits on how many reads | 35 ## Limits on how many reads |
40 ############################ | 36 ############################ |
41 #if str($rloc.reads_limits) == 'yes': | 37 #if str($rloc.reads_limits) == 'yes': |
105 --exclude "$aoc.exclude" | 101 --exclude "$aoc.exclude" |
106 ############################ | 102 ############################ |
107 ## MultiGPS priors | 103 ## MultiGPS priors |
108 ############################ | 104 ############################ |
109 #if str($umc.use_motif) == 'yes': | 105 #if str($umc.use_motif) == 'yes': |
106 #set mpc = $umc.multigps_priors_cond | |
107 #set rgc = $umc.reference_genome_cond | |
110 #if str($rgc.reference_genome_source) == "cached": | 108 #if str($rgc.reference_genome_source) == "cached": |
111 #set seq_dir = os.path.split($rgc.reference_genome.fields.path)[0] | 109 #set seq_dir = os.path.split($rgc.reference_genome.fields.path)[0] |
112 #else: | 110 #else: |
113 ## MultiGPS requires a directory containing reference files, so symlink the history dataset. | 111 ## MultiGPS requires a directory containing reference files, so symlink the history dataset. |
114 #import os | 112 #import os |
120 ln -f -s $tmp_seq_file $seq_file && | 118 ln -f -s $tmp_seq_file $seq_file && |
121 #end if | 119 #end if |
122 --seq $seq_dir | 120 --seq $seq_dir |
123 --geninfo ${chromInfo} | 121 --geninfo ${chromInfo} |
124 #if str($mpc.multigps_priors) == 'yes': | 122 #if str($mpc.multigps_priors) == 'yes': |
123 #set bmc = $mpc.both_motifs_cond | |
125 #if str($mpc.noposprior) == 'no': | 124 #if str($mpc.noposprior) == 'no': |
126 --noposprior | 125 --noposprior |
127 #end if | 126 #end if |
128 --probshared $mpc.probshared | 127 --probshared $mpc.probshared |
129 #if str($bmc.nomotifs) == 'yes': | 128 #if str($bmc.nomotifs) == 'yes': |
130 --memenmotifs $bmc.memenmotifs | 129 --memenmotifs $bmc.memenmotifs |
131 --mememinw $bmc.mememinw | 130 --mememinw $bmc.mememinw |
132 --mememaxw $bmc.mememaxw | 131 --mememaxw $bmc.mememaxw |
133 #else: | 132 #else: |
133 #set mfoc = $bmc.nomotifprior_cond | |
134 --nomotifs | 134 --nomotifs |
135 --nomotifprior $mfoc.nomotifprior | 135 --nomotifprior $mfoc.nomotifprior |
136 #if str($mfoc.nomotifprior) == 'yes': | 136 #if str($mfoc.nomotifprior) == 'yes': |
137 --memenmotifs $mfoc.memenmotifs | 137 --memenmotifs $mfoc.memenmotifs |
138 --mememinw $mfoc.mememinw | 138 --mememinw $mfoc.mememinw |