The list of dates returned in case of detection of NA in sim was not …

…correct
SticsRPacks · Apr 18, 2024 · ad395ff · ad395ff
1 parent 4603053
commit ad395ff
Showing 1 changed file with 4 additions and 2 deletions.
diff --git a/R/obsSim_consistency.R b/R/obsSim_consistency.R
@@ -236,8 +236,10 @@ is_sim_inf_or_na <- function(sim_list, obs_list, param_values) {
  var_list <- lapply(names(sim_list[[sit]]), function(x) {
  if (any(is.infinite(sim_list[[sit]][!is.na(obs_list[[sit]][,x]),][[x]])) ||
  any(is.na(sim_list[[sit]][!is.na(obs_list[[sit]][,x]),][[x]]))) {
- return(list(sim_list[[sit]]$Date[is.infinite(sim_list[[sit]][!is.na(obs_list[[sit]][,x]),][[x]]) |
- is.na(sim_list[[sit]][!is.na(obs_list[[sit]][,x]),][[x]])]))
+ return(list(sim_list[[sit]]$Date[(is.infinite(sim_list[[sit]][[x]] &
+ !is.na(obs_list[[sit]][,x]))) |
+ (is.na(sim_list[[sit]][[x]] &
+ !is.na(obs_list[[sit]][,x])))]))
  } else {
  return(NULL)
  }