[EMAIL PROTECTED] writes: > The problem is caused by exp(x) returning Inf when x is large. > > One way of fixing the problem [there may be better ways] would be to include the > line > > out[x > 709] <- 1 > > in inv.logit()
Yes, or replace the entire code with something like inv.logit2 <- function(x) ifelse(x > 0, 1 - 1/(1+exp(x)), 1/(1+exp(-x))) -- O__ ---- Peter Dalgaard Blegdamsvej 3 c/ /'_ --- Dept. of Biostatistics 2200 Cph. N (*) \(*) -- University of Copenhagen Denmark Ph: (+45) 35327918 ~~~~~~~~~~ - ([EMAIL PROTECTED]) FAX: (+45) 35327907 ______________________________________________ [EMAIL PROTECTED] mailing list http://www.stat.math.ethz.ch/mailman/listinfo/r-devel
