[R] plot.survfit and cloglog

[gb]

With R-1.0.1, Linux RH6.2 and survival5, I try

> x2.fit <- survfit(Surv(stopp, event) ~ x2)
> plot(x2.fit, fun = "cloglog")
Error in rep(2, n2 - 1) : invalid number of copies in "rep"
In addition: Warning message: 
1 x value <= 0 omitted from logarithmic plot in: xy.coords(x, y, xlabel,
ylabel, log) 

(an empty plot area pops up), also

> print(x2.fit)
Call: survfit(formula = Surv(stopp, event) ~ x2)

      n events mean se(mean) median 0.95LCL 0.95UCL
x2=0 52     23 3.19    0.514  2.277   0.839     Inf
x2=1 48     40 0.69    0.173  0.282   0.169    0.66
Warning message: 
no finite arguments to min/max; returning extreme. in: min(..., na.rm =
na.rm) 

seem to indicate something strange. However, 'fun="cumhaz"' works as
expected. Is this a bug, or am I doing something stupid (again)? 

I tried to replicate this with a small data set:
> t1
 [1] 1.233 0.391 0.219 0.420 0.710 1.658 1.487 2.888 0.175 0.045
> ev
 [1] 1 1 1 1 1 1 1 1 1 1

> fit <- survfit(Surv(t1, ev) ~ 1)
> plot(fit, fun = "cloglog")
Error in rep(2, n2 - 1) : invalid number of copies in "rep"
In addition: Warning message: 
1 x value <= 0 omitted from logarithmic plot in: xy.coords(x, y, xlabel,
ylabel, log) 

However, this time a plot (apparently correct) shows up.

Göran

-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-help-request at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._