R/clean.timeseries.R

In brisneve/THE: Thermal History Exploration

clean.timeseries <- function(x,y, start_date, end_date, by = "week"){

  if(!require(tidyverse)){install.packages("tidyverse"); library(tidyverse)}
  if(!require(imputeTS)){install.packages("imputeTS"); library(imputeTS)}
  if(!require(lubridate)){install.packages("lubridate"); library(lubridate)}

  ts.df <-as.data.frame(list(date = x, estimate = y))
  nd <- as.data.frame(seq.Date(start_date, end_date, by))
  colnames(nd) <- "date"
  ts.df1 <- full_join(ts.df,nd,by="date")
  ts.df1$date <- round_date(ts.df1$date, by)
  ts.df1 <- ts.df1[order(ts.df1$date),]

  ts.df2 <- ts.df1 %>%
    group_by(date) %>%
    summarise(estimate = mean(estimate))
  #ts.df2$estimate <-   ifelse(unique(is.na(ts.df2$estimate)) == FALSE, ts.df2$estimate, na.interpolation(ts.df2$estimate))
  ts.df2
}