R/nw_setup.R

In EvoNetHIV/Test_and_Treat: Evonet: Integrated Bio-Social Models for HIV Epidemiology

#' @export

nw_setup <- function(params){

nw <- setup_initialize_network(params)

#--------------------------------------------------------------

#estimate initial network (create argument list, then call fxn)
netest_arg_list <- list(
  nw            =  nw,
  formation     =  as.formula(params$nw_form_terms),
  target.stats  =  params$target_stats,
  coef.form     =  params$nw_coef_form,
  constraints   =  as.formula(params$nw_constraints),
  verbose       =  FALSE,
  coef.diss     =  dissolution_coefs( dissolution =  as.formula(params$dissolution),
                                      duration    =  params$relation_dur,
                                      d.rate      =  params$d_rate) )

sim <- do.call(EpiModel::netest, netest_arg_list)
if(params$hyak_par==T){
  save(sim,file="estimated_nw.RData")
}
return(sim)
}