Properly deal with non-zero 'mumiss' value in fsn().

R/fsn.r

@@ -92,7 +92,7 @@ fsn <- function(x, vi, sei, subset, data, type, alpha=.05, target,
          stop(mstyle$stop("Method only applicable to models without moderators."))
 
       if (!missing(type) && type != "General")
-         warning(mstyle$warning("Setting type='General' when using fsn() on a model object. "), call.=FALSE)
+         warning(mstyle$warning("Setting type='General' when using fsn() on a model object."), call.=FALSE)
 
       type <- "General"
 
@@ -365,6 +365,11 @@ fsn <- function(x, vi, sei, subset, data, type, alpha=.05, target,
       tau2 <- tmp$tau2
       pval <- tmp$pval
 
+      if (mumiss != 0 && sign(est) == sign(mumiss)) {
+         mumiss <- -mumiss
+         message(mstyle$message("Flipped the sign of 'mumiss'."))
+      }
+
       if (is.null(target)) {
 
          if (pval >= alpha) {

R/misc.func.hidden.fsn.r

@@ -38,6 +38,9 @@
          tau2.fsn <- tmp$tau2
          pval.fsn <- tmp$pval
 
+         if (mumiss != 0 && sign(est.fsn) == sign(mumiss))
+            pval.fsn <- 1
+
       } else {
 
          k <- length(yi)
@@ -65,8 +68,12 @@
             est.fsn  <- mean(yinew)
             zval.new <- (k + fsnum) * est.fsn / sqrt(sum(vi + tau2.fsn) + fsnum * (vt + tau2.fsn))
          }
+
          pval.fsn <- 2*pnorm(abs(zval.new), lower.tail=FALSE)
 
+         if (mumiss != 0 && sign(est.fsn * mumiss) == 1)
+            pval.fsn <- 1
+
       }
 
       if (newest) {

docs/pkgdown.yml

@@ -3,7 +3,7 @@ pkgdown: 2.0.7
 pkgdown_sha: ~
 articles:
   diagram: pkgdown/diagram.html
-last_built: 2023-09-10T15:37Z
+last_built: 2023-09-10T16:22Z
 urls:
   reference: https://wviechtb.github.io/metafor/reference
   article: https://wviechtb.github.io/metafor/articles