inst/scripts/choose_best_model.R

In joshwlambert/relaxedDAISIE: Reproduce Relaxed-rate DAISIE Paper

# script to collate all the DAISIE maximum likelihood parameter estimates for
# each island

# Canaries

canaries_cr_dd <- choose_best_model(
  data_name = "Canaries",
  model = "cr_dd"
)
canaries_rr_lac_dd <- choose_best_model(
  data_name = "Canaries",
  model = "rr_lac_dd"
)
canaries_rr_mu_dd <- choose_best_model(
  data_name = "Canaries",
  model = "rr_mu_dd"
)
canaries_rr_k <- choose_best_model(
  data_name = "Canaries",
  model = "rr_k"
)
canaries_rr_gam_dd <- choose_best_model(
  data_name = "Canaries",
  model = "rr_gam_dd"
)
canaries_rr_laa_dd <- choose_best_model(
  data_name = "Canaries",
  model = "rr_laa_dd"
)

# Comoros

comoros_cr_dd <- choose_best_model(
  data_name = "Comoros",
  model = "cr_dd"
)
comoros_rr_lac_dd <- choose_best_model(
  data_name = "Comoros",
  model = "rr_lac_dd"
)
comoros_rr_mu_dd <- choose_best_model(
  data_name = "Comoros",
  model = "rr_mu_dd"
)
comoros_rr_k <- choose_best_model(
  data_name = "Comoros",
  model = "rr_k"
)
comoros_rr_gam_dd <- choose_best_model(
  data_name = "Comoros",
  model = "rr_gam_dd"
)
comoros_rr_laa_dd <- choose_best_model(
  data_name = "Comoros",
  model = "rr_laa_dd"
)

# Galapagos

galapagos_cr_dd <- choose_best_model(
  data_name = "Galapagos",
  model = "cr_dd"
)
galapagos_rr_lac_dd <- choose_best_model(
  data_name = "Galapagos",
  model = "rr_lac_dd"
)
galapagos_rr_mu_dd <- choose_best_model(
  data_name = "Galapagos",
  model = "rr_mu_dd"
)
galapagos_rr_k <- choose_best_model(
  data_name = "Galapagos",
  model = "rr_k"
)
galapagos_rr_gam_dd <- choose_best_model(
  data_name = "Galapagos",
  model = "rr_gam_dd"
)
galapagos_rr_laa_dd <- choose_best_model(
  data_name = "Galapagos",
  model = "rr_laa_dd"
)
galapagos_cr_dd_2type_lac <- choose_best_model(
  data_name = "Galapagos_2type",
  model = "cr_dd_2type_lac"
)
galapagos_cr_dd_2type_mu <- choose_best_model(
  data_name = "Galapagos_2type",
  model = "cr_dd_2type_mu"
)
galapagos_cr_dd_2type_k <- choose_best_model(
  data_name = "Galapagos_2type",
  model = "cr_dd_2type_k"
)
galapagos_cr_dd_2type_lac_mu <- choose_best_model(
  data_name = "Galapagos_2type",
  model = "cr_dd_2type_lac_mu"
)
galapagos_cr_dd_2type_lac_k <- choose_best_model(
  data_name = "Galapagos_2type",
  model = "cr_dd_2type_lac_k"
)
galapagos_cr_dd_2type_mu_k <- choose_best_model(
  data_name = "Galapagos_2type",
  model = "cr_dd_2type_mu_k"
)
galapagos_cr_dd_2type_lac_mu_k <- choose_best_model(
  data_name = "Galapagos_2type",
  model = "cr_dd_2type_lac_mu_k"
)

# Hawaii

hawaii_cr_dd <- choose_best_model(
  data_name = "Hawaii",
  model = "cr_dd"
)
hawaii_rr_lac_dd <- choose_best_model(
  data_name = "Hawaii",
  model = "rr_lac_dd"
)
hawaii_rr_mu_dd <- choose_best_model(
  data_name = "Hawaii",
  model = "rr_mu_dd"
)
hawaii_rr_k <- choose_best_model(
  data_name = "Hawaii",
  model = "rr_k"
)
hawaii_rr_gam_dd <- choose_best_model(
  data_name = "Hawaii",
  model = "rr_gam_dd"
)
hawaii_rr_laa_dd <- choose_best_model(
  data_name = "Hawaii",
  model = "rr_laa_dd"
)
hawaii_cr_dd_2type_lac <- choose_best_model(
  data_name = "Hawaii_2type",
  model = "cr_dd_2type_lac"
)
hawaii_cr_dd_2type_mu <- choose_best_model(
  data_name = "Hawaii_2type",
  model = "cr_dd_2type_mu"
)
hawaii_cr_dd_2type_k <- choose_best_model(
  data_name = "Hawaii_2type",
  model = "cr_dd_2type_k"
)
hawaii_cr_dd_2type_lac_mu <- choose_best_model(
  data_name = "Hawaii_2type",
  model = "cr_dd_2type_lac_mu"
)
hawaii_cr_dd_2type_lac_k <- choose_best_model(
  data_name = "Hawaii_2type",
  model = "cr_dd_2type_lac_k"
)
hawaii_cr_dd_2type_mu_k <- choose_best_model(
  data_name = "Hawaii_2type",
  model = "cr_dd_2type_mu_k"
)
hawaii_cr_dd_2type_lac_mu_k <- choose_best_model(
  data_name = "Hawaii_2type",
  model = "cr_dd_2type_lac_mu_k"
)

# Marquesas

marquesas_cr_dd <- choose_best_model(
  data_name = "Marquesas",
  model = "cr_dd"
)
marquesas_rr_lac_dd <- choose_best_model(
  data_name = "Marquesas",
  model = "rr_lac_dd"
)
marquesas_rr_mu_dd <- choose_best_model(
  data_name = "Marquesas",
  model = "rr_mu_dd"
)
marquesas_rr_k <- choose_best_model(
  data_name = "Marquesas",
  model = "rr_k"
)
marquesas_rr_gam_dd <- choose_best_model(
  data_name = "Marquesas",
  model = "rr_gam_dd"
)
marquesas_rr_laa_dd <- choose_best_model(
  data_name = "Marquesas",
  model = "rr_laa_dd"
)
marquesas_cr_dd_2type_lac <- choose_best_model(
  data_name = "Marquesas_2type",
  model = "cr_dd_2type_lac"
)
marquesas_cr_dd_2type_mu <- choose_best_model(
  data_name = "Marquesas_2type",
  model = "cr_dd_2type_mu"
)
marquesas_cr_dd_2type_k <- choose_best_model(
  data_name = "Marquesas_2type",
  model = "cr_dd_2type_k"
)
marquesas_cr_dd_2type_lac_mu <- choose_best_model(
  data_name = "Marquesas_2type",
  model = "cr_dd_2type_lac_mu"
)
marquesas_cr_dd_2type_lac_k <- choose_best_model(
  data_name = "Marquesas_2type",
  model = "cr_dd_2type_lac_k"
)
marquesas_cr_dd_2type_mu_k <- choose_best_model(
  data_name = "Marquesas_2type",
  model = "cr_dd_2type_mu_k"
)
marquesas_cr_dd_2type_lac_mu_k <- choose_best_model(
  data_name = "Marquesas_2type",
  model = "cr_dd_2type_lac_mu_k"
)

# New Caledonia

new_caledonia_cr_dd <- choose_best_model(
  data_name = "New_Caledonia",
  model = "cr_dd"
)
new_caledonia_rr_lac_dd <- choose_best_model(
  data_name = "New_Caledonia",
  model = "rr_lac_dd"
)
new_caledonia_rr_mu_dd <- choose_best_model(
  data_name = "New_Caledonia",
  model = "rr_mu_dd"
)
new_caledonia_rr_k <- choose_best_model(
  data_name = "New_Caledonia",
  model = "rr_k"
)
new_caledonia_rr_gam_dd <- choose_best_model(
  data_name = "New_Caledonia",
  model = "rr_gam_dd"
)
new_caledonia_rr_laa_dd <- choose_best_model(
  data_name = "New_Caledonia",
  model = "rr_laa_dd"
)

# Sao Tome Principe

saotome_principe_cr_dd <- choose_best_model(
  data_name = "SaoTome_Principe",
  model = "cr_dd"
)
saotome_principe_rr_lac_dd <- choose_best_model(
  data_name = "SaoTome_Principe",
  model = "rr_lac_dd"
)
saotome_principe_rr_mu_dd <- choose_best_model(
  data_name = "SaoTome_Principe",
  model = "rr_mu_dd"
)
saotome_principe_rr_k <- choose_best_model(
  data_name = "SaoTome_Principe",
  model = "rr_k"
)
saotome_principe_rr_gam_dd <- choose_best_model(
  data_name = "SaoTome_Principe",
  model = "rr_gam_dd"
)
saotome_principe_rr_laa_dd <- choose_best_model(
  data_name = "SaoTome_Principe",
  model = "rr_laa_dd"
)