怎么用stata做限制性立方样条图？有没有代码可以参考？谢谢