R/modtiscon_mcs_fun.R

In SCCWRP/SQObioaccumulation: Run bioaccumulation model for SQO human health objectives

#' mcs function for modelled tissue concentration
#'
#' @param nsim number of simulations
#' @param meanse mean and se values from user input for each guild species and contaminant class
#' @param propseaf proportion of seafood diet
#'
#' @export
modtiscon_mcs_fun <- function(nsim, meanse, propseaf){
  
  # simulated tissue concentrations across guild species, all sims
  sims <- meanse %>%  
    group_by(MCSvar, contam) %>% 
    mutate(
      ests = purrr::pmap(list(nsim, X, SD), genlognorm_fun)
    ) %>% 
    dplyr::select(-SD, -X) %>% 
    unnest(ests)

  # weighted tissue concentrations across guilds for each contam, all sims
  out <- sims %>% 
    dplyr::group_by(i, contam) %>% 
    arrange(i, contam, MCSvar) %>% 
    mutate(
      sims = case_when(
        sum(is.na(sims)) == 9 ~ NaN,
        is.na(sims) ~ 0, 
        T ~ sims
      )
    ) %>% 
    summarise(
      wgtave = sims %*% propseaf
    )
 
  
  return(out)
  
}