scripts/dataset_comparison.R

In Neves-P/islandpaleoarea: Paleo-Area Influence on Island Evolutionary Models

data("dataset_comparison")

nature_old_vs_new_data <- (dataset_comparison[1:28, -1] - dataset_comparison[29:56, -1]) / (dataset_comparison[1:28, -1] + dataset_comparison[29:56, -1])
nature_comparison <- (dataset_comparison[1:28, -1] - dataset_comparison[57:84, -1]) / (dataset_comparison[1:28, -1] + dataset_comparison[57:84, -1])
new_vs_nature_new <- dataset_comparison[29:56, -1] - dataset_comparison[57:84, -1]

nature_old_vs_new_data <- dataset_comparison[1:19, -1] - dataset_comparison[29:47, -1]
nature_comparison <- dataset_comparison[1:28, -1] - dataset_comparison[57:84, -1]
new_vs_nature_new <- dataset_comparison[29:56, -1] - dataset_comparison[57:84, -1]

nature_old_vs_new_data_no_d0 <- nature_old_vs_new_data[-12]
nature_comparison_no_d0 <- nature_comparison[-12]
new_vs_nature_new_no_d0 <- new_vs_nature_new[-12]

View(nature_old_vs_new_data)
View(nature_comparison)
View(new_vs_nature_new)

nature_old_vs_new_data <- nature_old_vs_new_data |>
  tibble::rowid_to_column(var = "Model") |>
  tidyr::gather(key = "Parameter", value = "Difference", -1)
ggplot2::ggplot(nature_old_vs_new_data, ggplot2::aes(x = Model, y = Parameter, fill = Difference)) +
  ggplot2::geom_tile() +
  ggplot2::ggtitle("Nature data vs new data")

nature_comparison <- nature_comparison |>
  tibble::rowid_to_column(var = "Model") |>
  tidyr::gather(key = "Parameter", value = "Difference", -1)
ggplot2::ggplot(nature_comparison, ggplot2::aes(x = Model, y = Parameter, fill = Difference)) +
  ggplot2::geom_tile() +
  ggplot2::ggtitle("Nature data with new\nand old DAISIE")

new_vs_nature_new <- new_vs_nature_new |>
  tibble::rowid_to_column(var = "Model") |>
  tidyr::gather(key = "Parameter", value = "Difference", -1)
ggplot2::ggplot(new_vs_nature_new, ggplot2::aes(x = Model, y = Parameter, fill = Difference)) +
  ggplot2::geom_tile() +
  ggplot2::ggtitle("New data vs new nature")





nature_old_vs_new_data_no_d0 <- nature_old_vs_new_data_no_d0 |>
  tibble::rowid_to_column(var = "Model") |>
  tidyr::gather(key = "Parameter", value = "Difference", -1)
ggplot2::ggplot(nature_old_vs_new_data_no_d0, ggplot2::aes(x = Model, y = Parameter, fill = Difference)) +
  ggplot2::geom_tile() +
  ggplot2::ggtitle("Nature data vs new data")

nature_comparison_no_d0 <- nature_comparison_no_d0 |>
  tibble::rowid_to_column(var = "Model") |>
  tidyr::gather(key = "Parameter", value = "Difference", -1)
ggplot2::ggplot(nature_comparison_no_d0, ggplot2::aes(x = Model, y = Parameter, fill = Difference)) +
  ggplot2::geom_tile() +
  ggplot2::ggtitle("Nature data with new\nand old DAISIE")

new_vs_nature_new_no_d0 <- new_vs_nature_new_no_d0 |>
  tibble::rowid_to_column(var = "Model") |>
  tidyr::gather(key = "Parameter", value = "Difference", -1)
ggplot2::ggplot(new_vs_nature_new, ggplot2::aes(x = Model, y = Parameter, fill = Difference)) +
  ggplot2::geom_tile() +
  ggplot2::ggtitle("New data vs new nature")