R :Fitting survival trees with time-varying covariates in RandomForestSRC2019 Community Moderator ElectionFitting a fully parametric proportional hazard model with time-varying covariates in RAdding time varying covariates to survival data using 'tmerge' in 'survival' packageSurvival analysis on left truncated data with ipredbagg or pecCforestCox PH modeling and prediction with time varying dependencyHow to find the accuracy for Survival Analysis using RandomForestSRC package in RInterval censored data: Cox proportional hazard and surival difference in RHow to calculate median survival in randomForestSRCSurvival analysis with time varying covariates using coxph in RDeviance residuals diagnostics on a Cox model with 13 independent covariates - any way I can get plots for each and every covariate?time-varying covariates for Survival analysis

Should we release the security issues we found in our product as CVE or we can just update those on weekly release notes?

An Accountant Seeks the Help of a Mathematician

What has been your most complicated TikZ drawing?

Employee lack of ownership

Can elves maintain concentration in a trance?

When do we add an hyphen (-) to a complex adjective word?

Why must traveling waves have the same amplitude to form a standing wave?

Why are the outputs of printf and std::cout different

Does splitting a potentially monolithic application into several smaller ones help prevent bugs?

Sword in the Stone story where the sword was held in place by electromagnets

PTIJ: Who should pay for Uber rides: the child or the parent?

Current sense amp + op-amp buffer + ADC: Measuring down to 0 with single supply

Can anyone tell me why this program fails?

Latest web browser compatible with Windows 98

Why using two cd commands in bash script does not execute the second command

Could the Saturn V actually have launched astronauts around Venus?

Distribution of Maximum Likelihood Estimator

What are some nice/clever ways to introduce the tonic's dominant seventh chord?

Brexit - No Deal Rejection

Does this AnyDice function accurately calculate the number of ogres you make unconcious with three 4th-level castings of Sleep?

2D counterpart of std::array in C++17

Possible Leak In Concrete

How to generate globally unique ids for different tables of the same database?

Why do Australian milk farmers need to protest supermarkets' milk price?



R :Fitting survival trees with time-varying covariates in RandomForestSRC



2019 Community Moderator ElectionFitting a fully parametric proportional hazard model with time-varying covariates in RAdding time varying covariates to survival data using 'tmerge' in 'survival' packageSurvival analysis on left truncated data with ipredbagg or pecCforestCox PH modeling and prediction with time varying dependencyHow to find the accuracy for Survival Analysis using RandomForestSRC package in RInterval censored data: Cox proportional hazard and surival difference in RHow to calculate median survival in randomForestSRCSurvival analysis with time varying covariates using coxph in RDeviance residuals diagnostics on a Cox model with 13 independent covariates - any way I can get plots for each and every covariate?time-varying covariates for Survival analysis










1















In the package LTRCtrees one can fit a decision tree to the special format of the Surv function :Surv(time, time2, event) as per below example



set.seed(0)
library(survival)
library(LTRCtrees)
## Create the start-stop-event triplet needed for coxph and LTRC trees
first <- with(pbcseq, c(TRUE, diff(id) !=0)) #first id for each subject
last <- c(first[-1], TRUE) #last id
time1 <- with(pbcseq, ifelse(first, 0, day))
time2 <- with(pbcseq, ifelse(last, futime, c(day[-1], 0)))
event <- with(pbcseq, ifelse(last, status, 0))
event <- 1*(event==2)

pbcseq$time1 <- time1
pbcseq$time2 <- time2
pbcseq$event <- event

## Fit the Cox model and LTRC trees with time-varying covariates
fit.cox <- coxph(Surv(time1, time2, event) ~ age + sex + log(bili), pbcseq)
LTRCIT.fit <- LTRCIT(Surv(time1, time2, event) ~ age + sex + log(bili), pbcseq)
LTRCART.fit <- LTRCART(Surv(time1, time2, event) ~ age + sex + log(bili), pbcseq)


Is it possible to use the same function on a Random Forest using the RandomForestSRC library



library(randomForestSRC)
RF.fit <- rfsrc(Surv(time1, time2, event) ~ age + sex + log(bili), data=pbcseq, nsplit = 3, ntree = 100, importance = TRUE)


Which generates an error:



Error in parseFormula(formula, data, ytry) : 
Survival formula incorrectly specified.









share|improve this question


























    1















    In the package LTRCtrees one can fit a decision tree to the special format of the Surv function :Surv(time, time2, event) as per below example



    set.seed(0)
    library(survival)
    library(LTRCtrees)
    ## Create the start-stop-event triplet needed for coxph and LTRC trees
    first <- with(pbcseq, c(TRUE, diff(id) !=0)) #first id for each subject
    last <- c(first[-1], TRUE) #last id
    time1 <- with(pbcseq, ifelse(first, 0, day))
    time2 <- with(pbcseq, ifelse(last, futime, c(day[-1], 0)))
    event <- with(pbcseq, ifelse(last, status, 0))
    event <- 1*(event==2)

    pbcseq$time1 <- time1
    pbcseq$time2 <- time2
    pbcseq$event <- event

    ## Fit the Cox model and LTRC trees with time-varying covariates
    fit.cox <- coxph(Surv(time1, time2, event) ~ age + sex + log(bili), pbcseq)
    LTRCIT.fit <- LTRCIT(Surv(time1, time2, event) ~ age + sex + log(bili), pbcseq)
    LTRCART.fit <- LTRCART(Surv(time1, time2, event) ~ age + sex + log(bili), pbcseq)


    Is it possible to use the same function on a Random Forest using the RandomForestSRC library



    library(randomForestSRC)
    RF.fit <- rfsrc(Surv(time1, time2, event) ~ age + sex + log(bili), data=pbcseq, nsplit = 3, ntree = 100, importance = TRUE)


    Which generates an error:



    Error in parseFormula(formula, data, ytry) : 
    Survival formula incorrectly specified.









    share|improve this question
























      1












      1








      1








      In the package LTRCtrees one can fit a decision tree to the special format of the Surv function :Surv(time, time2, event) as per below example



      set.seed(0)
      library(survival)
      library(LTRCtrees)
      ## Create the start-stop-event triplet needed for coxph and LTRC trees
      first <- with(pbcseq, c(TRUE, diff(id) !=0)) #first id for each subject
      last <- c(first[-1], TRUE) #last id
      time1 <- with(pbcseq, ifelse(first, 0, day))
      time2 <- with(pbcseq, ifelse(last, futime, c(day[-1], 0)))
      event <- with(pbcseq, ifelse(last, status, 0))
      event <- 1*(event==2)

      pbcseq$time1 <- time1
      pbcseq$time2 <- time2
      pbcseq$event <- event

      ## Fit the Cox model and LTRC trees with time-varying covariates
      fit.cox <- coxph(Surv(time1, time2, event) ~ age + sex + log(bili), pbcseq)
      LTRCIT.fit <- LTRCIT(Surv(time1, time2, event) ~ age + sex + log(bili), pbcseq)
      LTRCART.fit <- LTRCART(Surv(time1, time2, event) ~ age + sex + log(bili), pbcseq)


      Is it possible to use the same function on a Random Forest using the RandomForestSRC library



      library(randomForestSRC)
      RF.fit <- rfsrc(Surv(time1, time2, event) ~ age + sex + log(bili), data=pbcseq, nsplit = 3, ntree = 100, importance = TRUE)


      Which generates an error:



      Error in parseFormula(formula, data, ytry) : 
      Survival formula incorrectly specified.









      share|improve this question














      In the package LTRCtrees one can fit a decision tree to the special format of the Surv function :Surv(time, time2, event) as per below example



      set.seed(0)
      library(survival)
      library(LTRCtrees)
      ## Create the start-stop-event triplet needed for coxph and LTRC trees
      first <- with(pbcseq, c(TRUE, diff(id) !=0)) #first id for each subject
      last <- c(first[-1], TRUE) #last id
      time1 <- with(pbcseq, ifelse(first, 0, day))
      time2 <- with(pbcseq, ifelse(last, futime, c(day[-1], 0)))
      event <- with(pbcseq, ifelse(last, status, 0))
      event <- 1*(event==2)

      pbcseq$time1 <- time1
      pbcseq$time2 <- time2
      pbcseq$event <- event

      ## Fit the Cox model and LTRC trees with time-varying covariates
      fit.cox <- coxph(Surv(time1, time2, event) ~ age + sex + log(bili), pbcseq)
      LTRCIT.fit <- LTRCIT(Surv(time1, time2, event) ~ age + sex + log(bili), pbcseq)
      LTRCART.fit <- LTRCART(Surv(time1, time2, event) ~ age + sex + log(bili), pbcseq)


      Is it possible to use the same function on a Random Forest using the RandomForestSRC library



      library(randomForestSRC)
      RF.fit <- rfsrc(Surv(time1, time2, event) ~ age + sex + log(bili), data=pbcseq, nsplit = 3, ntree = 100, importance = TRUE)


      Which generates an error:



      Error in parseFormula(formula, data, ytry) : 
      Survival formula incorrectly specified.






      r random-forest survival-analysis






      share|improve this question













      share|improve this question











      share|improve this question




      share|improve this question










      asked Mar 7 at 12:27









      aieduaiedu

      306




      306






















          0






          active

          oldest

          votes











          Your Answer






          StackExchange.ifUsing("editor", function ()
          StackExchange.using("externalEditor", function ()
          StackExchange.using("snippets", function ()
          StackExchange.snippets.init();
          );
          );
          , "code-snippets");

          StackExchange.ready(function()
          var channelOptions =
          tags: "".split(" "),
          id: "1"
          ;
          initTagRenderer("".split(" "), "".split(" "), channelOptions);

          StackExchange.using("externalEditor", function()
          // Have to fire editor after snippets, if snippets enabled
          if (StackExchange.settings.snippets.snippetsEnabled)
          StackExchange.using("snippets", function()
          createEditor();
          );

          else
          createEditor();

          );

          function createEditor()
          StackExchange.prepareEditor(
          heartbeatType: 'answer',
          autoActivateHeartbeat: false,
          convertImagesToLinks: true,
          noModals: true,
          showLowRepImageUploadWarning: true,
          reputationToPostImages: 10,
          bindNavPrevention: true,
          postfix: "",
          imageUploader:
          brandingHtml: "Powered by u003ca class="icon-imgur-white" href="https://imgur.com/"u003eu003c/au003e",
          contentPolicyHtml: "User contributions licensed under u003ca href="https://creativecommons.org/licenses/by-sa/3.0/"u003ecc by-sa 3.0 with attribution requiredu003c/au003e u003ca href="https://stackoverflow.com/legal/content-policy"u003e(content policy)u003c/au003e",
          allowUrls: true
          ,
          onDemand: true,
          discardSelector: ".discard-answer"
          ,immediatelyShowMarkdownHelp:true
          );



          );













          draft saved

          draft discarded


















          StackExchange.ready(
          function ()
          StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fstackoverflow.com%2fquestions%2f55043807%2fr-fitting-survival-trees-with-time-varying-covariates-in-randomforestsrc%23new-answer', 'question_page');

          );

          Post as a guest















          Required, but never shown

























          0






          active

          oldest

          votes








          0






          active

          oldest

          votes









          active

          oldest

          votes






          active

          oldest

          votes















          draft saved

          draft discarded
















































          Thanks for contributing an answer to Stack Overflow!


          • Please be sure to answer the question. Provide details and share your research!

          But avoid


          • Asking for help, clarification, or responding to other answers.

          • Making statements based on opinion; back them up with references or personal experience.

          To learn more, see our tips on writing great answers.




          draft saved


          draft discarded














          StackExchange.ready(
          function ()
          StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fstackoverflow.com%2fquestions%2f55043807%2fr-fitting-survival-trees-with-time-varying-covariates-in-randomforestsrc%23new-answer', 'question_page');

          );

          Post as a guest















          Required, but never shown





















































          Required, but never shown














          Required, but never shown












          Required, but never shown







          Required, but never shown

































          Required, but never shown














          Required, but never shown












          Required, but never shown







          Required, but never shown







          Popular posts from this blog

          Can't initialize raids on a new ASUS Prime B360M-A motherboard2019 Community Moderator ElectionSimilar to RAID config yet more like mirroring solution?Can't get motherboard serial numberWhy does the BIOS entry point start with a WBINVD instruction?UEFI performance Asus Maximus V Extreme

          Identity Server 4 is not redirecting to Angular app after login2019 Community Moderator ElectionIdentity Server 4 and dockerIdentityserver implicit flow unauthorized_clientIdentityServer Hybrid Flow - Access Token is null after user successful loginIdentity Server to MVC client : Page Redirect After loginLogin with Steam OpenId(oidc-client-js)Identity Server 4+.NET Core 2.0 + IdentityIdentityServer4 post-login redirect not working in Edge browserCall to IdentityServer4 generates System.NullReferenceException: Object reference not set to an instance of an objectIdentityServer4 without HTTPS not workingHow to get Authorization code from identity server without login form

          2005 Ahvaz unrest Contents Background Causes Casualties Aftermath See also References Navigation menue"At Least 10 Are Killed by Bombs in Iran""Iran"Archived"Arab-Iranians in Iran to make April 15 'Day of Fury'"State of Mind, State of Order: Reactions to Ethnic Unrest in the Islamic Republic of Iran.10.1111/j.1754-9469.2008.00028.x"Iran hangs Arab separatists"Iran Overview from ArchivedConstitution of the Islamic Republic of Iran"Tehran puzzled by forged 'riots' letter""Iran and its minorities: Down in the second class""Iran: Handling Of Ahvaz Unrest Could End With Televised Confessions""Bombings Rock Iran Ahead of Election""Five die in Iran ethnic clashes""Iran: Need for restraint as anniversary of unrest in Khuzestan approaches"Archived"Iranian Sunni protesters killed in clashes with security forces"Archived