[R] Plotting adjusted KM curve

[John Fox]

On 2021-04-04 10:45 p.m., John Fox wrote:
> Dear John,
> 
> I think that what you're looking for is
> 
> plot(survfit(fit1Cox, newdata=data.frame(age=rep(65, 2), 
> sex=factor("female", "male"))))

Whoops, that should be

plot(survfit(fit1Cox, newdata=data.frame(age=rep(65, 2), 
sex=factor(c("female", "male")))))

John


> 
> assuming, of course, that sex is a factor with levels "female" and "male".
> 
> I hope this helps,
> John
> 
> John Fox, Professor Emeritus
> McMaster University
> Hamilton, Ontario, Canada
> web: https://socialsciences.mcmaster.ca/jfox/
> 
> On 2021-04-04 9:34 p.m., Sorkin, John wrote:
>> Colleagues,
>> I am using the coxph to model survival time. How do I plot an adjusted 
>> Kaplan Meir plot resulting from coxph? The code I would like to run 
>> would start with:
>>
>> # run cox model
>> fit1Cox <- coxph(surv_object ~age+sex,data=mydata)
>>
>> I have no idea what would follow.
>>
>> I would like to plot adjusted KM curves for men vs. women at age 65.
>>
>> Thank you,
>> John
>>
>>
>>     [[alternative HTML version deleted]]
>>
>> ______________________________________________
>> R-help using 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.
>>
> 
> ______________________________________________
> R-help using 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.