|
70 | 70 | mkdir($checkpoints_dir) or die $!; |
71 | 71 | } |
72 | 72 |
|
73 | | - my $pipeliner = new Pipeliner(-verbose => 2); |
74 | | - |
| 73 | + my $pipeliner = new Pipeliner(-verbose => 2, -checkpoint_dir => $checkpoints_dir); |
| 74 | + |
75 | 75 | my $cmd = "$utildir/build_atgPWM_+-.pl " |
76 | 76 | . " --transcripts $transcripts_file " |
77 | 77 | . " --selected_orfs $selected_orfs_file " |
78 | 78 | . " --out_prefix $out_prefix" |
79 | 79 | . " --pwm_left $pwm_left --pwm_right $pwm_right "; |
80 | 80 |
|
81 | | - $pipeliner->add_commands(new Command($cmd, "$checkpoints_dir/built_init.ok")); |
| 81 | + $pipeliner->add_commands(new Command($cmd, "built_init.ok")); |
82 | 82 |
|
83 | 83 | $cmd = "$utildir/feature_scoring.+-.pl " |
84 | 84 | . " --features_plus $out_prefix.+.features " |
85 | 85 | . " --features_minus $out_prefix.-.features " |
86 | 86 | . " --atg_position $pwm_left " |
87 | 87 | . " > $out_prefix.feature.scores"; |
88 | 88 |
|
89 | | - $pipeliner->add_commands(new Command($cmd, "$checkpoints_dir/score_features.ok")); |
| 89 | + $pipeliner->add_commands(new Command($cmd, "score_features.ok")); |
90 | 90 |
|
91 | 91 | $cmd = "$utildir/feature_scores_to_ROC.pl $out_prefix.feature.scores > $out_prefix.feature.scores.roc"; |
92 | | - $pipeliner->add_commands(new Command($cmd, "$checkpoints_dir/roc_features.ok")); |
| 92 | + $pipeliner->add_commands(new Command($cmd, "roc_features.ok")); |
93 | 93 |
|
94 | 94 | $cmd = "$utildir/plot_ROC.Rscript $out_prefix.feature.scores.roc || :"; |
95 | | - $pipeliner->add_commands(new Command($cmd, "$checkpoints_dir/rocplot.ok")); |
| 95 | + $pipeliner->add_commands(new Command($cmd, "rocplot.ok")); |
96 | 96 |
|
97 | 97 | $cmd = "$utildir/compute_AUC.pl $out_prefix.feature.scores.roc"; |
98 | | - $pipeliner->add_commands(new Command($cmd, "$checkpoints_dir/aucplot.ok")); |
| 98 | + $pipeliner->add_commands(new Command($cmd, "aucplot.ok")); |
99 | 99 |
|
100 | 100 | $cmd = "$utildir/make_seqLogo.Rscript $out_prefix.+.pwm || :"; |
101 | | - $pipeliner->add_commands(new Command($cmd, "$checkpoints_dir/seqlogo.+.ok")); |
| 101 | + $pipeliner->add_commands(new Command($cmd, "seqlogo.+.ok")); |
102 | 102 |
|
103 | 103 | $cmd = "$utildir/make_seqLogo.Rscript $out_prefix.-.pwm || :"; |
104 | | - $pipeliner->add_commands(new Command($cmd, "$checkpoints_dir/seqlogo.-.ok")); |
| 104 | + $pipeliner->add_commands(new Command($cmd, "seqlogo.-.ok")); |
105 | 105 |
|
106 | 106 |
|
107 | 107 | ################# |
|
113 | 113 | . " --pwm_minus $out_prefix.-.pwm " |
114 | 114 | . " --out_prefix $out_prefix.enhanced"; |
115 | 115 |
|
116 | | - $pipeliner->add_commands(new Command($cmd, "$checkpoints_dir/enhance.ok")); |
| 116 | + $pipeliner->add_commands(new Command($cmd, "enhance.ok")); |
117 | 117 |
|
118 | 118 |
|
119 | 119 | $cmd = "$utildir/feature_scoring.+-.pl " |
|
122 | 122 | . " --atg_position $pwm_left " |
123 | 123 | . " > $out_prefix.enhanced.feature.scores"; |
124 | 124 |
|
125 | | - $pipeliner->add_commands(new Command($cmd, "$checkpoints_dir/enhanced_score_features.ok")); |
| 125 | + $pipeliner->add_commands(new Command($cmd, "enhanced_score_features.ok")); |
126 | 126 |
|
127 | 127 |
|
128 | 128 | $cmd = "$utildir/feature_scores_to_ROC.pl $out_prefix.enhanced.feature.scores > $out_prefix.enhanced.feature.scores.roc"; |
129 | | - $pipeliner->add_commands(new Command($cmd, "$checkpoints_dir/enhanced_roc_features.ok")); |
| 129 | + $pipeliner->add_commands(new Command($cmd, "enhanced_roc_features.ok")); |
130 | 130 |
|
131 | 131 | $cmd = "$utildir/plot_ROC.Rscript $out_prefix.enhanced.feature.scores.roc || :"; |
132 | | - $pipeliner->add_commands(new Command($cmd, "$checkpoints_dir/enhanced_rocplot.ok")); |
| 132 | + $pipeliner->add_commands(new Command($cmd, "enhanced_rocplot.ok")); |
133 | 133 |
|
134 | 134 | $cmd = "$utildir/compute_AUC.pl $out_prefix.enhanced.feature.scores.roc"; |
135 | | - $pipeliner->add_commands(new Command($cmd, "$checkpoints_dir/enhanced_aucplot.ok")); |
| 135 | + $pipeliner->add_commands(new Command($cmd, "enhanced_aucplot.ok")); |
136 | 136 |
|
137 | 137 | $cmd = "$utildir/make_seqLogo.Rscript $out_prefix.enhanced.+.pwm || :"; |
138 | | - $pipeliner->add_commands(new Command($cmd, "$checkpoints_dir/seqlogo.enhanced.+.ok")); |
| 138 | + $pipeliner->add_commands(new Command($cmd, "seqlogo.enhanced.+.ok")); |
139 | 139 |
|
140 | 140 |
|
141 | 141 | $pipeliner->run(); |
|
0 commit comments