计算多个分组间的的p值，可以各种推广

allclu<-read.csv("clinical_cluster_west_all.csv")
library(survival)
fit<-survfit(Surv(survival.time, as.numeric(survival.event))~cluster,data=allclu)
survdiff(Surv(survival.time, survival.event) ~ cluster,data=allclu)




c.ch <- as.character(allclu$cluster)
u.ch <- unique(c.ch)
combn <- combn(u.ch,2)
N <- ncol(combn)


p.all <- NULL
chisq.all <- NULL
for(i in 1:N)
{
  index <- c.ch %in% combn[,i]
  temp <- allclu[index,]
  t <- survdiff(Surv(survival.time, survival.event) ~ cluster,data=temp)
  
  
  ### get the chisq
  chisq <- t$chisq
  
  ### get the p.value
  p <- capture.output({print(t,digits=8)})
  p <- p[length(p)]
  start <- regexpr("p=",p)+2
  p <- as.numeric(substr(p,start,nchar(p)))
  
  p.all <- c(p.all,p)
  chisq.all <- c(chisq.all,chisq)
}




result <- data.frame(combn1=combn[1,],combn2=combn[2,],chisq=chisq.all,p=p.all)