生存时间的中位值不等于中位生存时间（一段错误的代码）

一段错误的代码记录

以下是自己分析数据过程中的一段错误代码，写完才发现自己算错了，但还是保留下来。

有没有义务给demo?

顺便说一下，我没有义务给大家测试数据，能看懂的就看看，对大家有帮助我开心。但有时候我就是会懒，但是请注意我这样没有错，这是自由，我相信你会认可这样的自由。

### group by subtype

test<-aggregate(x=mydata$Survival_months,by=list(mydata$Subtype),FUN=median,na.rm=T) 
test2<-aggregate(x=mydata$Survival_months,by=list(mydata$Subtype),FUN=quantile,na.rm=T)
test2<-cbind(test2,test2$x[,c(2,4)])## 
test2<-test2[,-c(1:2)]
test2
surv_x1<-cbind(test,test2)
colnames(surv_x1)<-c("Subtype","median","lower","upper")
surv_x1$IQR<-paste(surv_x1$lower,surv_x1$upper,sep="-")
surv_x1$median_IQR<-paste0(surv_x1$median,"(",surv_x1$IQR,")")


### Function to group
med_IQR<-function(x,data=mydata){
  
  test<-aggregate(x=data$Survival_months,by=list(data[[x]]),FUN=median,na.rm=T) 
  test2<-aggregate(x=data$Survival_months,by=list(data[[x]]),FUN=quantile,na.rm=T)
  test2<-cbind(test2,test2$x[,c(2,4)])## 
  test2<-test2[,-c(1:2)]
  test2
  surv_x1<-cbind(test,test2)
  colnames(surv_x1)<-c("Group","median","lower","upper")
  surv_x1$IQR<-paste(surv_x1$lower,surv_x1$upper,sep="-")
  surv_x1$median_IQR<-paste0(surv_x1$median,"(",surv_x1$IQR,")")
  return(surv_x1)
}
## test function
med_IQR("Subtype",data=mydata)

### 批量计算
vars=c("Age_group","Sex","Race","Marital_status","Insurance","Laterality","Primary_site","Surgery","Radiation","Chemotherapy","Histology","M_count","AJCC_T","AJCC_N","AJCC_Stage","Grade","Subtype")
med_IQR_res<-lapply(vars,FUN=med_IQR)
str(med_IQR_res)
med_IQR_res[[1]]
### 
med_IQR_res<-do.call(rbind,med_IQR_res)
head(med_IQR_res,10)
write.csv(med_IQR_res,file = "All_patients_med_IQR_res_data.csv")