掌握语言中生存数据分析技术.pdfVIP

#源代码#怎样利用R”survival”package做分析

瑶轩

1.安装RPackage:survival(蓝色字体为代码)

install.packages(survival)

2.加载library(survival)

该包可以包含以下几个函数：

Createasurvivalsubject:Surv

KanMeierEstimator:survfit

Mantel-HaenzelTest:survdiff

CoxModel:coxph

3.以包自带数据PBC数据为例

data(pbc,package=survival)

Surv(pbc$time,pbc$status2)

surdays-with(pbc,Surv(pbc$time,pbc$status2))

Ps.+代表截尾数据

survfit(Surv(pbc$time,pbc$status2)~1)

plot(survfit(Surv(pbc$time,pbc$status2)~1))

diff-survdiff(Surv(pbc$time,pbc$status2)~pbc$trt)

pvalue-1-pchisq(diff$chisq,length(diff$n)-1)

fit-survfit(Surv(pbc$time,pbc$status2)~pbc$trt)

plot(fit,col=c(green,red),xlab=Time(days),ylab=Overallsurvival)

本文吗？欢迎加入云生信跟大家一起进步~

rsurvivalpackage

1.R软件包：（）

install.packages（“”）2。library（）

创建一个：幸存

KanMeier估计器：幸存

Mantel-Haenzel测试：SurvdiffCox型号：Coxph3.PBC数据（PBC，软件包=“”）

surv（pbc$time，pbc$status=2）

Surdays-with（pbc，surch（pbc$time，pbc$status=2）

PS。+

SurvIt（surs（pbc$time，pbc$status=2）〜1）

情节（surviT（cove（pbc$time，pbc$status=2）〜1）

diff-survdiff（surv（pbc$time，pbc$status=2）〜pbc$trt）

PVA