Fix test on presence of NA in simulations

SticsRPacks · Dec 15, 2023 · 9ac1b88 · 9ac1b88
1 parent c4e166d
commit 9ac1b88
Showing 1 changed file with 4 additions and 4 deletions.
diff --git a/R/main_crit.R b/R/main_crit.R
@@ -397,10 +397,10 @@ main_crit <- function(param_values, crit_options) {
  # check presence of Inf/NA in simulated results where obs is not NA
  for (sit in names(obs_sim_list$sim_list)) {
  var_list <- lapply(names(obs_sim_list$sim_list[[sit]]), function(x) {
- if (any(is.infinite(obs_sim_list$sim_list[[sit]][!is.na(obs_sim_list$obs_list[[sit]][,x]),x])) ||
- any(is.na(obs_sim_list$sim_list[[sit]][!is.na(obs_sim_list$obs_list[[sit]][,x]),x]))) {
- return(list(obs_sim_list$sim_list[[sit]]$Date[is.infinite(obs_sim_list$sim_list[[sit]][!is.na(obs_sim_list$obs_list[[sit]][,x]),x]) |
- is.na(obs_sim_list$sim_list[[sit]][!is.na(obs_sim_list$obs_list[[sit]][,x]),x])]))
+ if (any(is.infinite(obs_sim_list$sim_list[[sit]][!is.na(obs_sim_list$obs_list[[sit]][,x]),][[x]])) ||
+ any(is.na(obs_sim_list$sim_list[[sit]][!is.na(obs_sim_list$obs_list[[sit]][,x]),][[x]]))) {
+ return(list(obs_sim_list$sim_list[[sit]]$Date[is.infinite(obs_sim_list$sim_list[[sit]][!is.na(obs_sim_list$obs_list[[sit]][,x]),][[x]]) |
+ is.na(obs_sim_list$sim_list[[sit]][!is.na(obs_sim_list$obs_list[[sit]][,x]),][[x]])]))
  } else {
  return(NULL)
  }