|
[ÇóÖú]
ÈçºÎʹÓÃRÓïÑÔcmprsk°ü×öÒÔageΪtime scaleµÄ¾ºÕù·çÏÕÉú´æ·ÖÎö
¸÷λ´óÉñ~
ÎÒÃÇÏ£ÍûÓÃ5ÄêµÄÉú´æ·ÖÎöÊý¾ÝÀ´½øÐÐlife time ·çÏÕµÄÔ¤²â£¬ÐèÒªÇóstrokeµÄ·çÏÕº¯Êý£¨ÒÔdeathΪ¾ºÕù·çÏÕʼþ£©£¬²¢ÄêÁä×÷Ϊtime scale£¨³£¹æµÄ·ÖÎöʹÓÃÉú´æʱ¼äÀ´·ÖÎö£©¡£
ÎÒÃÇÏ£Íû½«³ö×éÄêÁäfage×÷Ϊ·¢ÉúʼþµÄʱ¼ä£¬½«Èë×éÄêÁäage×÷Ϊ½øÈëÄ£Ð͵Äʱ¼ä£¬²»ÖªµÀÔÚ²ÅcmprskÀïÃæÔõôʵÏÖ£¿
use (a)
mcrr <- crr(fage,fstatus, as.factor(treatment group))
ÔÚcmprskÀïÃæÖ»ÄÜÉèÖóö×éÄêÁäΪftime£¬²»ÄÜ¿¼ÂÇÈë×éÄêÁäageµÄÉèÖã¿
mcoxph <- coxph(Surv(age, fage, fstatus==1) ~ as.factor(treatment group) , data=a)
ÓÃcoxph º¯ÊýËäÈ»Äܹ»ÉèÖÃÈë×éÄêÁä¡¢³ö×éÄêÁäΪfstartºÍfstop£¬µ«Ä£ÐÍÖ»¿¼ÂÇÁËfstatus==1¼´strokeµÄ·çÏÕ£¬ÆäËûµÄ1ºÍ0±àÂ붼ĬÈÏΪɾʧʼþ£¬Ëã³öÀ´µÄ²»ÊǾºÕù·çÏÕ¡£
ÎÒµÄÎÊÌâÊÇ£º
ÏëÎÊÎʸ÷λ´óÉñ£¬ÓÃcrrº¯ÊýµÄʱºò£¬ÓÐûÓа취ÏñcoxphÒ»Ñù£¬ÔÚÄ£ÐÍÀïÃ潫Èë×éÄêÁäageÉèÖÃΪ½øÈëÄ£Ð͵Äʱ¼ä£¿
»òÕß˵ÔÚcoxphº¯Êý»ù´¡ÉÏÄܲ»ÄܽøÐоºÕùʼþµÄÉèÖã¬À´×ö¾ºÕù·çÏÕ·ÖÎö£¿
»òÕß˵survialÓÐÆäËûµÄ°üÄܹ»ÊµÏÖÕâ¸öÄ¿µÄ£¿
»ðºü½Øͼ_2018-01-09T12-33-25.009Z.png |
|