library(dplyr)
data("Example_data", package = "Covid19Wastewater")
smoothing_df <- Example_data%>%
select(site, date, N1, N2)%>%
filter(N1 != 0, N2 != 0)%>%
mutate(N1 = log(N1), N2 = log(N2), N12_avg = (N1 + N2) / 2)
df_data <- loessSmoothMod(smoothing_df, "N12_avg", "N12_avg_loess", Filter = NULL)
head(df_data%>%
group_by(site)%>%
Flag_From_Trend( N12_avg, N12_avg_loess)%>%
select(site, date, N12_avg, flagged_outlier))
Run the code above in your browser using DataLab