Thanks Marc;
That provides the answer that I had missed (and glaringly obvious
now). The plot on p 128 of Therneau and Grambsch is correctly labeled
with -Log(-Log(Survival). I was expecting Log(-Log(S)) based on my
other references.
--
David Winsemius
On Mar 1, 2009, at 12:42 PM, Marc
Using Frank's survplot() which allows for user definable functions as
the 'fun' argument, one could also do something like this:
library(Design)
neg.ll <- function(x) -log(-log(x))
survplot(fit, fun = neg.ll, conf = "none")
and the subsequent examples seem to work as well, so if the output is
Here is a modification of some simplistic code that I had sent Bob
offlist for the gastric data specifically.
I created a function and made it a bit more generic, for multiple
'strata', though there is no real error checking, etc. which would be
needed to make it more robust for production along
I think what you want may be produced by this code for InvNormal(S) vs
log(time):
survplot(fit, fun=qnorm, logT=T, conf = "none")
That is not what you describe, however.
I am worried about the plot on the page you cite, because it is not
similar to other log(-log(S)) (complementary log-log)
I am hoping for some advice regarding how to obtain a log-log
survival plot that is not in the inverse. On page 128 of Modelling
survival data by Therneau & Grambsch there is the an example of the
type of desired plot, with a log of the survival curve by
years. Marc Schwartz has provided me w
5 matches
Mail list logo
