Hi,

I am using coxph from the survival package to fit a large model (100,000 observations, ~35 covariates) using both ridge regression (on binary covariates) and penalized splines (for continuous covariates).

In fitting, I get a strange error:

Error in if (abs((y[nx] - target)/(y[nx - 1] - target)) > 0.6) doing.well <- FALSE else doing.well <- TRUE :
  missing value where TRUE/FALSE needed

Unfortunately, I can't reproduce this error without handing over my entire dataset, but I thought it would be worth checking if anyone had any insight. I should note that the outcome that I'm using has almost everyone having an event (~98,000 events out of 100,000). I have fit other models like this with no problem, but on one particular dataset it fails.

Thanks!

Jessica Myers
Instructor in Medicine
Brigham and Women's Hospital


The information in this e-mail is intended only for the ...{{dropped:7}}

______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

Reply via email to