D_05_Figure_S12c.R

Libraries

Code

suppressPackageStartupMessages({
library(hstats)
library(here)
library(dplyr)
library(tidymodels)
library(ranger)
library(ggplot2)
library(ggthemes)
})

Interactions & their partial dependencies

Code

data <- 
  readRDS(here::here("Data/output/B_models/",
                     "B_01_list_all_results_rf.rds"))


fit_3 <- 
  data$res[[3]] %>% 
  extract_fit_parsnip()


train <- 
  data$list_split_data[[3]]$train

pdp_lnR_y_interactions <- 
  list()

All 20 interactions from the yearly rate:

Pairwise Interactions

1

Code

pdp_lnR_y_interactions[[1]] <- 
  partial_dep(
    object = fit_3,
    X = train,
    v = c("Mass", "AOO"),
    BY = "datasetID",
    grid_size = 100L
  ) %>%
  plot() +
  ggtitle("Mass ~ AOO ~ dataset")+
  ggthemes::theme_few()

2

Code

pdp_lnR_y_interactions[[2]] <- 
  partial_dep(
    object = fit_3,
    X = train,
    v = c("D_AOO_a", "AOO"),
    BY = "datasetID",
    grid_size = 100L
  ) %>%
  plot() +
  ggtitle("D ~ AOO ~ dataset")+
  ggthemes::theme_few()

3

Code

pdp_lnR_y_interactions[[3]] <- 
  partial_dep(
    object = fit_3,
    X = train,
    v = c("Mass", "joincount_delta"),
    BY = "datasetID",
    grid_size = 100L
  ) %>%
  plot() +
  ggtitle("Mass ~ jc ~ dataset")+
  ggthemes::theme_few()

4

Code

pdp_lnR_y_interactions[[4]] <- 
  partial_dep(
    object = fit_3,
    X = train,
    v = c("D_AOO_a", "mean_lnLac"),
    BY = "datasetID",
    grid_size = 100L
  ) %>%
  plot() +
  ggtitle("D ~ lac ~ dataset")+
  ggthemes::theme_few()

5

Code

pdp_lnR_y_interactions[[5]] <- 
  partial_dep(
    object = fit_3,
    X = train,
    v = c("AOO", "circNorm"),
    BY = "datasetID",
    grid_size = 100L
  ) %>%
  plot() +
  ggtitle("AOO ~ circNorm ~ dataset")+
  ggthemes::theme_few()

6

Code

pdp_lnR_y_interactions[[6]] <- 
  partial_dep(
    object = fit_3,
    X = train,
    v = c("Mass", "D_AOO_a"),
    BY = "datasetID",
    grid_size = 100L
  ) %>%
  plot() +
  ggtitle("D ~ D ~ dataset")+
  ggthemes::theme_few()

7

Code

pdp_lnR_y_interactions[[7]] <- 
  partial_dep(
    object = fit_3,
    X = train,
    v = c("joincount_delta", "circNorm"),
    BY = "datasetID",
    grid_size = 100L
  ) %>%
  plot() +
  ggtitle("jc ~ circ ~ dataset")+
  ggthemes::theme_few()

8

Code

pdp_lnR_y_interactions[[8]] <- 
  partial_dep(
    object = fit_3,
    X = train,
    v = c("D_AOO_a", "circNorm"),
    BY = "datasetID",
    grid_size = 100L
  ) %>%
  plot() +
  ggtitle("D ~ circ ~ dataset")+
  ggthemes::theme_few()

9

Code

pdp_lnR_y_interactions[[9]] <- 
  partial_dep(
    object = fit_3,
    X = train,
    v = c("Mass", "circNorm"),
    BY = "datasetID",
    grid_size = 100L
  ) %>%
  plot() +
  ggtitle("Mass ~ circ ~ dataset")+
  ggthemes::theme_few()

10

Code

pdp_lnR_y_interactions[[10]] <- 
  partial_dep(
    object = fit_3,
    X = train,
    v = c("AOO", "joincount_delta"),
    BY = "datasetID",
    grid_size = 100L
  ) %>%
  plot() +
  ggtitle("AOO ~ jc ~ dataset")+
  ggthemes::theme_few()

Threeway Interactions

11

Code

pdp_lnR_y_interactions[[11]] <- 
  partial_dep(
    object = fit_3,
    X = train,
    v = c("Mass", "D_AOO_a"),
    BY = "AOO",
    grid_size = 100L,
    strategy = "uniform"
  ) %>%
  plot() +
  ggtitle("Mass ~ D ~ AOO")+
  ggthemes::theme_few()

12

Code

pdp_lnR_y_interactions[[12]] <- 
  partial_dep(
    object = fit_3,
    X = train,
    v = c("D_AOO_a", "AOO"),
    BY = "circNorm",
    grid_size = 100L,
    strategy = "uniform"
  ) %>%
  plot() +
  ggtitle("D ~ AOO ~ circ")+
  ggthemes::theme_few()

13

Code

pdp_lnR_y_interactions[[13]] <- 
  partial_dep(
    object = fit_3,
    X = train,
    v = c("Mass", "D_AOO_a"),
    BY = "joincount_delta",
    grid_size = 100L,
    strategy = "uniform"
  ) %>%
  plot() +
  ggtitle("Mass ~ D ~ jc")+
  ggthemes::theme_few()

14

Code

pdp_lnR_y_interactions[[14]] <- 
  partial_dep(
    object = fit_3,
    X = train,
    v = c("Mass", "AOO"),
    BY = "circNorm",
    grid_size = 100L,
    strategy = "uniform"
  ) %>%
  plot() +
  ggtitle("Mass ~ AOO ~ circNorm")+
  ggthemes::theme_few()

15

Code

pdp_lnR_y_interactions[[15]] <- 
  partial_dep(
    object = fit_3,
    X = train,
    v = c("Mass", "AOO"),
    BY = "joincount_delta",
    grid_size = 100L,
    strategy = "uniform"
  ) %>%
  plot() +
  ggtitle("Mass ~ AOO ~ jc")+
  ggthemes::theme_few()

16

Code

pdp_lnR_y_interactions[[16]] <- 
  partial_dep(
    object = fit_3,
    X = train,
    v = c("Mass", "joincount_delta"),
    BY = "circNorm",
    grid_size = 100L,
    strategy = "uniform"
  ) %>%
  plot() +
  ggtitle("Mass ~ jc ~ circ")+
  ggthemes::theme_few()

17

Code

pdp_lnR_y_interactions[[17]] <- 
  partial_dep(
    object = fit_3,
    X = train,
    v = c("D_AOO_a", "AOO"),
    BY = "joincount_delta",
    grid_size = 100L,
    strategy = "uniform"
  ) %>%
  plot() +
  ggtitle("D ~ AOO ~ jc")+
  ggthemes::theme_few()

18

Code

pdp_lnR_y_interactions[[18]] <- 
  partial_dep(
    object = fit_3,
    X = train,
    v = c("Mass", "D_AOO_a"),
    BY = "circNorm",
    grid_size = 100L,
    strategy = "uniform"
  ) %>%
  plot() +
  ggtitle("Mass ~ D ~ circ")+
  ggthemes::theme_few()

19

Code

pdp_lnR_y_interactions[[19]] <- 
  partial_dep(
    object = fit_3,
    X = train,
    v = c("D_AOO_a", "joincount_delta"),
    BY = "circNorm",
    grid_size = 100L,
    strategy = "uniform"
  ) %>%
  plot() +
  ggtitle("D ~ jc ~ circ")+
  ggthemes::theme_few()

20

Code

pdp_lnR_y_interactions[[20]] <- 
  partial_dep(
    object = fit_3,
    X = train,
    v = c("AOO", "joincount_delta"),
    BY = "circNorm",
    grid_size = 100L,
    strategy = "uniform"
  ) %>%
  plot() +
  ggtitle("AOO ~ jc ~ circ")+
  ggthemes::theme_few()

Code

saveRDS(pdp_lnR_y_interactions, here("Data/output/temp/D_05_lnR_y_pdp_all.rds"))

Code

lapply(pdp_lnR_y_interactions, plot)

[[1]]

[[2]]

[[3]]

[[4]]

[[5]]

[[6]]

[[7]]

[[8]]

[[9]]

[[10]]

[[11]]

[[12]]

[[13]]

[[14]]

[[15]]

[[16]]

[[17]]

[[18]]

[[19]]

[[20]]