R语言2种方法绘制分类的RCS曲线

本期文章应该算是一个答疑。
限制性立方样条函数（RCS）在比较非线性关系中很常用。我们既往已经《实例演示R语言制作限制性立方条图》，后台有粉丝问如何绘制分类的RCS曲线，今天我们来演示一下，R语言绘制多分类的RCS曲线，并回答一下之前粉丝的部分问题。

继续使用我们的乳腺癌数据（公众号回复：乳腺癌，可以获得数据）。我们先导入数据和R包。

library(foreign)
library(rms)
bc <- read.spss("E:/r/test/Breast cancer survival agec.sav",
                use.value.labels=F, to.data.frame=T)
bc <- na.omit(bc)
names(bc)
attach(bc)

这是个乳腺癌患者的生存数据，age表示年龄，pathsize表示病理肿瘤大小（厘米），lnpos表示腋窝淋巴结阳性，histgrad表示病理组织学等级，er表示雌激素受体状态，pr表示孕激素受体状态，status结局事件是否死亡，pathscat表示病理肿瘤大小类别（分组变量），ln_yesno表示是否有淋巴结肿大，time是生存时间，后面的agec是我们自己设定的，不用管它。

本文为转载文章，原文章地址如下：

https://mp.weixin.qq.com/s?__biz=MzI1NjM3NTE1NQ==&mid=2247485765&idx=1&sn=c17ca7f01e0d84535834462e30681f73&chksm=ea26e959dd51604fb821ad19f91b074e9103fba8e20d4a2e806f6bbf319bee604a4f886dc67c&payreadticket=HLy4Nnf3c52CoCwY92i37tzhAeP4tkdhJU6zgZA-OzZFyjfqJPejSeJfZPlaMeAdOM1MHLk#rd