Thanks a lot for the reply! On Tue, Dec 31, 2013 at 1:59 AM, Göran Broström <goran.brost...@umu.se>wrote:
> On 12/30/2013 11:04 PM, Jieyue Li wrote: > >> Dear All, >> >> I want to have the cumulative incidence curves for 'mstate' data using >> Survival package in R. But I got some problems: >> I. Problem 1: >> 1. If I only use intercept without any covariates, I can have 'right' >> cumulative incidence curves (2 for 2 competing risks): >> library(Survival) >> > > That shouldn't work;) This is from an example from the Survival package... > > fitCI <- survfit(Surv(stop, status*as.numeric(event), type="mstate") ~ >> 1,data=mgus1, subset=(start==0)) >> plot(fitCI) >> 2. If I include one variate ('sex'), I get 4 curves (attached; I guess >> because there are two levels in 'sex' and 2 competing risks): >> fitCI <- survfit(Surv(stop, status*as.numeric(event), type="mstate") >> ~sex,data=mgus1, subset=(start==0)) >> plot(fitCI) >> However, I want to just have 2 cumulative incidence curves estimated from >> several covariates (such as 'sex', 'age', 'alb', etc. in mgus1). Could you >> please help me to do that? Thank you very much! >> > > I suggest that you check the Task Views, under 'Survival' and 'Multistate > Models', for instance the 'cmprsk' and 'timereg' packages. > > > II. Problem 2: >> I try using an example from sourcecode.pdf: >> fit <- survfit(Surv(time, status, type=mstate) ~ sex, data=mine) >> but where can I have the 'mine' data? Thank you! >> > > Where do you find 'sourcecode.pdf'? > It's from http://stat.ethz.ch/R-manual/R-patched/library/survival/doc/sourcecode.pdf > > Göran Broström > > >> Best, >> >> Jieyue >> >> > ______________________________________________ > R-help@r-project.org mailing list > 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]]
______________________________________________ R-help@r-project.org mailing list 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.