Site icon R-bloggers

What if there’s no censoring?…

[This article was first published on R | Gianluca Baio, and kindly contributed to R-bloggers]. (You can report issue about the content on this page here)
Want to share your content on R-bloggers? click here if you have a blog, or here if you don't.

Philip Cooney has alerted me to a slight inconsistency in survHE. Basically, he was trying to work in a setting where the data contain no censoring, using the HMC/rstan models generated by survHE. All works OK when using some of the models implemented, but failed when trying to fit the Gamma, Generalised Gamma and Generalised F.

The reason for this is that, for simplicity (and actually, our empirical tests show that computational run-time is much reduced) these models construct the log-likelihood by separating the set of individuals who are fully observed from the one of individuals who are censored. Basically, the model splits the data into $\boldsymbol{y}=(y^{\rm{obs}},y^{\rm{cens}})$, when $y^{\rm{cens}}$ is a latent parameter and then considers $$y^{\rm{obs}}_i \sim p(y^{\rm{obs}}\mid \mu_i,\alpha) \qquad {\sf and} \qquad g(\mu_i) = \boldsymbol{X}^{\rm{obs}}_i \boldsymbol\beta$$ and $$y^{\rm{cens}}_i \sim p(y^{\rm{cens}}\mid \mu_i,\alpha) \qquad {\sf and} \qquad g(\mu_i) = \boldsymbol{X}^{\rm{obs}}_i \boldsymbol\beta + \log(d_i),$$ where $d_i$ is the censored time (so the point at which the individuals have exited from the study without experiencing the event).

Because the two variables share the parameters, effectively the estimates account for censoring in this formulation. But, for these three models, survHE couldn’t handle in its default implementation, the case in which the number of censored data was 0.

Philip and I have fixed the code and I’ve now updated the GitHub devel version — I’ve upped the version number to 1.1.3 and included Philip as a contributor.

To leave a comment for the author, please follow the link and comment on their blog: R | Gianluca Baio.

R-bloggers.com offers daily e-mail updates about R news and tutorials about learning R and many other topics. Click here if you're looking to post or find an R/data-science job.
Want to share your content on R-bloggers? click here if you have a blog, or here if you don't.