사용하신 데이터를 올려봐주시겠어요?

네 제가 살펴보니 문제가 세가지 있었습니다.

1) 먼저 제가 mycph 함수를 만들때 formula가 윚의 경우와 같이 길어질 것을 생각못해 오류가 납니다. 그 부분을 고친 함수를 mycph2로 다시 만들었습니다. 

2) mysph2함수를 호출할 때 데이터에 TS변수를 새로 만들어야 합니다. 

3) formula에 DEATH_F_DATE가 있는데 데이터에는 DAETH_F_DATE 열이 있습니다.

위의 세가지 문제를 해결하니 정상적으로 작동됩니다

library(moonBook)

library(survival)

mycph2=function (formula, data, digits = 2) 

{

     temp=paste0(deparse(formula),collapse="")

     call = paste(temp, ", ", "data= ", substitute(data), 

                  sep = "")

     cat("\n mycph : perform coxph of individual expecting variables\n")

     cat("\n Call:", call, "\n\n")

     f = formula

     myt = terms(f, data = data)

     y = as.character(f[[2]])

     if (class(data[[y]]) != "Surv") {

          cat(y, "is not an object of class Surv")

          return(invisible())

     }

     myvar = attr(myt, "term.labels")

     count = length(myvar)

     var <- HR <- lcl <- ucl <- p.value <- c()

     for (i in 1:count) {

          s = paste(y, myvar[i], sep = "~")

          suppressWarnings(out <- summary(survival::coxph(as.formula(s), 

                                                          data)))

          if (any(is.infinite(out$conf.int))) {

               cat(dimnames(out$conf.int)[[1]], " was excluded : infinite\n")

               next

          }

          if (any(is.nan(out$coef))) {

               cat(dimnames(out$conf.int)[[1]], " was excluded : NaN\n")

               next

          }

          if (any(is.na(out$coef))) {

               cat(dimnames(out$conf.int)[[1]], " was excluded : NA\n")

               next

          }

          var = c(var, dimnames(out$conf.int)[[1]])

          HR = c(HR, out$coef[, 2])

          lcl = c(lcl, out$conf.int[, 3])

          ucl = c(ucl, out$conf.int[, 4])

          p.value = c(p.value, out$coef[, 5])

     }

     if (length(HR) < 1) 

          return(invisible())

     result = round(data.frame(HR, lcl, ucl), digits)

     rownames(result) = var

     result = cbind(result, round(p.value, max(3, digits)))

     colnames(result)[4] = "p"

     result

}