使用R软件分析竞争风险模型简明攻略

©1994-2010ChinaAcademicJournalElectronicPublishingHouse.Allrightsreserved.(100050)(GVHD),,((right2censorddata)),,,,1,SAS,SPSS,Rcmprsk,,,R,R2,335,ALL17,AML18SPSSbmt.sav,11dis0ALL,1AMLftimestatus0,1(TRM),2,,TRM,(R-2.6.1),R,R,R3R,,Rcmprsk,:bmt=read.spss(bmt.sav,to.data.frame=TRUE)attach(bmt)bmtR,,,(factor,),:dis=factor(dis,levels=c(0,1),labels=c(ALL,AML)),dis(0,1),(0=ALL,1=AML)21(CIF)Gray4,splines,survivalcmprsk,:x2cuminc(ftime,status,dis);xcumincGray,,x,(2)2,Gray(,Gray),,CIF,TRM;2220,40,60,,1550CIF,:836200812256©1994-2010ChinaAcademicJournalElectronicPublishingHouse.Allrightsreserved.(x,c(1,5,50)),CIF,22cumincTests:(Gray)stat()pv(P)df()1(TRM)1.3016350.25391487112()7.0822780.0077849981EstimatesandVariances:()est()204060ALL1(TRM)0.17647060.17647060.1764706AML1(TRM)0.36818580.36818580.3681858ALL2()0.58823530.70588240.7058824AML2()0.20570750.20570750.2057075var()204060ALL1(TRM)0.0093441690.0093441690.009344169AML1(TRM)0.0160110460.0160110460.016011046ALL2()0.0160464780.0145849060.014584906AML2()0.0133583790.0133583790.01335837931CIFR,:plot(x,xlab=,ylab=),CIF(1)1RCIFRCIF,,,Word,!11R,21,,R,R,231,,R,CIF,:plot(xAML2time,xAML2est,type=s,lty=1,ylim=c(0.0,1.0),xlab=(),ylab=,lwd=3)plot.new()lines(xALL2time,xALL2est,type=s,lty=2,lwd=3)legend(0,0.9,c(AML,ALL),lty=1:2,lwd=3)AML2();type,s;lty,1,2;ylim;lwd,123ALL2()4,(0,019)22CIF41RGray,crr(),,,R,,R11DavidHA,MoeschbergerML.TheTheoryofCompetingRisks.Lon2don:CharlesGriffin,1978.21:,2007,41257.41GrayRJ.AclassofK2sampletestsforcomparingthecumulativeinci2denceofacompetingrisk.AnnStat,1988,16(3):114121145.936ChineseJournalofHealthStatistics,Dec2008,Vol.25,No.6