在学习一篇达格列净文献的时候（下图）显示，无论患者BMI大小如何，与安慰剂相比，达格列净对于主要终点事件的获益一致。文章（doi:10.1002/ejhf.2308）方法学说是使用stata做的限制性立方样条图。

医咖会讲授的方法主要是 连续变量 之于终点事件的 RCS，而这个治疗相关 （treatment A vs. treatment B）的 RCS 都没找到正确的实现方法。想请教一下，这个图的实现方法 SAS、R、stata...都可以~ 非常感谢