r语言找不到cochrane函数_R语言可视化（三十九）：森林图绘制

39. 森林图绘制

清除当前环境中的变量

rm(list=ls())

设置工作目录

setwd("C:/Users/Dell/Desktop/R_Plots/39forest/")

使用survminer包中的ggforest函数绘制森林图

require("survival")

## Loading required package: survival

library(survminer)

## Warning: package 'survminer' was built under R version 3.6.3

## Loading required package: ggplot2

## Warning: package 'ggplot2' was built under R version 3.6.3

## Loading required package: ggpubr

## Loading required package: magrittr

# 查看内置示例数据

head(colon)

## id study rx sex age obstruct perfor adhere nodes status differ

## 1 1 1 Lev+5FU 1 43 0 0 0 5 1 2

## 2 1 1 Lev+5FU 1 43 0 0 0 5 1 2

## 3 2 1 Lev+5FU 1 63 0 0 0 1 0 2

## 4 2 1 Lev+5FU 1 63 0 0 0 1 0 2

## 5 3 1 Obs 0 71 0 0 1 7 1 2

## 6 3 1 Obs 0 71 0 0 1 7 1 2

## extent surg node4 time etype

## 1 3 0 1 1521 2

## 2 3 0 1 968 1

## 3 3 0 0 3087 2

## 4 3 0 0 3087 1

## 5 2 0 1 963 2

## 6 2 0 1 542 1

# 构建COX回归比例风险模型

model

data = colon )

# 查看cox回归模型结果

model

## Call:

## coxph(formula = Surv(time, status) ~ sex + rx + adhere, data = colon)

##

## coef exp(coef) se(coef) z p

## sex -0.04615 0.95490 0.06609 -0.698 0.484994

## rxLev -0.02724 0.97313 0.07690 -0.354 0.723211

## rxLev+5FU -0.43723 0.64582 0.08395 -5.208 1.91e-07

## adhere 0.29355 1.34118 0.08696 3.376 0.000736

##

## Likelihood ratio test=46.51 on 4 df, p=1.925e-09

## n= 1858, number of events= 920

# 使用ggforest()函数绘制基础森林图

ggforest(model)