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]]
Hi Google answered https://rdrr.io/bioc/survcomp/man/km.coxph.plot.html Is it what do you want? Cheers Petr> -----Original Message----- > From: R-help <r-help-bounces at r-project.org> On Behalf Of Sorkin, John > Sent: Monday, April 5, 2021 3:35 AM > To: r-help at r-project.org (r-help at r-project.org) <r-help at r-project.org> > Subject: [R] Plotting adjusted KM curve > > Colleagues, > I am using the coxph to model survival time. How do I plot an adjustedKaplan> Meir plot resulting from coxph? The code I would like to run would startwith:> > # 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 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.
On 2021-04-05 03:34, 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.You should look at ?survfit.coxph in the survival package, especially the 'newdata' argument.> > I would like to plot adjusted KM curves for men vs. women at age 65.Then I guess that you should stratify on sex: fit <- coxph(surv_object ~ age + strata(sex), data = mydata) sfit <- survfit(fit, newdata = data.frame(age = 65)) plot(sfit) HTW, G?ran> > Thank you, > John > > > [[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. >