bootbayes now supports NaN value for PROB input argument

gnu-octave · Jun 16, 2024 · 2ec0144 · 2ec0144
1 parent 6970bb2
commit 2ec0144
Showing 1 changed file with 14 additions and 11 deletions.
diff --git a/inst/bootbayes.m b/inst/bootbayes.m
@@ -386,17 +386,20 @@
   bias = mean (bootstat, 2) - original;
 
   % Compute credible intervals
-  if (prior > 0)
-    ci = credint (bootstat, prob);
-  else
-    stdnorminv = @(p) sqrt (2) * erfinv (2 * p - 1);
-    switch nprob
-      case 1
-        z = stdnorminv (1 - (1 - prob) / 2);
-        ci = original + std (bootstat, 1 , 2) * z * [-1, 1];
-      case 2
-        z = stdnorminv (prob);
-        ci = original + std (bootstat, 1 , 2) * z;
+  ci = nan (p, 2);
+  if (any (~ isnan (prob)))
+    if (prior > 0)
+      ci = credint (bootstat, prob);
+    else
+      stdnorminv = @(p) sqrt (2) * erfinv (2 * p - 1);
+      switch nprob
+        case 1
+          z = stdnorminv (1 - (1 - prob) / 2);
+          ci = original + std (bootstat, 1 , 2) * z * [-1, 1];
+        case 2
+          z = stdnorminv (prob);
+          ci = original + std (bootstat, 1 , 2) * z;
+      end
     end
   end