如何用stata做限制性立方样条图变量，并建模？