bug fix

rsetienne · Jul 16, 2024 · 2918b5a · 2918b5a
1 parent 9077b3e
commit 2918b5a
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/R/DAISIE_loglik_high_lambda.R b/R/DAISIE_loglik_high_lambda.R
@@ -4,14 +4,14 @@ DAISIE_loglik_high_lambda <- function(pars1, brts, stac) {
      brts <- brts[-lbrts]
      lbrts <- length(brts)
    }
+   gam <- pars1[4]
    if (stac == 0) {
      out <- -gam * brts[1]
      return(out)
    }
    if (stac == 2 | stac == 6) {
      N <- lbrts - 1 - (stac == 6)
      mu <- pars1[2]
-     gam <- pars1[4]
      brtsdiff <- brts - c(brts[2:lbrts], 0)
      out <- -gam * brtsdiff[1] +
        log(gam) * (stac == 2) +