[R] Survfit plots in trellis graphics

[Marco Chiarandini]

Dear all,

is there a way to produce survfit plots in a trellis environment?

I am trying this:

print(Ecdf(~time | size*type, groups=alg,data=B,subscripts=TRUE,
            panel=function(x,groups,subscripts)
            {
              t <- survfit(Surv(time[subscripts],event[subscripts])~groups[subscripts],data=B,conf.type="none")
              panel.xyplot(t[1]$time,1-t[1]$ssurv,type="s",lty=2)
              panel.xyplot(t[2]$time,1-t[2]$ssurv,type="s",lty=2)
            }
)



but things get messed up if I try log transformations and to plot
confidence intervals.


- Marco.


--
Marco Chiarandini                     http://www.imada.sdu.dk/~marco
Department of Mathematics	      Email: marco AT imada.sdu.dk
and Computer Science,		      Phone: +45 6550 4031
University of Southern Denmark        Fax: +45 6593 2691