PRT_MIT

library(dygraphs)
library(xts)

## Loading required package: zoo

## 
## Attaching package: 'zoo'

## The following objects are masked from 'package:base':
## 
##     as.Date, as.Date.numeric

main_dir = "processed_data/"
csv = "PRT_MIT"

wp_sheet = read.csv(paste0(main_dir,csv,"_sheet7.csv"))
met_sheet = read.csv(paste0(main_dir,csv,"_sheet10.csv"))

met_sheet$date = as.Date(as.character(met_sheet$date),format="%Y%m%d")
met_sheet$time = ifelse(is.na(met_sheet$time), "00:00:00", met_sheet$time)

met_sheet$datetime <- as.POSIXct(paste(met_sheet$date, met_sheet$time), format = "%Y-%m-%d %H:%M:%S")
wp_sheet$datetime <- as.POSIXct(paste(wp_sheet$date, wp_sheet$time), format = "%Y%m%d %H:%M:%S")

temp_met = met_sheet
temp_met$date = temp_met$datetime
temp_met$datetime = NULL
temp_met$time = NULL
temp_met$dataset_name = NULL
temp_wp = wp_sheet[,c("water_potential_mean","datetime")]
temp_wp$date = temp_wp$datetime
temp_wp$datetime = NULL

merged_data <- merge(temp_met, temp_wp, by = "date", all.x = TRUE)

dygraph(merged_data, main = "Interactive Plot", xlab = "Date", ylab = "Value") %>%
  dyRangeSelector()

data1 = merged_data[,c("date", "photosynthetically_active_radiation_ppfd_mumol_photon_m_2_s_1", "incident_shortwave_radiation", "net_radiation_m_s_1")]
data2 = merged_data[,c("date", "relative_humidity_percent", "windspeed_m_s")]
data3 = merged_data[,c("date", "precipitation_mm", "vapor_pressure_deficit_k_pa")]
data4 = merged_data[,c("date", "water_potential_mean")]

dygraph(data1, main = "All radiations", xlab = "Date", ylab = "Value") %>%
  dyRangeSelector()

dygraph(data2, main = "Humidity and windspeed", xlab = "Date", ylab = "Value") %>%
  dyRangeSelector()

dygraph(data3, main = "PPT and vapor pressure", xlab = "Date", ylab = "Value") %>%
  dyRangeSelector()

dygraph(data4, main = "Water potential", xlab = "Date", ylab = "Value") %>%
  dyRangeSelector()