Hi Rlisters:

I'm doing survival analysis comparing 8 populations. I've used the survdiff function for comparing survival curves and the result I obtained says that the survivorship curves are different (p=1.55e-15). Now I want to know if all the curves differ or if there are some of them that don't. What function should I use for doing this comparison? I thought about running the test again with every possible comparison between the populations, but this seems like the long way of doing this. Sorry if this is a silly question but I'm not too familiar using R.
Thanks in advance,

Adriana Humanes

Sobrevivencia <- read.table(file="Rec_surv_adri_all_tiles.txt", header=TRUE)
summary(Sobrevivencia)
attach(Sobrevivencia)
library(survival)
Sobrevivencia.surv <- Surv(age, status)
survdiff(Sobrevivencia.surv~reef)
adri.fit=survfit(Sobrevivencia.surv~reef)
survdiff(Sobrevivencia.surv~reef)    
summary(adri.fit)                     
plot(adri.fit, lty=1:8, col=1:8)
legend("bottomleft", c("Dunk", "Fitzroy", "Hannah", "Hay", "High", "MacDonalds", "Russell", "Wilkie"), lty=1:8, col=1:8) #leyenda

	[[alternative HTML version deleted]]

