Error message with survReg

[Tristan Lorino <tristan.lorino@lcpc.fr>]

Hi (again)...

I get an error message with survReg (S-Plus 7):

> survReg(Surv(cens04[, 2], cens04[, 3], type = "interval2") ~ cov04[, 2], dist 
> = "weibull", na.action = na.exclude)
Problem in survReg(Surv(cens04[, 2], cens04[, 3], type = "interval2")
~ cov04[, ..: Invalid survival times for this distribution

All my 5,316 times are good (i.e. cens04[,2] <= cens04[,3])... I tried
to find the error: the compilation is OK for the first 1,250 rows...
an error occures between line 1,250 and line 1,260:

survReg(Surv(cens04[1250:1260, 2], cens04[1250:1260, 3], type =
"interval2") ~ cov04[1250:1260, 2], dist = "weibull", na.action =
        na.exclude)
Problem in .C("S_coxph_wtest",: subroutine S_coxph_wtest: Missing values in 
argument 3

The problematic data are at the end of the mail... If anyone has an
idea, it would be great...

Thank you,
Tristan Lorino


> cens04[1250:1260,  ]
       IdS cg04 cd04 
1250 19176    3    6
1251 19177    3    9
1252 19178    3    6
1253 19179    3    6
1254 19180    3    6
1255 19181    3    6
1256 19182    6    6
1257 19183    6    6
1258 19184    6    6
1259 19185    3    6
1260 19186    3    6
> cov04[1250:1260,  ]
     Granulo MJAPL 
1250      10  1083
1251      10  1083
1252      10  1083
1253      10  1083
1254      10  1083
1255      10  1083
1256      10   943
1257      10   943
1258      10   943
1259      10  1083
1260      10  1083



-- 
Laboratoire Central des Ponts et Chaussées
[Division ESAR ? Section AGR]
Route de Bouaye BP 4129
44341 Bouguenais Cedex
France
Tél 33 (0)2 40 84 56 18
Fax 33 (0)2 40 84 59 92