Thanks a lot for the reply!

On Tue, Dec 31, 2013 at 1:59 AM, Göran Broström <goran.brostrom@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]]


