Liangint function

Fit a semi-parametric joint model, incorporating intercept estimation