Wollschlaeger, Daniel
2016-Aug-22 10:50 UTC
[R] problem concernig Survsplit, package survival
Thanks for bringing this issue in the book's description of survSplit() to my attention. It seems the change to the behavior of survSplit() was introduced in survival version 2.39-2. Up to (including) version 2.38-3, no formula was required if arguments 'end' and 'event' were specified. A corrected pdf of the chapter is available at http://dwoll.de/r/gddmr/09_survival.pdf Daniel Wollschlaeger On August 19, 2016 4:21:31 AM PDT, "Vinzenz V?lkel via R-help" <r-help at r-project.org> wrote:>Dear R-help-community, > > > >I hope, that?s the appropiate channel to post a quastion? > > > >For some days I have been struggling with a problem concerning the >?survSplit?-function of the package ?survival?. Searching the internet >I >have found a pretty good -German- description of Daniel Wollschl?ger >describing how to use survSplit: > > > >(https://books.google.de/books?id=9G4oBgAAQBAJ ><https://books.google.de/books?id=9G4oBgAAQBAJ&pg=PA345&lpg=PA345&dq=daniel+ >wollschl%C3%A4ger+survsplit&source=bl&ots=ajPgInB83d&sig=noQIAXMWe6oQJ-LwsTC >kroWdwHQ&hl=de&sa=X&ved=0ahUKEwiqwoX8q83OAhVlM5oKHSMGDrEQ6AEIHjAA#v=onepage& >q=daniel%20wollschl%C3%A4ger%20survsplit&f=false> >&pg=PA345&lpg=PA345&dq=daniel+wollschl%C3%A4ger+survsplit&source=bl&ots=ajPg >InB83d&sig=noQIAXMWe6oQJ-LwsTCkroWdwHQ&hl=de&sa=X&ved=0ahUKEwiqwoX8q83OAhVlM >5oKHSMGDrEQ6AEIHjAA#v=onepage&q=daniel%20wollschl%C3%A4ger%20survsplit&f=fal >se > > > >Mr. Wollschl?ger also provides the dataset, all of the code he used and >an >excerpt of his output: > > > > > >> library(survival) # f?r survSplit() > >> dfSurvCP <- survSplit(dfSurv, cut=seq(30, 90, by=30), end="obsT", > >+ event="status", start="start", id="ID", zero=0) > ># sortiere nach Beobachtungsobjekt und linken Intervallgrenzen > >> idxOrd <- order(dfSurvCP$ID, dfSurvCP$start) > >> head(dfSurvCP[idxOrd, ], n=7) > >obsT status sex X IV start ID > >1 30 0 f -1.3130607 A 0 >1 > >181 60 0 f -1.3130607 A 30 >1 > >361 63 1 f -1.3130607 A 60 >1 > >10 11 1 f -1.2282824 A 0 >10 > >100 27 1 m -0.1018403 B 0 >100 > >101 30 0 m -0.4079027 B 0 >101 > >281 42 1 m -0.4079027 B 30 >101 > > > > > >Unfortunaltely when I tried to copy Mr. Wollschl?gers operations using >his >very dataset, RStudio just returned the following error-message: > > > >> dfSurvCP2 <- >survSplit(dfSurv,cut=seq(30,90,by=30),end="obsT",event="status",start="start >",id="ID",zero=0) > >Error in formula.default(eval(parse(text = x, keep.source >FALSE)[[1L]])) : > > > invalid formula > > > >I tried to figure out what?s wrong, but so far didn?t find any >solutions? I >would be very glad if anyone could help me out here. > > > >vhfv > > > > > [[alternative HTML version deleted]] > > > >------------------------------------------------------------------------ > >______________________________________________ >R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see >https://stat.ethz.ch/mailman/listinfo/r-help >PLEASE do read the posting guide >http://www.R-project.org/posting-guide.html >and provide commented, minimal, self-contained, reproducible code.