实例演示R语言制作限制性立方条图

限制性立方样条函数（RCS）在比较非线性关系中很常用。既往我们已经讲过R语言制作限制性立方条图，但是讲得比较简单，中间有些环节没写出来，我也不是很满意，今天重新来说一下。主要是要用到rms包的rcs函数来绘制，继续用我们的乳腺癌数据来绘制，因为目前也没找到什么好的数据。
我们先导入包，导入数据

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)

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

转换分类变量并且抽取一部分变量等下来拟合模型

attach(bc)
be<-data.frame(age,status,time,ln_yesno)
be$ln_yesno<-as.factor(be$ln_yesno)

对数据进行打包，整理

dd <- datadist(be) #为后续程序设定数据环境
options(datadist='dd') #为后续程序设定数据环境

拟合模型

fit<- cph(Surv(time,status) ~ rcs(age,4) + ln_yesno,data=be)
an<-anova(fit)

生成预测值并且做图，我们要注意一下，这里的exp等于把它返回了HR值

plot(Predict(fit, age,fun=exp), anova=an, pval=T)

Nonlinear值为0.33，提示不是非线性，
论文中我们常可以看见有一条等于1的对比线，如下图这样，说明线上的值大于1，线下面小于1。我们可以通过HR等于1的age对它进行设定。

我们进行查看预测表，HR等于1，age为50

dd$limits$age[2] <-50
fit1=update(fit)

最后画出图形

HR<-Predict(fit, age,fun=exp,ref.zero = TRUE)
ggplot(HR)
ggplot()+geom_line(data=HR, aes(age,yhat),linetype=1,size=1,alpha = 0.9,colour="red")+geom_ribbon(data=HR, aes(age,ymin = lower, ymax = upper),alpha = 0.3,fill="red")+geom_hline(yintercept=1, linetype=2,size=1)+theme_classic()+ labs(title = "RCS", x="age", y="HR (95%CI)")

其他模型做法也是一样的，只是在构建模型这一块不同。
原创不易，求关注，求分享，求点击广告，各种求

实例演示R语言制作限制性立方条图相关推荐

R语言绘制限制立方条图2（基于logistic回归和线性回归）
临床上,因变量和临床的结局有时候不是线性关系,而回归模型有一个重要的假设就是自变量和因变量呈线性关联,因此非线性关系模型用回归分析来拟合受到限制.因此,一个更好的解决方法是拟合自变量与因变量之间的非线 ...
stata两种方法制作限制立方条图
我们前面已经讲了怎么使用R语言制作限制立方条图,今天来讲讲怎么使用stata制作限制立方条图,参考了网上常见的两种做法,顺便比较一下做图出来的差异,然后在和R语言比较一下. 首先要导入stata的xb ...
R语言----制作数据分布图（直方图+概率密度曲线）
R语言--制作group1和group2的分布图并输出tiff格式300dpi保存制作分布图我选择的是制作直方图+概率密度曲线~ library(ggplot2) data<-read.csv ...
手把手教你使用R语言制作临床决策曲线
DCA(Decision Curve Analysis)临床决策曲线是一种用于评价诊断模型诊断准确性的方法,在2006年由AndrewVickers博士创建,我们通常判断一个疾病喜欢使用ROC曲线的A ...
r语言 html gif,推荐：6步教你用R语言制作动图
原标题:推荐 :6步教你用R语言制作动图引言数据可视化可能是数据科学领域最重要却通常最少被提及的部分. 我这么说是因为创造数据故事和可视化对你的顾客最终怎么看待你的工作有巨大的影响.数据科学不仅仅 ...
R语言绘制环形树状图
R语言绘制环形树状图 1.主要用到dendextend和circlize包绘图: library(dendextend) library(circlize)# 距离矩阵 d <- dist(US ...
R语言绘制二维密度图
R语言绘制二维密度图二维密度图显示了两个数值变量之间的关系,一个在x轴上表示,另一个在Y轴上表示,与散点图类似,然后计算二维空间中特定区域内的观测数,并用颜色梯度表示.二维密度图有几种类型,以下主要 ...
R语言ggplot2可视化分面图使用facet_wrap函数和facet_grid函数实战
R语言ggplot2可视化分面图使用facet_wrap函数和facet_grid函数实战目录 R语言ggplot2可视化分面图使用facet_wrap函数和facet_grid函数实战
R语言ggplot2可视化树状图、层次聚类系统树图、树状图根据给定的距离度量将相似点分组在一起、并根据点的相似性将它们组织成树状图链接起来（Hierarchical Dendrogram）
R语言ggplot2可视化树状图.层次聚类系统树图.树状图根据给定的距离度量将相似点分组在一起.并根据点的相似性将它们组织成树状图链接起来(Hierarchical Dendrogram) 目录

实例演示R语言制作限制性立方条图

实例演示R语言制作限制性立方条图相关推荐

最新文章

热门文章