-----Original Message-----
From: r-help-bounces at r-project.org [mailto:r-help-bounces at r-project.org] On Behalf Of R. Michael
Weylandt
Sent: Monday, October 17, 2011 8:51 AM
To: Cem Girit
Cc: r-help at r-project.org
Subject: Re: [R] Extracting results from a function output
Assuming I'm reading your post correctly -- fit is the output object
and the lines following >fit are the console output, which implicitly
calls print -- it looks like the median isn't stored in fit but is
calculated on the fly in the print call. But this doesn't seem
consistent with your claim that:
"Print" function does not reveal them either.
Could you perhaps clarify your statement about print(fit)?
Cheers from Princeton!
Michael Weylandt
On Mon, Oct 17, 2011 at 11:13 AM, Cem Girit <girit at biopticon.com> wrote:
Hello,
I am having hard time obtaining a value from a function. "fit" is a survival
function that produces some results, such as "median", "confidence
intervals" etc. But str() function does not list these values. How can I
extract these to be able use them? For example, I need "median" value for
the group DrugA which is 48. "Print" function does not reveal them either.
Thank you.
Call: survfit(formula = Surv(tT, dT) ~ gT, conf.type = "log-log")
records n.max n.start events median 0.95LCL 0.95UCL
gT=DrugA 9 9 9 6 48 32 NA
gT=DrugB 9 9 9 3 NA 42 NA
gT=DrugC 9 9 9 4 NA 42 NA
gT=Vehicle 9 9 9 8 40 37 45
List of 14
$ n : int [1:4] 9 9 9 9
$ time : num [1:28] 32 43 45 46 48 54 55 60 62 42 ...
$ n.risk : num [1:28] 9 8 7 6 5 4 3 2 1 9 ...
$ n.event : num [1:28] 1 1 1 1 1 0 1 0 0 2 ...
$ n.censor : num [1:28] 0 0 0 0 0 1 0 1 1 0 ...
$ surv : num [1:28] 0.889 0.778 0.667 0.556 0.444 ...
$ type : chr "right"
$ strata : Named int [1:4] 9 6 6 7
..- attr(*, "names")= chr [1:4] "gT=DrugA" "gT=DrugB" "gT=DrugC"
"gT=Vehicle"
$ std.err : num [1:28] 0.118 0.178 0.236 0.298 0.373 ...
$ upper : num [1:28] 0.984 0.939 0.878 0.805 0.719 ...
$ lower : num [1:28] 0.433 0.365 0.282 0.204 0.136 ...
$ conf.type: chr "log-log"
$ conf.int : num 0.95
$ call : language survfit(formula = Surv(tT, dT) ~ gT, conf.type =
"log-log")
- attr(*, "class")= chr "survfit"
Cem
Cem Girit, PhD
Biopticon Corporation
182 Nassau Street, Suite 204
Princeton, NJ 08542
Tel: (609)-853-0231
Email: ?<mailto:cgirit at biopticon.com> girit at biopticon.com
This email message and any attachments are confidential ...{{dropped:14}}