With model matrix. Takes 20–60 mins to fit.