Adrian,

You can take a look at the R package wild1.

Best,

Manuel Spínola

2012/9/6 AdRiAnA HuMaNeS <adrihumanes@yahoo.com>

> 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]]
>
>
> _______________________________________________
> R-sig-ecology mailing list
> R-sig-ecology@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-ecology
>
>


-- 
*Manuel Spínola, Ph.D.*
Instituto Internacional en Conservación y Manejo de Vida Silvestre
Universidad Nacional
Apartado 1350-3000
Heredia
COSTA RICA
mspinola@una.ac.cr
mspinola10@gmail.com
Teléfono: (506) 2277-3598
Fax: (506) 2237-7036
Personal website: Lobito de río <https://sites.google.com/site/lobitoderio/>
Institutional website: ICOMVIS <http://www.icomvis.una.ac.cr/>

	[[alternative HTML version deleted]]

