Fix: TableSubgroupMultiCox p-value extraction for survey Cox models

- Fixed p-value column index issue in forestcox.R
- Survey Cox models have 6 columns (including 'robust se'), regular Cox models have 5
- Now using ncol(coef_matrix) to get the last column (always p-value) instead of hardcoded index 5
- Resolves issue where p-values > 1 were returned for survey Cox models with 3+ category variables
- Version bump to 1.3.18

Author: jinseob2kim

.gitignore

@@ -1,4 +1,4 @@
-inst/doc
+/docs
 .Rproj.user
 .Rhistory
 .RData

DESCRIPTION

@@ -1,7 +1,7 @@
 Package: jstable
 Title: Create Tables from Different Types of Regression
-Version: 1.3.17
-Date: 2025-08-29
+Version: 1.3.18
+Date: 2025-09-07
 Authors@R: c(person("Jinseob", "Kim", email = "jinseob2kim@gmail.com", role = c("aut", "cre"),  comment = c(ORCID = "0000-0002-9403-605X")),
              person("Zarathu", role = c("cph", "fnd")),
              person("Yoonkyoung","Jeon", role = c("aut")),

NEWS.md

@@ -1,3 +1,6 @@
+# jstable 1.3.18
+* Fix: `TableSubgroupMultiCox` now correctly extracts p-values from survey Cox models (svycoxph) by using the last column index instead of hardcoded column 5
+
 # jstable 1.3.17
 * Fix: `cox2.display` now properly handles variables with no variation in complete cases during univariate analysis
 * Fix: `cox2.display` correctly filters out NA p-values when using `pcut.univariate` option

R/forestcox.R

@@ -490,7 +490,9 @@ TableSubgroupCox <- function(formula, var_subgroup = NULL, var_cov = NULL, data,
             lapply(function(x) {
               tryCatch(
                 {
-                  round(summary(x)$coefficients[1:ncoef, 5], decimal.pvalue)
+                  # P-value is always in the last column
+                  coef_matrix <- summary(x)$coefficients
+                  round(coef_matrix[1:ncoef, ncol(coef_matrix)], decimal.pvalue)
                 },
                 error = function(e) {
                   return(rep(NA, length(xlev[[1]]) - 1))