diff --git a/src/model_fitting/abundance_model.R b/src/model_fitting/abundance_model.R index 63a5a7f0cdd742eeb85332ed96db89954acde80e..638f606ee07566056dafad6618b47b3892a0f746 100644 --- a/src/model_fitting/abundance_model.R +++ b/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