RobinHankin / hyper2

https://robinhankin.github.io/hyper2/
5 stars 3 forks source link

zero-weighted brackets in hyper3 objects #195

Open RobinHankin opened 1 year ago

RobinHankin commented 1 year ago
suppressMessages(library("hyper2"))
(H <-  hyper3(list(c(a=1,b=2),c(a=0,b=0)),powers=c(3,4)))
#> log( ()^4 * (a=1, b=2)^3)
loglik(c(a=0.2,b=0.8),H)
#> [1] -Inf

Created on 2023-06-26 with reprex v2.0.2

It is not clear to me what the correct behaviour should be. Returning -Inf is at least consistent but is perhaps undesirable.