Statalist


[Date Prev][Date Next][Thread Prev][Thread Next][Date index][Thread index]

st: case-cohort analysis: stcascoh updated


From   Enzo Coviello <enzo.coviello@tin.it>
To   statalist@hsphsun2.harvard.edu
Subject   st: case-cohort analysis: stcascoh updated
Date   Wed, 27 Aug 2008 19:41:31 +0200

Dear Stata Users,

thanks to Matthias Mohner a bug has been found and fixed in -stcascoh-,
a command used for case-cohort analysis.

The new version is available for download on the SSC Archive.



Here is some more detail for people interested in this topic
The problem occurs when the data need not to be sampled. i.e when they
already consist of subcohort sample and of non-subcohort cases.

In this cases -stcascoh- can still be applied to appropriately generate
the observations for the next Cox model eventually with Self-Prentice
variance estimator (-stselpre- command).

The steps are as follows:

1) Divide the data set in two files: the first one for subcohort
observations, the second for additional cases.

2) Appropriately -stset- the new files.

3) Process the subcohort file with -stcascoh- setting alpha just below 1
(i.e. .999) and save this file.
Do some assert to verify that _subco _wSelpre _wBarlow are coded as you
expect.

4) Process the additional cases file with -stcascoh- setting alpha = 0
and save this second file.
Also for this file some -assert- is appropriate to verify that _subco
_wSelpre _wBarlow are coded as you expect.

5) Append the two files you generated.

6) Adapt the Barlow weights.
Suppose you sampled 1000 out of 16000 subjects.
You should type
replace _wBarlow=ln(16000/1000) if _wBarlow>0


Well, after using these codes, the old version of -stcascoh- does not
allow to use -stselpre-. This has been fixed.



Best wishes

Enzo

--
Enzo Coviello
Unità di Epidemiologia e Statistica ASL BARI
piazza V. Emanuele 14
70054 Giovinazzo (Bari)
Italy
mobile +39 345 5016016
tel +39 080 3357883 - fax +39 080 3357867
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
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/



© Copyright 1996–2014 StataCorp LP   |   Terms of use   |   Privacy   |   Contact us   |   What's new   |   Site index