Merge pull request #183 from keaven/180-allow-multiple-alpha-levels-to-gsboundsummary-table

180 allow multiple alpha levels to `gsBoundSummary` table

Author: nanxstats

R/gsCP.R

@@ -228,8 +228,11 @@ gsCP <- function(x, theta = NULL, i = 1, zi = 0, r = 18) {
 #  } else if (test.type > 1 && zi < x$lower$bound[i]) {
 #    stop("gsCP must have x$lower$bound[i]<=zi<=x$upper$bound[i]")
 #  }
-
-  if (is.null(theta)) theta <- c(zi/sqrt(x$n.I[i]), 0, x$delta)
+# 20250208 KA added check for that x$delta is there, if needed
+# This is needed when input is generated by gsProbability
+  
+# Default is to use estimated theta based on interim z, theta = 0 (conditional error) and delta or, if none, last theta
+  if (is.null(theta)) theta <- c(zi/sqrt(x$n.I[i]), 0, ifelse(!is.null(x$delta), x$delta, x$theta[length(x$theta)]))
 
   knew <- x$k - i
   Inew <- x$n.I[(i + 1):x$k] - x$n.I[i]
@@ -374,7 +377,7 @@ gsPI <- function(x, i = 1, zi = 0, j = 2, level = .95, theta = c(0, 3), wgts = c
 # gsBoundCP function [sinew] ----
 gsBoundCP <- function(x, theta = "thetahat", r = 18) {
   if (!(inherits(x, "gsProbability") || inherits(x, "gsDesign"))) {
-    stop("gsPI: class(x) must be gsProbability or gsDesign")
+    stop("gsBoundCP: class(x) must be gsProbability or gsDesign")
   }
   checkScalar(r, "integer", c(1, 70))
   if (!is.character(theta)) {