This should solve your problem:
mod <- survreg((Surv(as.numeric(Time), event=Status)) ~ Prison+Dose+Clinic,
data = meth,
dist = "lognormal")
as you reported to me that:
names(meth)
# [1] "ID" "Clinic" "Status" "Time" "Prison" "Dose"
Note, it is Time
, not time
; Also, it is Status
, not status
. In fact, all variables names start with a capital letter!
5
solved Survival analysis [closed]