Displaying 1 result from an estimated 1 matches for "lower59".
2003 Feb 27
2
interval-censored data in survreg()
...plnorm(0)==0 and
pnorm(-Inf)==0 are well defined. Below is a short example to
reproduce the problem.
Does anyone know why survreg() must behave that way?
Is there an alternate solution to this problem?
Sincerely,
Jerome Asselin
library(survival)
data(ovarian)
newovarian <- ovarian
newovarian$lower59 <- newovarian$futime-59
newovarian$time59 <- Surv(newovarian$lower59,newovarian$futime,
event=rep(3,nrow(newovarian)),type="interval")
survreg(time59~ecog.ps+rx,data=newovarian,dist="lognormal")
#THIS DOES NOT WORK BECAUSE ONE OF THE LOWER BOUNDS IS ZERO
#Error in surv...