s_data_plot.R

library(tidyverse)
library(lubridate)
s_col_names <- c("sample_id_idx", "Visit_no", "patch_id", "visit_duration", "exit_time", "entry_time", "Between_visit_duration",
                 "ID", "Species", "Patch_centroid_x", "Patch_centroid_y", "revisit_ID", "Time_since_last_visit", "trip_vedba",
                 "trip_vedba_mean", "total_patch_vedba", "mean_patch_vedba", "vedba_above_35", "vedba_below_35", "VID",
                 "Weight", "Sex", "Year", "week", "Study_day", "Week_year")

s_data <- read_csv("data/Revisits_BTH_edits_included.csv",
                   col_names = s_col_names,skip = 1) %>% 
  filter(Species == "Cebus capucinus") %>% 
  mutate(entry_time = with_tz(entry_time,tzone = "America/Panama"),
         exit_time = with_tz(exit_time,tzone = "America/Panama"),
         visit_mid = map2(.x = entry_time, .y = exit_time, .f = function(x,y){mean(c(x, y))}) %>%
         unlist() %>% as_datetime() %>% lubridate::with_tz(tzone = "America/Panama"),
         visit_duration_k = as.duration(seconds(visit_duration)),
         visit_mid_date = as_date(visit_mid),
         season_start = season_named[as.character(Year)],
         season_week = (as_datetime(visit_mid_date) - season_start)/ dweeks(1),
         season_week_bin = floor(season_week / week_bin_width) * week_bin_width,
         season_bin = str_c(Year,"-", str_pad(season_week_bin, width = 2, pad = 0)))

s_data %>%
  dplyr::select(ID, season_week_bin, visit_duration, visit_duration_k) %>% 
  ggplot(aes(x = season_week_bin,
             y = log10(visit_duration_k / as.duration(hours(1))))) +
  geom_density_2d_filled() +
  geom_hline(yintercept = log10(c(4,8)/60), linetype = 3, color = rgb(1,0,0,.8)) +
  labs(title = "Shauin") +
  facet_wrap(ID ~ ., nrow = 1) +
  coord_cartesian(xlim = c(0, 15), ylim = c(-3,1.2))
  
s_data %>%
  dplyr::select(ID, season_week_bin, visit_duration, visit_duration_k) %>% 
  filter(visit_duration_k > (8 * 60)) %>% 
  ggplot(aes(x = season_week_bin,
             y = log10(visit_duration_k / as.duration(hours(1))))) +
  geom_density_2d_filled() +
  geom_hline(yintercept = log10(c(4,8)/60), linetype = 3, color = rgb(1,0,0,.8)) +
  labs(title = "Shauin (filtered)") +
  facet_wrap(ID ~ ., nrow = 1) +
  coord_cartesian(xlim = c(0, 15), ylim = c(-3,1.2))