Notice: On April 23, 2014, Statalist moved from an email list to a forum, based at statalist.org.
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
st: SV: stcompet and cumulative incidence at a specific point in time
From
Katja Maretty Nielsen <[email protected]>
To
"[email protected]" <[email protected]>
Subject
st: SV: stcompet and cumulative incidence at a specific point in time
Date
Sun, 16 Jun 2013 06:42:42 +0000
Dear Enzo.
I've encountered the same problem, but was advised to use the psudovalue approach with the -stpci-, were you generate pseudo values at the time you are interested in:
gen cevent=(dssstatus==2)
label variable cevent "competing risk event"
stpci cevent, at(1 2 3 4 5)
glm pseudo1, fam(gauss) link(id) vce(robust)
glm pseudo2, fam(gauss) link(id) vce(robust)
glm pseudo3, fam(gauss) link(id) vce(robust)
glm pseudo4, fam(gauss) link(id) vce(robust)
glm pseudo5, fam(gauss) link(id) vce(robust)
Hope this helps.
Best regards
Katja Maretty Nielsen
MD, Phd-student
Dept. of Experimental Clinical Oncology
Aarhus University Hospital
Noerrebrogade 44, building 5
DK- 8000 Aarhus C
Tel.: +45 78462621 / 50904640
Mail: [email protected]
________________________________________
Fra: [email protected] [[email protected]] på vegne af Enzo Coviello [[email protected]]
Sendt: 16. juni 2013 08:14
Til: [email protected]
Emne: st: stcompet and cumulative incidence at a specific point in time
Hi All,
Is there a Stata command that can be used to list (at specific time
points e.g. -sts list, at (0 (1) 5)-)
the cumulative incidence function with confidence intervals from a
competing risk analysis (after the ?stcompet- command)?
As far as I know there is not a command.
The following code should approximately obtain what you need supposing
that the -ci- is the name of the variable containing the cumulative
incidence:
forval i = 1(1)5 {
qui su _t if _t<=`i'
local tobs = `r(max)'
qui su ci if _t<=`i'
local ciat = `r(max)'
if `i'==10 di _col(5) "at" _col(18) "Obs Time" _col(33) "Cum Inc"
di _col(5) `i' _col(20) `tobs' _col(35) `ciat'
}
Best
Enzo
--
Enzo Coviello
Epidemiology Unit - Cancer Registry ASL BT
Piazza Umberto 1
76121 BARLETTA (BT)
Italy
mobile +39 347 5016016
tel +39 0883 577329
fax +39 0883 577288
Home +39 0883 695055
--
Enzo Coviello
Epidemiology Unit - Cancer Registry ASL BT
Piazza Umberto 1
76121 BARLETTA (BT)
Italy
mobile +39 347 5016016
tel +39 0883 577329
fax +39 0883 577288
Home +39 0883 695055
*
* For searches and help try:
* http://www.stata.com/help.cgi?search
* http://www.stata.com/support/faqs/resources/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/faqs/resources/statalist-faq/
* http://www.ats.ucla.edu/stat/stata/