fit
fit
rstime<-rnorm(50,30,4)
rsentry<-rstime-runif(50,2,5)
rstatus<-as.numeric(rnorm(50)>0.1)
trt<-as.factor(as.numeric(rnorm(50)>0.15))
str<-as.factor(sample(1:3,50,replace=T))

    ### This is the one works
fit<-cph(Surv(rstime,rstatus)~trt+strat(str),x=T, y=T)
survest(fit, expand.grid(trt=c(0,1),str=as.factor(1:3)),
        times=c(20,25,30), conf.int=.95)

    ### This one doesn't
fit<-cph(Surv(rsentry,rstime,rstatus)~trt+strat(str),x=T, y=T)
survest(fit, expand.grid(trt=c(0,1),str=as.factor(1:3)),
        times=c(1,2,3), conf.int=.95)
q()
a
a
f
names(a)
setdiff(names(a), 'crap')
setdiff(names(a), 'TOTAL NONLINEAR')
a
A
a
reshape(a)
args(reshape)
?reshape
a
unname(a)
A
dotplot
A
args(dotplot)
?dotplot
s
a
    a <- plot(anova(f),
              what='partial R2', pl=FALSE, rm.totals=FALSE)

a
subset(a, select=-4)
args(subset)
args(subset.data.frame)
?subset.data.frame
subset(as.data.frame(a), select=-4)
as.data.frame(a)
row.names(a)
as.data.frame(a)
row.names(as.data.frame(a))
q()
