Complete_GC = GC_survival_data[, c("sname",
                                   "age_start_time","age_end_time",
                                   "status","age_sample",
                                   "grp","hydroyear",
                                   "mean_density",
                                   "density_squared",
                                   "prep_time",
                                   "rnk_alpha_new",
                                   "proprank_new",
                                   "extract_time",
                                   "state",
                                   "mean_tmax_30d",
                                   "season",
                                   "rain_anom_yr",
                                   "n_adult_mat_kin","gc",
                                   "adv_rain","adv_density","adv_sib",
                                   "adv_mom","adv_mom_dsi","adv_mom_rank")]



##Standard
Complete_GC[, c(
  "rain_anom_yr",
  "mean_density",
  "mean_tmax_30d",
  "density_squared",
  "n_adult_mat_kin",
  "proprank_new",
  "prep_time",
  "extract_time"
)] =
  scale(Complete_GC[, c(
    "rain_anom_yr",
    "mean_density",
    "mean_tmax_30d",
    "density_squared",
    "n_adult_mat_kin",
    "proprank_new",
    "prep_time",
    "extract_time"
  )])
complete_id=which(!apply(Complete_GC,1,FUN=function(x){any(is.na(x))}))
Complete_GC=Complete_GC[complete_id,]
Complete_GC$gc=log(Complete_GC$gc)


Complete_GC$adv_rain = as.integer(Complete_GC$adv_rain)
Complete_GC$adv_density = as.integer(Complete_GC$adv_density)
Complete_GC$adv_sib = as.integer(Complete_GC$adv_sib)
Complete_GC$adv_mom = as.integer(Complete_GC$adv_mom)
Complete_GC$adv_mom_dsi = as.integer(Complete_GC$adv_mom_dsi)
Complete_GC$adv_mom_rank = as.integer(Complete_GC$adv_mom_rank)



Complete_GC$adv_use=Complete_GC[,adverse_name[i]]


Complete_GC.id=aggregate(Complete_GC[,c("mean_density","rain_anom_yr","n_adult_mat_kin",
                                        "adv_use")],
                         by=list(Complete_GC$sname),mean)

Complete_GC.id=cbind(Complete_GC.id,aggregate(Complete_GC[,c("age_end_time","status")],
                                              by=list(Complete_GC$sname),max)[,c(2,3)])


Complete_GC = GC_survival_data[, c("sname",
                                   "age_start_time","age_end_time",
                                   "status","age_sample",
                                   "grp","hydroyear",
                                   "mean_density",
                                   "density_squared",
                                   "prep_time",
                                   "rnk_alpha_new",
                                   "proprank_new",
                                   "extract_time",
                                   "state",
                                   "mean_tmax_30d",
                                   "season",
                                   "rain_anom_yr",
                                   "n_adult_mat_kin","gc",adverse_name[i])]



##Standard
Complete_GC[, c(
  "rain_anom_yr",
  "mean_density",
  "mean_tmax_30d",
  "density_squared",
  "prep_time",
  "rnk_alpha_new",
  "proprank_new",
  "extract_time",
  "n_adult_mat_kin"
)] =
  scale(Complete_GC[, c(
    "rain_anom_yr",
    "mean_density",
    "mean_tmax_30d",
    "density_squared",
    "prep_time",
    "rnk_alpha_new",
    "proprank_new",
    "extract_time",
    "n_adult_mat_kin"
  )])
complete_id=which(!apply(Complete_GC,1,FUN=function(x){any(is.na(x))}))
Complete_GC=Complete_GC[complete_id,]
Complete_GC$gc=log(Complete_GC$gc)






Complete_GC$adv_use=Complete_GC[,adverse_name[i]]


Complete_GC.id=aggregate(Complete_GC[,c("mean_density","rain_anom_yr","n_adult_mat_kin",
                                        "adv_use")],
                         by=list(Complete_GC$sname),mean)

Complete_GC.id=cbind(Complete_GC.id,aggregate(Complete_GC[,c("age_end_time","status")],
                                              by=list(Complete_GC$sname),max)[,c(2,3)])