在网上找到了一套算法，但是只适用于二分类变量，没找到多分类变量的求解方法想请问下R语言大神，我的问题有没有算法可以解决，谢谢

以下是二分类变量可行的time-independent HR作图算法
mod_tvc <- rstpm2::stpm2(Surv(time, status==2) ~ sex, #需要比较的变量
data = lung, #数据集名称
df=3,
tvc = list(sex = 1))#变量平滑水平

mod_ph <- rstpm2::stpm2(Surv(time, status==2) ~ sex,
                        data = lung)#不考虑时间变量影响的HR

plot(mod_tvc, #刚刚构建的模型
newdata = data.frame(sex = 2), #比对数据
type = "hr",
var = "sex",#需要比较的变量，同模型
ci = TRUE,#置信区间
rug = FALSE,
main = "Time-dependent HR",
ylab = "Hazard ratio", xlab = "Time")

plot(mod_ph,  newdata = data.frame(sex = 2), type = "hr",
     var = "sex", add = TRUE, ci = FALSE, lty = 2)