- 临床研究论文 COX 比例风险回归模型
- Fit Proportional Hazards Regression Model
- 如果如下内容有错误,请留言,谢谢,感激不尽!
1 写在前面
- 本文需要的R包:"survival" 、"tableone"、"broom"、"forestplot"和 "survminer".
- 本代码使用survival自带的lung数据集,lung数据集中出现的变量解释:
![caa134a0a1f6184beee0837400a074da.png](https://i-blog.csdnimg.cn/blog_migrate/52bf6da48e5e360aa4be213d5ba90bd4.jpeg)
2 安装拓展包
install.packages("survival")install.packages("survminer")install.packages("forestplot")install.packages("tableone")install.packages("broom")# 不需要的包可以不安装
3 调用拓展包
library(survival) # 生存分析需要library(survminer) # 作图使用ggplot画library(forestplot) # 画森林图library(tableone) library(broom) # 后面tidy函数需要data("lung") head(lung) # 预览 lung 数据集信息
4 Kaplan-Meier曲线
4.1 统计分析
fit survdiff(Surv(time,status) ~ sex,data = lung, rho = 0) # 比较分析两组的生存时间分布是否不同# rho = 0 表示使用long-rank检验或者Mantel-Haenszel 检验)