According to the guidelines (if I'm not mistaken), the code below is sufficient (without any specification) to give Kaplan-Meier curves with censored data markings on Kaplan-Meier curves. But in my case censored data don't appears on the curves?! library(survival) mydata<-read.csv (file="C:/mydata/mydata.csv", header=TRUE, sep=";" ) # Sic! The separator in my csv file is ";" dput (mydata, "dputmydata.r") #attached Y <- Surv (mydata$time, mydata$status == 2) # 2 -- encodes event km <- survfit (Y~mydata$stage) plot (km)
Have you consulted ?plot.survfit ? There are examples for KM plots there. Also, obvious question: Have you specfied the censoring properly in your data and fit? 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 Thu, Jan 17, 2019 at 11:39 AM Medic <mailiPadpost at gmail.com> wrote:> According to the guidelines (if I'm not mistaken), the code below is > sufficient (without any specification) to give Kaplan-Meier curves with > censored data markings on Kaplan-Meier curves. But in my case censored data > don't appears on the curves?! > > library(survival) > mydata<-read.csv (file="C:/mydata/mydata.csv", header=TRUE, sep=";" ) > # Sic! The separator in my csv file is ";" > > dput (mydata, "dputmydata.r") > #attached > > Y <- Surv (mydata$time, mydata$status == 2) > # 2 -- encodes event > > km <- survfit (Y~mydata$stage) > > plot (km) > ______________________________________________ > 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. >[[alternative HTML version deleted]]
Hi, Just to emphasize Bert's e-mail with a hint, the 'mark.time' argument for plot.survfit() is FALSE by default. Regards, Marc Schwartz> On Jan 17, 2019, at 2:45 PM, Bert Gunter <bgunter.4567 at gmail.com> wrote: > > Have you consulted ?plot.survfit ? There are examples for KM plots there. > > Also, obvious question: Have you specfied the censoring properly in your > data and fit? > > > 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 Thu, Jan 17, 2019 at 11:39 AM Medic <mailiPadpost at gmail.com> wrote: > >> According to the guidelines (if I'm not mistaken), the code below is >> sufficient (without any specification) to give Kaplan-Meier curves with >> censored data markings on Kaplan-Meier curves. But in my case censored data >> don't appears on the curves?! >> >> library(survival) >> mydata<-read.csv (file="C:/mydata/mydata.csv", header=TRUE, sep=";" ) >> # Sic! The separator in my csv file is ";" >> >> dput (mydata, "dputmydata.r") >> #attached >> >> Y <- Surv (mydata$time, mydata$status == 2) >> # 2 -- encodes event >> >> km <- survfit (Y~mydata$stage) >> >> plot (km)