#
# Licensed to the Apache Software Foundation (ASF) under one or more
# contributor license agreements.  See the NOTICE file distributed with
# this work for additional information regarding copyright ownership.
# The ASF licenses this file to You under the Apache License, Version 2.0
# (the "License"); you may not use this file except in compliance with
# the License.  You may obtain a copy of the License at
#
#    http://www.apache.org/licenses/LICENSE-2.0
#
# Unless required by applicable law or agreed to in writing, software
# distributed under the License is distributed on an "AS IS" BASIS,
# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
# See the License for the specific language governing permissions and
# limitations under the License.


ams_tukeys <- function(train_data, test_data, n) {

#  res <- get_data(url)
#  data <- data.frame(as.numeric(names(res$metrics[[1]]$metrics)), as.numeric(res$metrics[[1]]$metrics))
#  names(data) <- c("TS", res$metrics[[1]]$metricname)
#  train_data <- data[which(data$TS >= train_start & data$TS <= train_end), 2]
#  test_data <- data[which(data$TS >= test_start & data$TS <= test_end), ]

  anomalies <- data.frame()
  quantiles <- quantile(train_data[,2])
  iqr <- quantiles[4] - quantiles[2]
  niqr <- 0

  for ( i in 1:length(test_data[,1])) {
    x <- test_data[i,2]
    lb <- quantiles[2] - n*iqr
    ub <- quantiles[4] + n*iqr
    if ( (x < lb)  || (x > ub) ) {
      if (iqr != 0) {
        if (x < lb) {
          niqr <- (quantiles[2] - x) / iqr
        } else {
          niqr <- (x - quantiles[4]) / iqr
        }
      }
        anomaly <- c(test_data[i,1], x, niqr)
        anomalies <- rbind(anomalies, anomaly)
      }
  }
  if(length(anomalies) > 0) {
    names(anomalies) <- c("TS", "Value", "niqr")
  }
  return (anomalies)
}