कोई भी जानता है कि जीवित रहने के विश्लेषण में ggplot या जाली का लाभ कैसे लें? ट्रेली या पहलू-जैसे अस्तित्व ग्राफों को करना अच्छा लगेगा।ggplot या जाली में सर्व ऑब्जेक्ट का उपयोग
तो अंत में मैंने चारों ओर खेला और एक कपलान-मेयर प्लॉट के लिए एक समाधान मिला। सूची तत्वों को डेटाफ्रेम में लेने में गन्दा कोड के लिए मैं क्षमा चाहता हूं, लेकिन मैं एक और तरीका नहीं समझ सका।
नोट: यह केवल स्तर के दो स्तरों के साथ काम करता है। यदि कोई जानता है कि मैं ऐसा करने के लिए x<-length(stratum)
का उपयोग कैसे कर सकता हूं तो कृपया मुझे बताएं (स्टेटा में मैं एक मैक्रो-अनिश्चितता में शामिल हो सकता हूं कि यह आर में कैसे काम करता है)।
ggkm<-function(time,event,stratum) {
m2s<-Surv(time,as.numeric(event))
fit <- survfit(m2s ~ stratum)
f$time <- fit$time
f$surv <- fit$surv
f$strata <- c(rep(names(fit$strata[1]),fit$strata[1]),
rep(names(fit$strata[2]),fit$strata[2]))
f$upper <- fit$upper
f$lower <- fit$lower
r <- ggplot (f, aes(x=time, y=surv, fill=strata, group=strata))
+geom_line()+geom_ribbon(aes(ymin=lower,ymax=upper),alpha=0.3)
return(r)
}
रेमन Saccilotto लिखा एक ggplot2 ट्यूटोरियल कि ggplot2 में के.एम. भूखंडों के लिए कार्य शामिल हैं: http://www.ceb-institute.org/bbs/wp-content/uploads/2011/09/handout_ggplot2.pdf – MattBagg