changing the maxit in abundance_model

src/model_fitting/abundance_model.R

@@ -319,7 +319,7 @@ model_terms <- names(glm.nb(as.formula(paste("nr_nests~", paste(m_terms,
                                              "+ offset(offset_term)",
                                              sep = "")),
                               data = predictors_obs,
-                            control = glm.control(maxit = 500))$coefficients)
+                            control = glm.control(maxit = 1000))$coefficients)
 # prediction estimates
 intercept <- rep(1, nrow(predictors_obs))
 predictor_estimates <- cbind( intercept,
@@ -343,7 +343,7 @@ model <- as.formula(
   paste("nr_nests ~", full_model, "+ offset(offset_term)"))
 
 res_full <- glm.nb(model, data = predictors_obs,
-                   control = glm.control(maxit = 500))
+                   control = glm.control(maxit = 1000))
 
 # HERE I CAN NOW USE THE OTHER FUNCTION
 dfbeta_frame <- data.frame(slope=res_full$coefficients, res_full$coefficients+
@@ -386,7 +386,7 @@ results_res <- foreach(i = 1:nrow(all_model_terms),
 	    paste(m_terms[all_model_terms[i, ] == 1], collapse = "+"),
             "+ offset(offset_term)"))
             res <- glm.nb(model, data = predictors_obs,
-                          control = glm.control(maxit = 500))
+                          control = glm.control(maxit = 1000))
 
 
 # model