1 / 7
文档名称:

用R作生存分析.doc

格式:doc   大小:74KB   页数:7页
下载后只包含 1 个 DOC 格式的文档,没有任何的图纸或源代码,查看文件列表

如果您已付费下载过本站文档,您可以点这里二次下载

分享

预览

用R作生存分析.doc

上传人:pppccc8 2019/1/25 文件大小:74 KB

下载得到文件列表

用R作生存分析.doc

文档介绍

文档介绍:()library(survival)library(help=survival)data(aml)ami#查看你的系统中所有己装的libraries#加载一个library,也口J使用鼠标来做#seethelistofavailablefunctionsanddatasets.#加载数据ami#seethedata#《生存分析》需要的两个librarieslibrary(survival)#包含我们教材中所有的数据library(KMsurv) #包"#若想使数据中的变量能够在R中作为一个变量使用,useattach():attach(aids)infectdetach(aids)—.SurvivalObjects#由函数Surv产生一个survivalobjectSurv(time,time2,event,type)time:survivaltimetime2:endingtimeifitisintervalcensoredevent:,O=rightcensored,l=eventattime,2=leftcensored,3=:indicatingwhetheritisrightcensored,#假设一个观测值为[2,3],则R语句为Surv(time=2,time2=3,event=3,type="interval")#右删失数据attach(aml)Surv(time,status)detach(aml)#左截断右删失数据data(psych)attach(psych)<-Surv(age,age+time,death)[1](51,52](58,59](55,57](28,50](21,51+](19,47](25,57][22](29,63+](35,65+](32,67](36,76](32,71+]detach(psych)三TheKaplan-MeierEstimates#计算生存曲线:survfit,有三个重要的参数:formula,,(formula,=,="log")"plai #〃log〃 #对H(t)"logJog" #log-transformedinourbook#在新版的R里,这个"是必要的。#在新版的R里,这个T是必要的。myfit=survfit(Surv(time,status)~l,data=kidney)attach(kidney)myfit=survfit(Surv(time,status)~1)a=summary(survfit(Surv(time,status