Merge pull request #18 from cyk0315/master

Fix: error in anova

DESCRIPTION

@@ -1,7 +1,7 @@
 Package: jstable
 Title: Create Tables from Different Types of Regression
-Version: 1.1.5
-Date: 2024-02-02
+Version: 1.1.6
+Date: 2024-02-15
 Authors@R: c(person("Jinseob", "Kim", email = "[email protected]", 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,7 @@
+# jstable 1.1.6
+
+* Bugfix `TableSubgroupGLM`: thanks for `weisx2022`
+
 # jstable 1.1.5
 
 * Bugfix `TableSubgroupCox`: thanks for `ciciing`

R/forestglm.R

@@ -136,9 +136,13 @@ TableSubgroupGLM <- function(formula, var_subgroup = NULL, var_cov = NULL, data,
         coefficients()
       pv_int <- round(pvs_int[nrow(pvs_int), ncol(pvs_int)], decimal.pvalue)
       # if (!is.null(xlev) & length(xlev[[1]]) != 2) stop("Categorical independent variable must have 2 levels.")
-      model.int <- survey::svyglm(as.formula(gsub(xlabel, paste(xlabel, "*", var_subgroup, sep = ""), deparse(formula))), design = data, family = family.svyglm)
-
-
+      data.design<-data
+      if(family=='binomial'){
+        model.int <- survey::svyglm(as.formula(gsub(xlabel, paste(xlabel, "*", var_subgroup, sep = ""), deparse(formula))), design = data.design, family = quasibinomial())
+      }
+      else{
+        model.int <- survey::svyglm(as.formula(gsub(xlabel, paste(xlabel, "*", var_subgroup, sep = ""), deparse(formula))), design = data.design, family = gaussian())
+      }
       if (sum(grepl(":", names(coef(model.int)))) > 1) {
         pv_anova <- anova(model.int, method = "Wald")
         pv_int <- round(pv_anova[[length(pv_anova)]][[7]], decimal.pvalue)