Dear Enzo and Phil,
Many thanks for your quick replies, most helpful.
I have now downloaded the paper in Stata journal and will read it carefully.
Would it be possible that Enzo sends me -stcompadj- to try?
Many thanks again.
With best wishes,
Natalie
On Sun, Apr 5, 2009 at 7:29 AM, Enzo Coviello <[email protected]> wrote:
> Hi Natalie,
>
> the Phil's answer is complete about your problem.
>
> I would add that you can download the paper where the use of -stcompet- is
> illustrated from
> http://www.statajournal.com/sjpdf.html?articlenum=st0059
>
> Here some practical example can help you in using the output of -stcompet-
> for making graphs.
>
> Note that -stcompet- does not estimate the cause specific mortality rates
> but the cumulative incidence function in the presence of the competing
> risks, a quite different function.
>
> I wrote on this topic a further macro -stcompadj- that computes the
> covariate-adjusted cumulative incidence according to a method proposed by
> Rosthoj and Andersen. The macro is not yet available on the SSC Archive
> because I would add a delta-method estimate of the variance of this
> function. At present -stcompadj- computes the bootstrap confidence intervals
> of the adjusted cumulative incidence that seem quite good but I am not
> completely convinced about their validity.
> I would pleased to send privately to other Stata users -stcompadj- for
> further checks and, eventually, for adding the delta-method variance
> estimate of the covariate-adjusted cumulative incidence.
>
> Best wishes.
>
> Enzo
>
>
>
> natalie chan ha scritto:
>>
>> Dear Statalister,
>>
>> I am trying to use Stcompet to determine the cause specific mortality
>> rate on a dataset with about 3500 subjects, 30000 observations and 3
>> possible causes.
>>
>> I use stset first:
>>
>> stset time1, id(id) time0(time0) failure(path==1)
>>
>> the results are as follows:
>>
>> id: id
>> failure event: path == 1
>> obs. time interval: (time0, time1]
>> exit on or before: failure
>>
>>
>> ------------------------------------------------------------------------------
>> 30367 total obs.
>> 2671 obs. begin on or after (first) failure
>>
>> ------------------------------------------------------------------------------
>> 27696 obs. remaining, representing
>> 3521 subjects
>> 563 failures in single failure-per-subject data
>> 27696 total analysis time at risk, at risk from t = 0
>> earliest observed entry t =1
>> last observed exit t =15
>>
>>
>> then I use stcompet:
>>
>> stcompet test1=ci test2=se, compet1(2) compet2(3)
>>
>> In the new variables generated, there are numbers when the failure
>> variable (=1, 2 or 3), but in other rows, the value is missing.
>>
>> I have searched Statalist Archive and found that someone had the same
>> problem but I did not find the answer. I would be grateful if anyone
>> could help me on this.
>>
>> Many thanks.
>>
>> Natalie
>> *
>> * For searches and help try:
>> * http://www.stata.com/help.cgi?search
>> * http://www.stata.com/support/statalist/faq
>> * http://www.ats.ucla.edu/stat/stata/
>>
>>
>
>
> --
> Enzo Coviello
> Unità di Epidemiologia e Statistica ASL Bari
> Piazza V. Emanuele, 14 - 70054 GIOVINAZZO (BA)
> Italy
> tel +39 080 3357883 - fax +39 080 3357867
> mobile 347 5016016 ^^^^
> Home
> via Imbriani, 186 - 70055 MINERVINO MURGE (BA)
> tel +39 0883 695055
>
> *
> * For searches and help try:
> * http://www.stata.com/help.cgi?search
> * http://www.stata.com/support/statalist/faq
> * http://www.ats.ucla.edu/stat/stata/
>
*
* For searches and help try:
* http://www.stata.com/help.cgi?search
* http://www.stata.com/support/statalist/faq
* http://www.ats.ucla.edu/stat/stata/
