Hello everyone.
I am working on cure models using stpm2.
I want to estimate the cure point and wrote the program as below, but the figure is something odd...
Does my program have some mistake? Does someone have any good idea?
or Is this model inappropriate for this annalysis?

-------------------------------------------------------------------------------------------------------
use http://pauldickman.com/data/colon.dta, clear

stset surv_mm, failure(status=1 2) scale(12) exit(time 10*12)
gen _age = min(int(age + _t),99)
gen _year = int(yydx + mmdx/12 + _t)
merge m:1 _year sex _age using http://pauldickman.com/data/popmort
keep if _merge==3

stpm2 if stage==1, scale(hazard) df(5) bhazard(rate) cure

predict rs_all, survival
predict cure1, cure
predict uncure, survival uncured
twoway line rs_all cure1 uncure surv_mm, sort


Array