[R] Teasing out logrank differences *between* groups using survdiff or something else?

2009-09-15 Thread Bryan Hanson
R Folk: Please forgive what I'm sure is a fairly naïve question; I hope it's clear. A colleague and I have been doing a really simple one-off survival analysis, but this is an area with which we are not very familiar, we just happen to have gathered some data that needs this type of analysis.

Re: [R] Teasing out logrank differences *between* groups using survdiff or something else?

2009-09-15 Thread David Winsemius
On Sep 15, 2009, at 9:58 PM, Bryan Hanson wrote: R Folk: Please forgive what I'm sure is a fairly naïve question; I hope it's clear. A colleague and I have been doing a really simple one-off survival analysis, but this is an area with which we are not very familiar, we just happen to

Re: [R] Teasing out logrank differences *between* groups using survdiff or something else?

2009-09-15 Thread Thomas Lumley
I think you do in fact want to just run the analysis for the four groups you are interested in. The logrank chisquared test would then be of the hypothesis that these four groups have the same survival and censoring distributions, with the greatest power for detecting proportional-hazards

Re: [R] Teasing out logrank differences *between* groups using survdiff or something else?

2009-09-15 Thread Bryan Hanson
Thomas, thanks for your comments. We weren't entirely sure we we even framing the question right, your comments are encouraging. Here are our results: Call: coxph(formula = Surv(lifespan, status) ~ group, data = four) n= 573 coef exp(coef) se(coef) z Pr(|z|) groupT1