Hi!

I am having a few problems with coxph function, I had the same problem with
the use of survfit. Here it is:

when calling 'M22<-coxphw(Surv(V1,V2,V4)~V5, data=XTDV, id=XTDV$V3, *
prentice*= ~V5, robust=TRUE, censcorr=TRUE)'

knowning that my data is:
> XTDV[1:10,1:3]
   V1       V2   V3 V4        V5
1   0 36.39706   36  1  343.2224
2   0 36.39706  943  0  161.5931
3   0 36.39706   45  0  790.3714

I get this message: "Error in w[i] <- w[i] * w.sorted[fit$time ==
obj$resp[i, 2]] : replacement has length zero". When I perform a debug I
find that the problem is that the time returned by obj$resp[i,2] is for a
certain i different from fit$time so all the vector is false....but in fact
the difference is minimal and it doesn't do that with precedent values...

Browse[2]> w.sorted[fit$time==obj$resp[2051,2]]
numeric(0)
Browse[2]> obj$resp[2051,2]
[1] 49.53238
Browse[2]> fit$time[23]
[1] 49.53238
Browse[2]> obj$resp[2051,2]-fit$time[23]
[1] -7.105427e-15

maybe it's because the time I use is not integer? Like I said, I had the
same problem using results from survfit but I fixed it by rounding values
before comparing them. Is there a way to change the function coxphw?

Thank you.

	[[alternative HTML version deleted]]

