Skip to content
Prev 310831 / 398503 Next

Jackknife in Logistic Regression

untested, but something like this should get you what you want:

no.it <- 5
out <- vector("list", length=no.it)

      for(i in 1:no.it){
mydata2 <- mydata[ sample(1:nrow(mydata),  76000/no.it) ,]

out[[i]]  <-  coef(  
glm(f_ocur~altitud+UTM_X+UTM_Y+j_sin+j_cos+temp_res+pp+offset(log(1/off)),
data=mydata2, family='binomial')   )  

      }

do.call(c, out) 





Lucas wrote

            
--
View this message in context: http://r.789695.n4.nabble.com/Jackknife-in-Logistic-Regression-tp4649520p4649563.html
Sent from the R help mailing list archive at Nabble.com.