Hi everyone,
I am using -stpm2cr- to estimate the cumulative incidence function of cardiovascular outcomes in the presence of death as a competing risk. On the whole, I've been able to work through most of the syntax challenges and have most calculations working. However, I'm still running into a problem with the Cumulative Incidence Function Difference (-cifdiff- in the code). First, here's my model code that runs without issue:
However, when I try to predict the CIF difference at 30 days with the following code, I get an error I don't understand: -Maximum number of iterations exceeded.-
My time variable is set to 1 month and I want to predict the differences for CV outcomes and mortality. What is curious is that when I plot this for one year or more, I can get estimates and plotted results. I'd welcome any ideas as to what this error means and how I can work around it!
Thanks!
I am using -stpm2cr- to estimate the cumulative incidence function of cardiovascular outcomes in the presence of death as a competing risk. On the whole, I've been able to work through most of the syntax challenges and have most calculations working. However, I'm still running into a problem with the Cumulative Incidence Function Difference (-cifdiff- in the code). First, here's my model code that runs without issue:
Code:
stpm2cr ///
[CVOutcome: study_grp, scale(hazard) df(5) tvc(revasc_grp3) dftvc(5)] ///
[Death: study_grp, scale(hazard) df(5) tvc(revasc_grp3) dftvc(5)], ///
events(statusCVODTH3yrCR) cause(1 2) censvalue(0) eform level(95)
Code:
predict cifdiff_CVO_1m_un in 1, cifdiff1(study_grp 0) cifdiff2(study_grp 1) ci timevar(t01cr)
Thanks!
