Added cumincglm quasi glm class option

gmelloni · Jan 11, 2024 · 6e60465 · 6e60465
1 parent 4ee8a99
commit 6e60465
Show file tree

Hide file tree

Showing 2 changed files with 12 additions and 1 deletion.
diff --git a/R/toyexample.R b/R/toyexample.R
@@ -2,6 +2,17 @@
 # Testing multiple knots #
 ##########################
 
+# library(eventglm)
+# library(survival)
+# library(rms)
+# data(cancer)
+# cancer$sex <- ifelse(cancer$sex==2 , 0 , cancer$sex)
+# myformula <- Surv(time, status) ~ age*sex
+# cuminc <- cumincglm(myformula, data=cancer, time = 365, link = "identity", model.censoring="coxph")
+# cumincINT <- intEST( var2values = 40:80
+#                      , model = cuminc , data = cancer , var1 ="sex", var2="age"
+#                      , ci=TRUE , conf = 0.95 , ci.method = "delta" , R=100)
+# plotINT(cumincINT)
 # library(rms)
 # library(mlbench)
 # data(PimaIndiansDiabetes)

diff --git a/R/wrapper_int.R b/R/wrapper_int.R
@@ -66,7 +66,7 @@ intEST <- function(var2values , model , data , var1 , var2
     }
   } else if(any(c("Glm" , "glm") %in% modelClass)){
     family <- model$family$family
-    if(family == "gaussian"){
+    if(family %in% c("gaussian","quasi")){
       if(isRCS){
         argg[[1]] <- as.name("rcsLIN")
       } else {