Model segmental linear regression

Thanks! There is a typo in “mu=pm.math.switch(y<bp1,mu1,pm.math.switch(y<bp2,mu2,mu3))” . “y” should be “t”. I didn’t notice that. After changing it, it works perfectly! I am so confident about pymc