From b6b8acc5faa5c907e64c87ad26ce056063ab92c5 Mon Sep 17 00:00:00 2001
From: Ben Bolker <bbolker@gmail.com>
Date: Sat, 18 May 2024 20:44:41 -0400
Subject: [PATCH] fix qKumIW bug

---
 R/dKumIW.R | 6 ++++--
 1 file changed, 4 insertions(+), 2 deletions(-)

diff --git a/R/dKumIW.R b/R/dKumIW.R
index 76a93d36..586dbda2 100644
--- a/R/dKumIW.R
+++ b/R/dKumIW.R
@@ -93,8 +93,10 @@ qKumIW <- function(p, mu, sigma, nu, lower.tail=TRUE, log.p=FALSE){
   else p <- 1 - p
   if (any(p < 0) | any(p > 1)) 
     stop(paste("p must be between 0 and 1", "\n", ""))
-  
-  q <- (-sigma / (log(1 - (1 - p)^(1/nu))))
+
+  ## FIXME: might be able to improve precision/stability via
+  ## logspace computation, e.g. see ?statnet.common::logspace.utils
+  q <- (-sigma / (log(1 - (1 - p)^(1/nu))))^(1/mu)
   q
 }
 #' @importFrom stats runif