search for: aml3

...u.dk/~bxc/SPE/library/, but built under 1.9. I recently tried to build it under 2.0, but it crashed for me and I was advised to wait till 2.1 before I had another go at it. A little further exploration of what goes on in survSplit gives: > data(aml) > aml$id <- letters[1:dim(aml)[1]] > aml3<-survSplit(aml,cut=c(5,10,50),end="time",start="start", + event="status",episode="i") > # Sort the rows sensibly > aml3 <- aml3[order(aml3$id,aml3$start),] > aml3$expand <- 1:(dim(aml3)[1]) > aml3[aml3$id=="k&quot...

...the modified example from the help page of survSplit. The original works, as expected. If, however, a Surv object is added to the data.frame, each record gets doubled. Is there some solution other than avoiding Surv objects in data.frames? Thanks, Heinz require(survival) ## from the help page aml3<-survSplit(aml,cut=c(5,10,50),end="time",start="start", event="status",episode="i") summary(aml) summary(aml3) coxph(Surv(time,status)~x,data=aml) ## the same coxph(Surv(start,time,status)~x,data=aml3) ## added to show doubling of records aml.so...

Hi, I am new to time-dependent Cox model to estimate time dependent hazard ratios. Let me use aml dataset from survival package: > aml3<-survSplit(aml2,cut=c(5,10,20),end="time",start="start", event="status",episode="i") If I want to esimate hazard ratio for each of the time intervals 0-5, 5-10, 10-20 and >=20, would the following calculate such HRs? > coxph(Surv(start,time,s...