Hello, I am using the function survfit in the 'survival' package. Calling the function produces the median survival time automatically, as below. sleepfit <- survfit(Surv(timeb, death)~1)> sleepfitCall: survfit(formula = Surv(timeb, death) ~ 1) records n.max n.start events median 0.95LCL 0.95UCL 10761.0 10761.0 10761.0 5289.0 15.5 15.5 15.5 I was wondering, is it possible to also output the quartiles Q1 and Q3 survival times? Regards, Lucy [[alternative HTML version deleted]]
I don't have access to a computer that will run R at the moment but if I were tackling this, I would first look at the code for summary.survfit and adapt the code that extracts the median ( the 50th percentile ). -- David Sent from my iPhone> On Mar 15, 2014, at 1:43 PM, Lucy Leigh <lucy.leigh at newcastle.edu.au> wrote: > > Hello, > > I am using the function survfit in the 'survival' package. Calling the function produces the median survival > time automatically, as below. > sleepfit <- survfit(Surv(timeb, death)~1) >> sleepfit > Call: survfit(formula = Surv(timeb, death) ~ 1) > > records n.max n.start events median 0.95LCL 0.95UCL > 10761.0 10761.0 10761.0 5289.0 15.5 15.5 15.5 > > I was wondering, is it possible to also output the quartiles Q1 and Q3 survival times? > > Regards, > Lucy > > [[alternative HTML version deleted]] > > ______________________________________________ > R-help at r-project.org mailing list > 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.
This may be of interest: http://stackoverflow.com/questions/20724872/median-from-survfit-object-and-textconnection Sent from my iPhone> On Mar 15, 2014, at 1:43 PM, Lucy Leigh <lucy.leigh at newcastle.edu.au> wrote: > > Hello, > > I am using the function survfit in the 'survival' package. Calling the function produces the median survival > time automatically, as below. > sleepfit <- survfit(Surv(timeb, death)~1) >> sleepfit > Call: survfit(formula = Surv(timeb, death) ~ 1) > > records n.max n.start events median 0.95LCL 0.95UCL > 10761.0 10761.0 10761.0 5289.0 15.5 15.5 15.5 > > I was wondering, is it possible to also output the quartiles Q1 and Q3 survival times? > > Regards, > Lucy > > [[alternative HTML version deleted]] > > ______________________________________________ > R-help at r-project.org mailing list > 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.