x1(二分类变量)针对y(二分类变量)得or/hr随着x2(连续变量)得变量和切点

老师,我想看x1(二分类变量)针对y(二分类变量)得or/hr随着x2(连续变量)得变量和切点,这样的图可以用限制性立方样条回归做出来嘛?还是说其它方法,或者只能对x2分层后求每一层得参数。用什么r包可以完成呢?