SAMOS Day 3

1 Start Day 3 lectures

library(tidyverse)

Read in the mangled data:

data_path <- "../data/SAMOS/other_data/SACTN_mangled.RData"
load(data_path)

ggplot(SACTN2, aes(x = date, y = temp)) +
  geom_line()

ggplot(data = SACTN1, aes(x = date, y = temp)) +
  geom_line(aes(colour = site, group = paste(site, "-", src))) +
  labs(x = "Date", y = "Temperature (°C)", colour = "Site")

paste(SACTN1$site, "-", SACTN1$src)[1:2]

[1] "Port Nolloth - DEA" "Port Nolloth - DEA"

SACTN2_tidy <- pivot_longer(SACTN2, cols = c("DEA", "KZNSB", "SAWS"),
                            names_to = "src",
                            values_to = "temp")
head(SACTN2_tidy)

# A tibble: 6 × 6
  site         date                depth type   src    temp
  <fct>        <dttm>              <dbl> <chr>  <chr> <dbl>
1 Port Nolloth 1998-01-01 00:00:00     0 thermo DEA    NA  
2 Port Nolloth 1998-01-01 00:00:00     0 thermo KZNSB  NA  
3 Port Nolloth 1998-01-01 00:00:00     0 thermo SAWS   12.7
4 Port Nolloth 1998-01-01 00:00:00     5 UTR    DEA    11.8
5 Port Nolloth 1998-01-01 00:00:00     5 UTR    KZNSB  NA  
6 Port Nolloth 1998-01-01 00:00:00     5 UTR    SAWS   NA  

SACTN3_tidy1 <- SACTN3 %>% 
  pivot_wider(names_from = "var", values_from = "val")

# unique(SACTN3$var)

SACTN4a_tidy <- SACTN4a |> 
  separate_wider_delim(index, names = c("site", "src"), delim = "/ ")

SACTN_tidy2 <- SACTN4a %>% 
  separate_wider_delim(index, names = c("site", "src"), delim = "/ ") %>% 
  mutate(year = lubridate::year(date))

SACTN4b_tidy <- SACTN4b |> 
  unite(col = "date", year, month, day, sep = "-") |> 
  select(site, src, date, depth, type)

SACTN4_tidy <- left_join(SACTN4a_tidy, SACTN4b_tidy)