adverse.name <-
c(
  "adv_rain",
  "adv_density",
  "adv_sib",
  "adv_mom",
  "adv_mom_sci",
  "adv_mom_rank",
  "adv_01",
  "adv_12",
  "adv_cumulative_sci"
)

adverse.index <-
c(
  "Rain Adversity",
  "Sibling Adversity",
  "Density Adversity",
  "Maternal Death",
  "Maternal SCI",
  "Maternal Rank",
  "Cumulative 0/1 Compared",
  "Cumulative 1/2 Compared",
  "Cumulative Adversity"
)


covariate.index = c(
  "mean_density",
  "density_squared",
  "mean_prop_rank",
  "n_adult_mat_kin",
  "rain_anom_yr"
  # "hyb_score"
)
if(hybrid.score){
  covariate.index = c(covariate.index,"hyb_score")
}




if(m==1){
covariate.index.dsi = c(
  "mean_density",
  "n_adult_mat_kin",
  "mean_prop_rank",
  "percent_days_with_infant"
  #"hyb_score"
)
}
if(m==2){
covariate.index.dsi = c(
  "mean_density",
  "mean_prop_rank",
 # "hyb_score"
  "percent_days_cycling"
)
}

if(m==3){
covariate.index.dsi = c(
  # "state",
  "mean_density",
  "density_squared",
  "mean_tmax_30d",
  # "season",
  "rain_anom_yr",
  "extract_time",
  "prep_time",
  "rnk_alpha", 
  "mean_prop"
  #"hyb_score"
)
}

if(m==4){
covariate.index.dsi = c(
  "mean_density"
  #"hyb_score"
)
}

if(hybrid.score){
  covariate.index.dsi = c(covariate.index.dsi,"hyb_score")
}


mediator.name=c("DSI_F","DSI_M","gc","mean_prop_rank")
mediator.index=c("DSI_F","DSI_M","GC","mean_prop_rank")

# mediation.index <- c("DSI_F", "DSI_M", "proprank")
# response.index <- c("gc")
group_id=c("grp")


covariate.collect.index=union(covariate.index.dsi,covariate.index)