Merge pull request #5 from belleau/master

move u initialisation outside of the if-else
ArnaudDroitLab · Jun 9, 2016 · 7849252 · 7849252
2 parents 74629a4 + 0e18345
commit 7849252
Showing 1 changed file with 2 additions and 3 deletions.
diff --git a/R/rjmcmcMethod.R b/R/rjmcmcMethod.R
@@ -212,10 +212,10 @@ rjmcmc <- function(startPosForwardReads, startPosReverseReads,
     for (i in 2:nbrIterations) {
         ## List of current values that is going to be passed to sub-functions
         varTilde <- list()
-
+        u <- runif(1)
         if (kValue == 1L) {
             ## CASE : Number of nucleosomes equal to 1
-            u <- runif(1)
+
 
             if (u <= 0.5) {
                 ## Birth move in case k=1
@@ -233,7 +233,6 @@ rjmcmc <- function(startPosForwardReads, startPosReverseReads,
         }  ## end CASE : Number of nucleosomes equal to 1
         else {
             ## CASE : Number of nucleosomes larger than 1
-            u<-runif(1)
 
             if (u <= Dk(kValue, lambda, kMax)) {
                 ## Death move