From 41f386afb38667a7b9a0cd4e52f576c5e8946708 Mon Sep 17 00:00:00 2001
From: Maria Voigt <maria.voigt@idiv.de>
Date: Fri, 16 Jun 2017 18:38:01 +0200
Subject: [PATCH] fixing bug with bin_id
---
src/model_fitting/abundance_model.R | 6 +++---
1 file changed, 3 insertions(+), 3 deletions(-)
diff --git a/src/model_fitting/abundance_model.R b/src/model_fitting/abundance_model.R
index 9ccafdc..a2492f7 100644
--- a/src/model_fitting/abundance_model.R
+++ b/src/model_fitting/abundance_model.R
@@ -338,8 +338,8 @@ if (!is.na(exclude_rand)){
ids_to_exclude <- sample(predictors_obs$id,
size = nrow(predictors_obs)/100 * exclude_rand_perc,
replace = FALSE)
- predictors_excluded <- predictors_obs[predictors_obs$bin_id %in% ids_to_exclude, ]
- predictors_obs <- predictors_obs[!predictors_obs$bin_id %in% ids_to_exclude, ]
+ predictors_excluded <- predictors_obs[predictors_obs$id %in% ids_to_exclude, ]
+ predictors_obs <- predictors_obs[!predictors_obs$id %in% ids_to_exclude, ]
nr_excluded <- nrow(predictors_excluded)
}
@@ -516,7 +516,7 @@ results_res <- foreach(i = 1:nrow(all_model_terms),
predictors_excluded_pred <- predictors_excluded
predictors_excluded_pred$offset_term <- 0
prediction_transect_excluded <- predict.glm(res,
- newdata = predictors_excluded,
+ newdata = predictors_excluded_pred,
type = "response")
cross_lm = lm(log(predictors_excluded$ou_dens+ 1) ~
log(prediction_transect_excluded + 1))
--
GitLab