Fitting illness-death model for interval-censored data