[Date Prev][Date Next][Thread Prev][Thread Next][Date index][Thread index]
Re: st: stcascoh.ado - CASE-cohort
At 15.35 17/04/2003 +0200, you wrote:
I am using the stcascoh.ado to model a case-cohort dataset. To figure out
how things work and if things go along properly, I tried to reproduce the
example of BArlow (Barlow WE, Ichicawa L, Rosner D, and Izumi S: Analysis
of Case-Cohort Designs; Journal Clinical Epidemiology 1999; 52: 1165-1172)
from the Breslow and Day nickel dataset.
My problem is that I am not getting there and cannot reproduce the results
(I tried it in SAS and it worked).
In this statement you are using as time axis time since first employment.
At examining the dataset after - stcascoh - I had the impression that
weights were wrongly assigned in the programm in lines 179
and 184. But I am not very familar with the programming language.
stset age_stop,failure(icd==160) id(id) enter( age_strt) exit( age_stop)
Are you sure that in SAS codes time axis is the same?
If SAS' time axis is age you should omit "origin( age_emp)"
If this check confirms that Stata's results do not match SAS' results,
please send me privately the output you obtained using SAS and I'll do some
further inspection in stcascoh.ado.
* For searches and help try: