Fixed bugs in the Gibbs sampler. Fixed a bug which will set probF to 0 if the protoco...

[rsem.git] / rsem-calculate-expression

diff --git a/rsem-calculate-expression b/rsem-calculate-expression
index 0f22282dfc8f0b4ea259d057b7aca3ebdc2528ce..7826ec72680bc551f7c4a167d53e2d9d0b830905 100755 (executable)
--- a/rsem-calculate-expression
+++ b/rsem-calculate-expression
@@ -173,7 +173,7 @@ if (!$is_sam && !$is_bam) {
     if ($read_type == 2 || $read_type == 3) { $command .= " -I $minL -X $maxL"; }
     
     if ($strand_specific || $probF == 1.0) { $command .= " --norc"; }
-    elsif ($probF = 0.0) { $command .= " --nofw"; }
+    elsif ($probF == 0.0) { $command .= " --nofw"; }
 
     $command .= " -p $nThreads -a -m $maxHits -S";
     if ($quiet) { $command .= " --quiet"; }