Fit joint model for event and censoring with cluster-level - View it on GitHub