[R] How to get CI from surfit object in survival
Adrian Johnson
oriolebaltimore at gmail.com
Sun Aug 27 05:30:21 CEST 2017
Thank you. I checked everything else but my own argument. :-) that
is silly. Thanks
On Sat, Aug 26, 2017 at 8:38 PM, Bert Gunter <bgunter.4567 at gmail.com> wrote:
> Did you not notice the conf.type = "none" argument to your survfit
> call and the associated documentation in the survfit help?
>
> -- Bert
>
>
> Bert Gunter
>
> "The trouble with having an open mind is that people keep coming along
> and sticking things into it."
> -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip )
>
>
> On Sat, Aug 26, 2017 at 5:18 PM, Adrian Johnson
> <oriolebaltimore at gmail.com> wrote:
>> Dear Bert,
>> thank you for suggestion. I am aware of R-help function. I must
>> apologize, my earlier question could lead to assumptions otherwise.
>>
>> As you can see below, I only get Std. error but not lower and Upper CIs.
>>
>> I was wondering if there is another argument or method, could give CIs
>> which I cannot find anywhere.
>>
>> Thanks
>> Adrian
>>
>>
>>> Nsurv <- survfit(Surv(NOSM,NStat)~as.factor(expcat1),conf.type='none',type='kaplan-meier')
>>
>>
>>> summary(Nsurv,time=c(12,24))
>> Call: survfit(formula = Surv(NOSM, NStat) ~ as.factor(expcat1),
>> conf.type = "none", type = "kaplan-meier")
>>
>> 1 observation deleted due to missingness
>> as.factor(expcat1)=DN
>> time n.risk n.event survival std.err
>> 12 36 15 0.711 0.0629
>> 24 21 8 0.535 0.0722
>>
>> as.factor(expcat1)=UP
>> time n.risk n.event survival std.err
>> 12 20 9 0.704 0.0833
>> 24 9 9 0.359 0.0931
>>
>>> sessionInfo()
>> R version 3.2.3 (2015-12-10)
>> Platform: i386-w64-mingw32/i386 (32-bit)
>> Running under: Windows 8.1 x64 (build 9600)
>>
>> locale:
>> [1] LC_COLLATE=English_United States.1252 LC_CTYPE=English_United
>> States.1252 LC_MONETARY=English_United States.1252 LC_NUMERIC=C
>> [5] LC_TIME=English_United States.1252
>>
>> attached base packages:
>> [1] stats graphics grDevices utils datasets methods base
>>
>> other attached packages:
>> [1] survival_2.39-5
>>
>> loaded via a namespace (and not attached):
>> [1] Matrix_1.2-6 splines_3.2.3 grid_3.2.3 lattice_0.20-33
>>
>> On Sat, Aug 26, 2017 at 7:24 PM, Bert Gunter <bgunter.4567 at gmail.com> wrote:
>>> ???
>>>
>>> Both ?survit.object (linked in the see also section of ?survfit) and
>>> ?summary.survfit give you this information. Do you not know how to use
>>> R's help faciities -- in which case you should learn them now; see
>>> ?help -- or have I misunderstood your query?
>>>
>>> Cheers,
>>> Bert
>>>
>>>
>>> Bert Gunter
>>>
>>> "The trouble with having an open mind is that people keep coming along
>>> and sticking things into it."
>>> -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip )
>>>
>>>
>>> On Sat, Aug 26, 2017 at 4:03 PM, Adrian Johnson
>>> <oriolebaltimore at gmail.com> wrote:
>>>> Hi
>>>> I am not sure hiw to get thr CI from summary function on the survfit object.
>>>> I can get the percent survival for 2 years from survfit object but I dont
>>>> get Confidence intervals
>>>> Could anyone suggest a hint
>>>> Thanks
>>>> Adrian
>>>>
>>>> [[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.
More information about the R-help
mailing list