跪求COX模型限制性立方样条实现代码（http://www.360doc.com/content/19/0421/14/47874886_830331378.shtml这个网址中的怎么实现不了）