限制性立方样条图：适用于线性回归、logistic回归、cox回归

xujunzju • 2024年6月24日 20:27 • 统计分析 • 阅读 227

官网教程：https://github.com/KunHuo/plotRCS，安装指令：

# 安装 devtools（如果您还没有安装）
install.packages("devtools")

devtools::install_github("kunhuo/plotRCS")

线性回归模型的 RCS 曲线

library(plotRCS)

# 查看数据
head(cancer)

# 线性回归模型的 RCS 曲线
rcsplot(data = cancer,
        outcome = "size",
        exposure = "age",
        covariates = c("sex", "race", "metastasis"))
        
res.line <- rcsplot(data = cancer,
                    outcome = "size",
                    exposure = "age",
                    covariates = c("sex", "race", "metastasis"))

res. line$data

逻辑回归模型的 RCS 曲线

# RCS curves for a logistic regression model
rcsplot(data = cancer,
        positive = "Dead",#可以在逻辑回归中人为指定事件
        outcome = "status",
        exposure = "age",
        covariates = c("sex", "race", "size", "metastasis"))

Cox 回归模型的 RCS 曲线

rcsplot(data = cancer,
        outcome = "status",
        time = "time",
        exposure = "age",
        covariates = c("sex", "race", "size", "metastasis"))