Displaying 1 result from an estimated 1 matches for "timebad".
2009 Jun 22
1
How to make try to catch warnings in logistic glm
...g(23.9)-0.06^2/2,sd=0.06))
?fatbad<-exp(rnorm(nbad, mean=log(27.4)-0.09^2/2, sd=0.09))
?fat<-c(fatgood, fatbad)
?thgood<-sample(thgood, ngood)
?thbad<-sample(thbad, nbad)
?btun<-c(thgood, thbad)
?btun<-as.factor(btun)
?timegood<-exp(rnorm(ngood, mean=log(443)-0.5^2/2, sd=0.5))
?timebad<- exp(rnorm(nbad, mean=log(555)-0.6^2/2, sd=0.6))
?time<-c(timegood, timebad)
?DF<-data.frame(outcome, ald, fat, hstop, btun, time)
?my.glm<-try(glm(outcome~ald + fat+ hstop + btun + time, data=DF,
family=binomial))
?test<-!is.null(attr(my.glm, "try-error"))
?if(test | !my...