Survcomp包,y z变量

时间:2016-01-03 18:30:25

标签: survival-analysis

我正在运行survcomp包并想知道y和z值。我有多个临床数据:

> colnames(ClinicalDataHep)
 [1] "follow_upTime"                                              
 [2] "RecurrenceTime"                                             
 [3] "Age"                                                        
 [4] "OS"                                                         
 [5] "Survival_dead0_alive1"                                      
 [6] "Tumour_size"                                                
 [7] "HVB_preop"                                                  
 [8] "HCV_preop"                                                  
 [9] "HBD_preop"                                                  
[10] "Cirrhosis_preop"                                            
[11] "Status:_no_recurrence-0._recurrence-1_"                     
[12] "Surgery:_resection-1._tx-2;_rfa-3;_resection+rfa-4;tx+rfa-5"
[13] "new_time"                                                   
[14] "new_death"                                                  
[15] "death_event" 

使用Overall Survival作为y-variabledead/alive作为z variable是否正确?

cindexall.Hep.serum <- as.data.frame(t(apply(X=matrix_cpm, MARGIN=1, function(x, y, z) { 
  tt <- concordance.index(x=x, surv.time=y, surv.event=z, method="noether", na.rm=TRUE);  
  return(c("cindex"=tt$c.index, "cindex.se"=tt$se, "lower"=tt$lower, "upper"=tt$upper,"p.value"=tt$p.value)); }, 
  y=ClinicalData$OS,  z=ClinicalData$Survival_dead0_alive1)))

0 个答案:

没有答案