diff --git a/DESCRIPTION b/DESCRIPTION
index 7f0fbffd..e6ba893a 100644
--- a/DESCRIPTION
+++ b/DESCRIPTION
@@ -1,6 +1,6 @@
 Package: gsDesign2
 Title: Group Sequential Design with Non-Constant Effect
-Version: 1.1.3.3
+Version: 1.1.3.4
 Authors@R: c(
     person("Keaven", "Anderson", email = "keaven_anderson@merck.com", role = c("aut")),
     person("Yilong", "Zhang", email = "elong0527@gmail.com", role = c("aut")),
diff --git a/NAMESPACE b/NAMESPACE
index 40b95269..d5b4bff5 100644
--- a/NAMESPACE
+++ b/NAMESPACE
@@ -50,6 +50,7 @@ export(gs_update_ahr)
 export(ppwe)
 export(pw_info)
 export(s2pwe)
+export(text_summary)
 export(to_integer)
 export(wlr_weight_1)
 export(wlr_weight_fh)
diff --git a/R/gs_design_ahr.R b/R/gs_design_ahr.R
index 42df836d..86f87f3c 100644
--- a/R/gs_design_ahr.R
+++ b/R/gs_design_ahr.R
@@ -406,5 +406,7 @@ gs_design_ahr <- function(
   )
 
   ans <- add_class(ans, if (!binding) "non_binding", "ahr", "gs_design")
+  attr(ans, 'uninteger_is_from') <- "gs_design_ahr"
+
   return(ans)
 }
diff --git a/R/gs_design_rd.R b/R/gs_design_rd.R
index 64226761..3ea80872 100644
--- a/R/gs_design_rd.R
+++ b/R/gs_design_rd.R
@@ -278,5 +278,7 @@ gs_design_rd <- function(p_c = tibble::tibble(stratum = "All", rate = .2),
   )
 
   ans <- add_class(ans, if (!binding) "non_binding", "rd", "gs_design")
+  attr(ans, 'uninteger_is_from') <- "gs_design_rd"
+
   return(ans)
 }
diff --git a/R/gs_design_wlr.R b/R/gs_design_wlr.R
index cc370ee9..18fce32d 100644
--- a/R/gs_design_wlr.R
+++ b/R/gs_design_wlr.R
@@ -340,6 +340,8 @@ gs_design_wlr <- function(
     analysis = analysis)
 
   ans <- add_class(ans, if (!binding) "non_binding", "wlr", "gs_design")
+  attr(ans, 'uninteger_is_from') <- "gs_design_wlr"
+
   return(ans)
 }
 
diff --git a/R/gs_power_ahr.R b/R/gs_power_ahr.R
index b769a57a..364980ec 100644
--- a/R/gs_power_ahr.R
+++ b/R/gs_power_ahr.R
@@ -122,7 +122,7 @@
 #' # Events = max(events, calculated events for targeted analysis_time)
 #' \donttest{
 #' gs_power_ahr(
-#'   analysis_time = c(12, 24, 36), 
+#'   analysis_time = c(12, 24, 36),
 #'   event = c(30, 40, 50), h1_spending = FALSE,
 #'   binding = TRUE,
 #'   upper = gs_spending_bound,
@@ -208,14 +208,14 @@ gs_power_ahr <- function(
     x$info <- x$info * q
   }
 
-  if (h1_spending) { 
+  if (h1_spending) {
       theta1 <- x$theta
       info1 <- x$info
   } else {
       theta1 <- 0
       info1 <- x$info0
   }
-    
+
   # Given the above statistical information, calculate the power ----
   y_h1 <- gs_power_npe(
     theta = x$theta, theta0 = 0, theta1 = theta1,
@@ -281,10 +281,12 @@ gs_power_ahr <- function(
     enroll_rate = enroll_rate, fail_rate = fail_rate,
     event = event, analysis_time = analysis_time,
     info_scale = info_scale,
+    alpha = if (identical(upper, gs_spending_bound)) {upar$total_spend} else {NULL},
     upper = upper, upar = upar,
     lower = lower, lpar = lpar,
     test_lower = test_lower, test_upper = test_upper,
-    ratio = ratio, binding = binding, info_scale = info_scale, r = r, tol = tol
+    ratio = ratio, binding = binding, h1_spending = h1_spending,
+    info_scale = info_scale, r = r, tol = tol
   )
 
   ans <- list(
@@ -296,5 +298,7 @@ gs_power_ahr <- function(
   )
 
   ans <- add_class(ans, if (!binding) "non_binding", "ahr", "gs_design")
+  attr(ans, 'uninteger_is_from') <- "gs_power_ahr"
+
   return(ans)
 }
diff --git a/R/gs_power_rd.R b/R/gs_power_rd.R
index a0f2c9a5..6d16a50e 100644
--- a/R/gs_power_rd.R
+++ b/R/gs_power_rd.R
@@ -349,5 +349,7 @@ gs_power_rd <- function(
   )
 
   ans <- add_class(ans, if (!binding) "non_binding", "rd", "gs_design")
+  attr(ans, 'uninteger_is_from') <- "gs_power_rd"
+
   return(ans)
 }
diff --git a/R/gs_power_wlr.R b/R/gs_power_wlr.R
index 8eb3b6bd..e6c93e32 100644
--- a/R/gs_power_wlr.R
+++ b/R/gs_power_wlr.R
@@ -305,6 +305,8 @@ gs_power_wlr <- function(enroll_rate = define_enroll_rate(duration = c(2, 2, 10)
   )
 
   ans <- add_class(ans, if (!binding) "non_binding", "wlr", "gs_design")
+  attr(ans, 'uninteger_is_from') <- "gs_power_wlr"
+
   return(ans)
 }
 
diff --git a/R/text_summary.R b/R/text_summary.R
new file mode 100644
index 00000000..3a792e7e
--- /dev/null
+++ b/R/text_summary.R
@@ -0,0 +1,227 @@
+#  Copyright (c) 2024 Merck & Co., Inc., Rahway, NJ, USA and its affiliates.
+#  All rights reserved.
+#
+#  This file is part of the gsDesign2 program.
+#
+#  gsDesign2 is free software: you can redistribute it and/or modify
+#  it under the terms of the GNU General Public License as published by
+#  the Free Software Foundation, either version 3 of the License, or
+#  (at your option) any later version.
+#
+#  This program is distributed in the hope that it will be useful,
+#  but WITHOUT ANY WARRANTY; without even the implied warranty of
+#  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+#  GNU General Public License for more details.
+#
+#  You should have received a copy of the GNU General Public License
+#  along with this program.  If not, see <http://www.gnu.org/licenses/>.
+
+#' Generates a textual summary of a group sequential design using the AHR method.
+#' @param x A design object created by [gs_design_ahr()] with or without [to_integer()].
+#' @param information A logical value indicating whether to include statistical information in the textual summary. Default is FALSE.
+#' @param time_unit A character string specifying the time unit used in the design. Options include "days", "weeks", "months" (default), and "years".
+#' @return A character string containing a paragraph that summarizes the design.
+#'
+#' @export
+#' @examples
+#' library(gsDesign)
+#'
+#' # Text summary of a 1-sided design
+#' x <- gs_design_ahr(info_frac = 1:3/3, test_lower = FALSE)
+#' x |> text_summary()
+#' x |> to_integer() |> text_summary()
+#' gs_power_ahr(event = c(10, 20, 30), test_lower = FALSE) |> text_summary()
+#' gs_power_ahr(event = c(10, 20, 30), test_lower = FALSE) |> to_integer() |> text_summary()
+#'
+#' # Text summary of a 2-sided symmetric design
+#' x <- gs_design_ahr(info_frac = 1:3/3,
+#'                    upper = gs_spending_bound, lower = gs_spending_bound,
+#'                    upar = list(sf = sfLDOF, total_spend = 0.025),
+#'                    lpar = list(sf = sfLDOF, total_spend = 0.025),
+#'                    binding = TRUE, h1_spending = FALSE) |> to_integer()
+#' x |> text_summary()
+#'
+#' # Text summary of a asymmetric 2-sided design with beta-spending and non-binding futility bound
+#' x <- gs_design_ahr(info_frac = 1:3/3, alpha = 0.025, beta = 0.1,
+#'                    upper = gs_spending_bound, lower = gs_spending_bound,
+#'                    upar = list(sf = sfLDOF, total_spend = 0.025),
+#'                    lpar = list(sf = sfHSD, total_spend = 0.1, param = -4),
+#'                    binding = FALSE, h1_spending = TRUE) |> to_integer()
+#' x |> text_summary()
+#'
+#' # Text summary of a asymmetric 2-sided design with fixed non-binding futility bound
+#' x <- gs_design_ahr(info_frac = 1:3/3, alpha = 0.025, beta = 0.1,
+#'                    upper = gs_spending_bound, lower = gs_b,
+#'                    upar = list(sf = sfLDOF, total_spend = 0.025),
+#'                    test_upper = c(FALSE, TRUE, TRUE),
+#'                    lpar = c(-1, -Inf, -Inf),
+#'                    test_lower = c(TRUE, FALSE, FALSE),
+#'                    binding = FALSE, h1_spending = TRUE) |> to_integer()
+#' x |> text_summary()
+#'
+#' # If there are >5 pieces of HRs, we provide a brief summary of HR.
+#' gs_design_ahr(
+#'   fail_rate = define_fail_rate(duration = c(rep(3, 5), Inf),
+#'                                hr = c(0.9, 0.8, 0.7, 0.6, 0.5, 0.4),
+#'                                fail_rate = log(2) / 10, dropout_rate = 0.001),
+#'   info_frac = 1:3/3, test_lower = FALSE) |>
+#'   text_summary()
+text_summary <- function(x, information = FALSE, time_unit = "months") {
+
+  n_analysis <- nrow(x$analysis)
+  is_gs_design_ahr <- attributes(x)$uninteger_is_from == "gs_design_ahr"
+  is_gs_power_ahr <- attributes(x)$uninteger_is_from == "gs_power_ahr"
+  is_gs_design_wlr <- attributes(x)$uninteger_is_from == "gs_design_wlr"
+  is_gs_power_wlr <- attributes(x)$uninteger_is_from == "gs_power_wlr"
+  is_gs_design_rd <- attributes(x)$uninteger_is_from == "gs_design_rd"
+  is_gs_power_rd <- attributes(x)$uninteger_is_from == "gs_power_rd"
+
+  # ---------------------------------------- #
+  # Check if it is two-sided design or not
+  # ---------------------------------------- #
+  if ((identical(x$input$lower, gs_b) && (!is.list(x$input$lpar))) || all(!x$input$test_lower)) {
+    if (all(x$input$test_lower == FALSE)) {
+      two_sided <- FALSE
+    } else {
+      two_sided <- ifelse(identical(x$input$lpar, rep(-Inf, n_analysis)), FALSE, TRUE)
+    }
+  } else {
+    two_sided <- TRUE
+  }
+
+  # ---------------------------------------- #
+  # Initialize the output
+  # ---------------------------------------- #
+  out <- NULL
+
+  # ---------------------------------------- #
+  # Add the test type
+  # ---------------------------------------- #
+  ## test_type = 1
+  if (!two_sided) {
+    out <- paste(out, "One-sided group sequential design with ", sep = "")
+  } else {
+  # test_type = 2
+    if (identical(x$input$upper, gs_spending_bound) && identical(x$input$lower, gs_spending_bound) &&
+        identical(x$input$upar, x$input$lpar) && x$input$binding && !x$input$h1_spending) {
+      out <- paste(out, "Symmetric two-sided group sequential design with ", sep = "")
+  # test_type = 3, 4, 5, 6
+    } else {
+      out <- paste(out, "Asymmetric two-sided group sequential design with ", sep = "")
+      if (x$input$binding) {
+        out <- paste(out, "binding futility bound, ", sep = "")
+      } else {
+        out <- paste(out, "non-binding futility bound, ", sep = "")
+      }
+    }
+  }
+
+  # ---------------------------------------- #
+  # Add the number of analyses, sample size, events
+  # ---------------------------------------- #
+  out <- paste(out, n_analysis, " analyses, ", sep = "")
+
+  out <- paste(out,
+               "time-to-event outcome with sample size ",
+               ifelse(is_wholenumber(max(x$analysis$n)), max(x$analysis$n), max(x$analysis$n) |> round(1)),
+               " and ",
+               ifelse(is_wholenumber(max(x$analysis$event)), max(x$analysis$event), max(x$analysis$event) |> round(1)),
+               " events, ",
+               sep = "")
+
+  # ---------------------------------------- #
+  # Add information, power and type I error
+  # ---------------------------------------- #
+  if (information) {
+    out <- paste(out, " total information ", round(x$analysis$info[n_analysis], 2), ", ", sep = "")
+  }
+
+  # ---------------------------------------- #
+  # Add power and type I error
+  # ---------------------------------------- #
+  # if it is a gs_design_ahr object...
+  if (is_gs_design_ahr) {
+    out <- paste(out, 100 * round(x$bound$probability[x$bound$bound == "upper" & x$bound$analysis == n_analysis], 2), " percent power, ", 100 * x$input$alpha, " percent (1-sided) Type I error", sep = "")
+  } else if (is_gs_power_ahr) {
+    out <- paste(out, 100 * x$input$alpha, " percent (1-sided) Type I error", sep = "")
+  }
+
+  # ---------------------------------------- #
+  # Add HR assumption
+  # ---------------------------------------- #
+  if (nrow(x$fail_rate) == 1) {
+    temp <- paste("a hazard ratio of ", round(x$fail_rate$hr, 2), sep = "")
+  } else if (nrow(x$fail_rate) == 2) {
+    temp <- paste("hazard ratio of ",
+                  round(x$fail_rate$hr[1], 2), " during the first ", round(x$fail_rate$duration[1], 2), " ", time_unit,
+                  " and ", round(x$fail_rate$hr[2], 2), " thereafter", sep = "")
+  } else if (nrow(x$fail_rate) <= 5) {
+    temp <- paste(x$fail_rate$hr[1:(nrow(x$fail_rate) - 1)] |> round(2),
+                  c(" during the first ", rep(" during the next ", nrow(x$fail_rate) - 2)),
+                  c(x$fail_rate$duration[1:(nrow(x$fail_rate) - 1)] |> round(2)), time_unit) |>
+      paste(collapse = ", ") |>
+      paste(" and ", x$fail_rate$hr[nrow(x$fail_rate)] |> round(2), " thereafter")
+    temp <- paste("hazard ratio of ", temp, sep = "")
+  } else {
+    temp <- "piecewise hazard ratio"
+  }
+
+  if (is_gs_design_ahr) {
+    out <- paste(out, " to detect ", temp, sep = "")
+  } else if (is_gs_power_ahr) {
+    out_end <- paste(" With ", temp,
+                     ", the power is ",
+                     100 * round(x$bound$probability[x$bound$analysis == n_analysis & x$bound$bound == "upper"], 2),
+                     " percent.", sep = "")
+  }
+  # ---------------------------------------- #
+  # Add enrollment and study duration
+  # ---------------------------------------- #
+  out <- paste(out, ". Enrollment and total study durations are assumed to be ", round(sum(x$enroll_rate$duration), 1),
+               " and ", round(max(x$analysis$time), 1), " ", time_unit, ", respectively",
+               sep = "")
+
+  # ---------------------------------------- #
+  # Add upper bounds derivation
+  # ---------------------------------------- #
+  if (identical(x$input$upper, gs_spending_bound) && identical(x$input$lower, gs_spending_bound) &&
+      identical(x$input$upar, x$input$lpar) && x$input$binding && !x$input$h1_spending) {
+    out <- paste(out, ". Bounds derived using a ", sep = "")
+  } else {
+    out <- paste(out, ". Efficacy bounds derived using a", sep = "")
+  }
+
+  analysis_seq <- c(paste("IA", 1:(n_analysis - 1), sep = ""), "FA")
+  upper_text <- x$input$upar$sf(alpha = x$input$upar$total_spend, t = x$analysis$info_frac, param = x$input$upar$param)
+  upper_tested <- if (!all(x$input$test_upper)) {
+    paste(", tested at", paste("tested at", paste(analysis_seq[x$input$test_upper], collapse = ", ")))
+  }
+  out <- paste(out, " ", summary(upper_text), upper_tested, ".", sep = "")
+
+  # ---------------------------------------- #
+  # Add lower bounds derivation
+  # ---------------------------------------- #
+  if (any(x$bound$bound == "lower") &&
+      !(identical(x$input$upper, gs_spending_bound) && identical(x$input$lower, gs_spending_bound) &&
+        identical(x$input$upar, x$input$lpar) && x$input$binding && !x$input$h1_spending)) {
+    lower_tested <- if (!all(x$input$test_lower)) {
+      paste(", tested at", paste("tested at", paste(analysis_seq[x$input$test_lower], collapse = ", ")))
+    }
+
+    if (identical(x$input$lower, gs_spending_bound)) {
+      lower_text <- x$input$lpar$sf(alpha = x$input$lpar$total_spend, t = x$analysis$info_frac, param = x$input$lpar$param)
+      out <- paste(out, " Futility bounds derived using a ", summary(lower_text), lower_tested, ".", sep = "")
+    } else if (identical(x$input$lower, gs_b)) {
+      out <- paste(out, " Futility bounds is fixed as ", paste0(x$input$lpar, collapse = ", ") , lower_tested, ".", sep = "")
+    }
+  }
+
+  # ---------------------------------------- #
+  # Add power for gs_power_ahr object
+  # ---------------------------------------- #
+  if (is_gs_power_ahr) {
+    out <- paste(out, out_end, sep = "")
+  }
+
+  return(out)
+}
diff --git a/R/to_integer.R b/R/to_integer.R
index 2582406f..20d3f42e 100644
--- a/R/to_integer.R
+++ b/R/to_integer.R
@@ -570,5 +570,20 @@ to_integer.gs_design <- function(x, round_up_final = TRUE, ratio = x$input$ratio
   x_new$analysis$n <- round(x_new$analysis$n)
   if (!is_rd) x_new$analysis$event <- round(x_new$analysis$event)
 
+  # Add attributes to x_new to identify whether it is a gs_design_ahr orbject or gs_power_ahr object
+  if ("analysis_time" %in% names(x$input) && "info_frac" %in% names(x$input) && "ahr" %in% class(x)) {
+    attr(x_new, 'uninteger_is_from') <- "gs_design_ahr"
+  } else if ("analysis_time" %in% names(x$input) && "event" %in% names(x$input) && "ahr" %in% class(x)) {
+    attr(x_new, 'uninteger_is_from') <- "gs_power_ahr"
+  } else if ("analysis_time" %in% names(x$input) && "info_frac" %in% names(x$input) && "wlr" %in% class(x)) {
+    attr(x_new, 'uninteger_is_from') <- "gs_design_wlr"
+  } else if ("analysis_time" %in% names(x$input) && "event" %in% names(x$input) && "wlr" %in% class(x)) {
+    attr(x_new, 'uninteger_is_from') <- "gs_power_wlr"
+  } else if (!("n" %in% names(x$input)) && "rd" %in% class(x)) {
+    attr(x_new, 'uninteger_is_from') <- "gs_design_rd"
+  } else if ("n" %in% names(x$input) && "rd" %in% class(x)) {
+    attr(x_new, 'uninteger_is_from') <- "gs_power_rd"
+  }
+
   return(x_new)
 }
diff --git a/_pkgdown.yml b/_pkgdown.yml
index d7f2e830..3f8f70d9 100644
--- a/_pkgdown.yml
+++ b/_pkgdown.yml
@@ -90,6 +90,7 @@ reference:
   contents:
     - summary.fixed_design
     - summary.gs_design
+    - text_summary
     - as_gt
     - as_gt.fixed_design
     - as_gt.gs_design
diff --git a/man/gs_power_ahr.Rd b/man/gs_power_ahr.Rd
index 29b6ae87..aa85bc3e 100644
--- a/man/gs_power_ahr.Rd
+++ b/man/gs_power_ahr.Rd
@@ -207,8 +207,8 @@ gs_power_ahr(
 \donttest{
 gs_power_ahr(
   analysis_time = c(12, 24, 36),
-  event = c(30, 40, 50),
-  binding = TRUE, h1_spending = FALSE,
+  event = c(30, 40, 50), h1_spending = FALSE,
+  binding = TRUE,
   upper = gs_spending_bound,
   upar = list(sf = gsDesign::sfLDOF, total_spend = 0.025),
   lower = gs_spending_bound,
diff --git a/man/text_summary.Rd b/man/text_summary.Rd
new file mode 100644
index 00000000..c4e41195
--- /dev/null
+++ b/man/text_summary.Rd
@@ -0,0 +1,65 @@
+% Generated by roxygen2: do not edit by hand
+% Please edit documentation in R/text_summary.R
+\name{text_summary}
+\alias{text_summary}
+\title{Generates a textual summary of a group sequential design using the AHR method.}
+\usage{
+text_summary(x, information = FALSE, time_unit = "months")
+}
+\arguments{
+\item{x}{A design object created by \code{\link[=gs_design_ahr]{gs_design_ahr()}} with or without \code{\link[=to_integer]{to_integer()}}.}
+
+\item{information}{A logical value indicating whether to include statistical information in the textual summary. Default is FALSE.}
+
+\item{time_unit}{A character string specifying the time unit used in the design. Options include "days", "weeks", "months" (default), and "years".}
+}
+\value{
+A character string containing a paragraph that summarizes the design.
+}
+\description{
+Generates a textual summary of a group sequential design using the AHR method.
+}
+\examples{
+library(gsDesign)
+
+# Text summary of a 1-sided design
+x <- gs_design_ahr(info_frac = 1:3/3, test_lower = FALSE)
+x |> text_summary()
+x |> to_integer() |> text_summary()
+gs_power_ahr(event = c(10, 20, 30), test_lower = FALSE) |> text_summary()
+gs_power_ahr(event = c(10, 20, 30), test_lower = FALSE) |> to_integer() |> text_summary()
+
+# Text summary of a 2-sided symmetric design
+x <- gs_design_ahr(info_frac = 1:3/3,
+                   upper = gs_spending_bound, lower = gs_spending_bound,
+                   upar = list(sf = sfLDOF, total_spend = 0.025),
+                   lpar = list(sf = sfLDOF, total_spend = 0.025),
+                   binding = TRUE, h1_spending = FALSE) |> to_integer()
+x |> text_summary()
+
+# Text summary of a asymmetric 2-sided design with beta-spending and non-binding futility bound
+x <- gs_design_ahr(info_frac = 1:3/3, alpha = 0.025, beta = 0.1,
+                   upper = gs_spending_bound, lower = gs_spending_bound,
+                   upar = list(sf = sfLDOF, total_spend = 0.025),
+                   lpar = list(sf = sfHSD, total_spend = 0.1, param = -4),
+                   binding = FALSE, h1_spending = TRUE) |> to_integer()
+x |> text_summary()
+
+# Text summary of a asymmetric 2-sided design with fixed non-binding futility bound
+x <- gs_design_ahr(info_frac = 1:3/3, alpha = 0.025, beta = 0.1,
+                   upper = gs_spending_bound, lower = gs_b,
+                   upar = list(sf = sfLDOF, total_spend = 0.025),
+                   test_upper = c(FALSE, TRUE, TRUE),
+                   lpar = c(-1, -Inf, -Inf),
+                   test_lower = c(TRUE, FALSE, FALSE),
+                   binding = FALSE, h1_spending = TRUE) |> to_integer()
+x |> text_summary()
+
+# If there are >5 pieces of HRs, we provide a brief summary of HR.
+gs_design_ahr(
+  fail_rate = define_fail_rate(duration = c(rep(3, 5), Inf),
+                               hr = c(0.9, 0.8, 0.7, 0.6, 0.5, 0.4),
+                               fail_rate = log(2) / 10, dropout_rate = 0.001),
+  info_frac = 1:3/3, test_lower = FALSE) |>
+  text_summary()
+}
diff --git a/tests/testthat/_snaps/independent_as_gt.md b/tests/testthat/_snaps/independent_as_gt.md
index 232077ae..f2898cfa 100644
--- a/tests/testthat/_snaps/independent_as_gt.md
+++ b/tests/testthat/_snaps/independent_as_gt.md
@@ -1,6 +1,6 @@
 # Snapshot test for fixed_design summary as_gt
 
-    \begin{table}[!t]
+    \begin{table}[t]
     \caption*{
     {\large Fixed Design under AHR Method\textsuperscript{\textit{1}}}
     } 
@@ -19,7 +19,7 @@
 
 # Snapshot test for fixed_design summary as_gt with custom title and footnote
 
-    \begin{table}[!t]
+    \begin{table}[t]
     \caption*{
     {\large Custom Title\textsuperscript{\textit{1}}}
     } 
@@ -38,7 +38,7 @@
 
 # Snapshot test for fixed_design_fh summary as_gt
 
-    \begin{table}[!t]
+    \begin{table}[t]
     \caption*{
     {\large Fixed Design under Fleming-Harrington Method\textsuperscript{\textit{1}}}
     } 
@@ -57,7 +57,7 @@
 
 # Snapshot test for gs_design_ahr summary as_gt
 
-    \begin{table}[!t]
+    \begin{table}[t]
     \caption*{
     {\large Bound summary for AHR design} \\ 
     {\small AHR approximations of \textasciitilde{}HR at bound}
@@ -82,7 +82,7 @@
 
 # Snapshot test for gs_power_ahr summary as_gt
 
-    \begin{table}[!t]
+    \begin{table}[t]
     \caption*{
     {\large Bound summary for AHR design} \\ 
     {\small AHR approximations of \textasciitilde{}HR at bound}
@@ -118,7 +118,7 @@
 
 # Snapshot test for gs_design_wlr summary as_gt
 
-    \begin{table}[!t]
+    \begin{table}[t]
     \caption*{
     {\large Bound summary for WLR design} \\ 
     {\small WLR approximation of \textasciitilde{}wHR at bound}
@@ -144,7 +144,7 @@
 
 # Snapshot test for gs_power_wlr summary as_gt
 
-    \begin{table}[!t]
+    \begin{table}[t]
     \caption*{
     {\large Bound summary for WLR design\textsuperscript{\textit{1}}} \\ 
     {\small WLR approximation of \textasciitilde{}wHR at bound}
@@ -182,7 +182,7 @@
 
 # Snapshot test for gs_power_combo summary as_gt
 
-    \begin{table}[!t]
+    \begin{table}[t]
     \caption*{
     {\large Bound summary for MaxCombo design} \\ 
     {\small MaxCombo approximation}
@@ -218,7 +218,7 @@
 
 # Snapshot test for gs_design_rd summary as_gt
 
-    \begin{table}[!t]
+    \begin{table}[t]
     \caption*{
     {\large Bound summary of Binary Endpoint} \\ 
     {\small measured by risk difference}
@@ -242,7 +242,7 @@
 
 # Snapshot test for gs_power_rd summary as_gt
 
-    \begin{table}[!t]
+    \begin{table}[t]
     \caption*{
     {\large Bound summary of Binary Endpoint} \\ 
     {\small measured by risk difference}
@@ -276,7 +276,7 @@
 
 # Snapshot test for gs_power_wlr summary as_gt with custom title and subtitle
 
-    \begin{table}[!t]
+    \begin{table}[t]
     \caption*{
     {\large Bound Summary} \\ 
     {\small from gs\_power\_wlr}
@@ -313,7 +313,7 @@
 
 # Snapshot test for gs_power_wlr summary as_gt with colname_spanner and colname_spannersub
 
-    \begin{table}[!t]
+    \begin{table}[t]
     \caption*{
     {\large Bound summary for WLR design} \\ 
     {\small WLR approximation of \textasciitilde{}wHR at bound}
@@ -350,7 +350,7 @@
 
 # Snapshot test for gs_power_wlr summary as_gt with custom footnotes
 
-    \begin{table}[!t]
+    \begin{table}[t]
     \caption*{
     {\large Bound summary for WLR design\textsuperscript{\textit{1}}} \\ 
     {\small WLR approximation of \textasciitilde{}wHR at bound}
@@ -388,7 +388,7 @@
 
 # Snapshot test for gs_power_wlr summary as_gt with display_bound
 
-    \begin{table}[!t]
+    \begin{table}[t]
     \caption*{
     {\large Bound summary for WLR design} \\ 
     {\small WLR approximation of \textasciitilde{}wHR at bound}
@@ -422,7 +422,7 @@
 
 # Snapshot test for gs_power_wlr summary as_gt with display_columns
 
-    \begin{table}[!t]
+    \begin{table}[t]
     \caption*{
     {\large Bound summary for WLR design} \\ 
     {\small WLR approximation of \textasciitilde{}wHR at bound}