brm_pathfinder_inits.R

library(cmdstanr)
library(brms)

brm_pathfinder_inits <- function(..., psis_resample=TRUE, max_lbfgs_iters = 2500, history_size = 100){
  brm(empty = T, ...) -> 
    brm_empty
  
  cmdstanr::write_stan_file(brm_empty$model) %>% 
    cmdstan_model(cpp_options = list(stan_threads = TRUE)) -> 
    model_cmdstan
  
  model_cmdstan$pathfinder(data = make_standata(...), 
                           init = 2, 
                           history_size = history_size, 
                           max_lbfgs_iters = max_lbfgs_iters,
                           num_threads = list(...)$cores, psis_resample=psis_resample,
                           num_paths = list(...)$chains) ->
    pathfinder_inits
  
  
  brm(init = pathfinder_inits, ...) ->
    output_pathfinder_inits
  
  return(output_pathfinder_inits)
}