search

andybega / spduration

Split-population duration model in R
https://andybeger.com/spduration/
4 stars 0 forks source link

predict weibull >1 #15

Closed andybega closed 6 years ago

andybega commented 9 years ago

Happens around lines 76-78:

if (stat == "conditional hazard") res <- atrisk.t * ft/pmax(1e-10, (cure.t + atrisk.t * st))

from spduration-paper example debugging:

Browse[2]> s0[2014] 2562 0.007962334 Browse[2]> st[2014] 2562 0.001401938 Browse[2]> ft[2014]/st[2014] 2562 1.752906