scripts/analysis.R

results <- read_results(file.path("results", "archipelagos41_paleo"))

ordered_results <- results |>
  dplyr::mutate(bic = calc_bic(loglik, df, 1000)) |>
  dplyr::group_by(age) |>
  dplyr::arrange(dplyr::desc(bic), .by_group = TRUE) |> dplyr::slice_min(bic)

View(ordered_results)
Neves-P/islandpaleoarea documentation built on Feb. 10, 2024, 8:46 a.m.