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

          Thal And Out Agency railway station See also References External links Navigation menuOfficial Web Site of Pakistan RailwaysArchivedOfficial Web Site of Pakistan Railwayseeexpanding ite

          Posting a File and Associated Data to a RESTful WebService preferably as JSON2019 Community Moderator ElectionREST URI convention - Singular or plural name of resource while creating itUploading both data and files in one form using Ajax?How do I upload a file with metadata using a REST web service?REST API - file (ie images) processing - best practicesSpring: JSON data and file in the same requestREST API Put large dataPosting a file and JSON data to Spring rest serviceRESTful Web Service Upload/Download Large Data With JSONREST API Design sending JSON data and a file to the api in same requestPost JSONArray to REST servicePUT vs. POST in RESTParsing values from a JSON file?JavaScript/jQuery to download file via POST with JSON dataHow do I upload a file with metadata using a REST web service?How to POST JSON data with Curl from Terminal/Commandline to Test Spring REST?Separate REST JSON API server and client?what's the correct way to send a file from REST web service to client?How do I write JSON data to a file?What is “406-Not Acceptable Response” in HTTP?REST API - file (ie images) processing - best practices

          How can I change the color of pagination dots of UIPageControl?How to change UIPageControl dotsIs there a way to change page indicator dots colorCustomize dot with image of UIPageControl at index 0 of UIPageControlNo visible @interface for 'NSObject<PageControlDelegate>' declares the selector 'pageControlPageDidChange:'How to change the color of pagination dots in UIPageControl with a different color per pagepagecontrol indicator custom image instead of DefaultChanging the colour of UIPageControl dots in MonoTouchpagecontrol selectable page visibility color?Alternative way to load ViewControllers on a UIPageControlHow to set only layer.border-color for UIpage control dots in swiftHow can I develop for iPhone using a Windows development machine?How to change the name of an iOS app?UITableView - change section header coloruipagecontrol indicator(dot)issueCustom UIPageControl dots color not changingchange the interspace between UIPageControl dotsHow to change Status Bar text color in iOSHow can I change image tintColor in iOS and WatchKitUIPageControl dots with larger space in between each dotsHow to change the color of pagination dots in UIPageControl with a different color per page