Fix the condition to compute the loglikelihood

Although the pmf is defined in this interval the tweedie::dtweedie only works for power greater than or equal to one.
jreduardo · Apr 9, 2019 · fdc9aa8 · fdc9aa8
1 parent e237367
commit fdc9aa8
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/R/utils.R b/R/utils.R
@@ -22,7 +22,7 @@
 # Compute the loglikelihood for Poisson-Tweedie models
 `_ptw_loglik` <- function(y, mu, omega, power, npts = 100L) {
   # Conditions
-  if (!(omega >= 0 & (power >= 1 | power <= 0))) {
+  if (!(omega >= 0 & (power >= 1))) {
     ll <- NA
   } else {
     lli <- vapply(seq(mu), function(i)