written my understanding of Neighbors2; TODO: sameSeasonSimilarIndices

[talweg.git] / pkg / R / utils.R

diff --git a/pkg/R/utils.R b/pkg/R/utils.R
index 712a4f8..f79c300 100644 (file)
--- a/pkg/R/utils.R
+++ b/pkg/R/utils.R
@@ -58,9 +58,10 @@ integerIndexToDate = function(index, data)
 #' @param index Day index (numeric or date)
 #' @param limit Maximum number of indices to return
 #' @param same_season Should the indices correspond to day in same season?
+#' @param data Dataset is required for a search in same season
 #'
 #' @export
-getSimilarDaysIndices = function(index, limit, same_season)
+getSimilarDaysIndices = function(index, limit, same_season, data=NULL)
 {
        index = dateIndexToInteger(index)
 
@@ -72,7 +73,9 @@ getSimilarDaysIndices = function(index, limit, same_season)
                return ( rep(index,nb_days) - 7*seq_len(nb_days) )
        }
 
-       #Look for similar days in similar season (+/- 30 days)
+
+       #TODO: use data... 12-12-1-2 CH, 3-4-9-10 EP et le reste NP
+       #Look for similar days in similar season
        days = c()
        i = index
        while (i >= 1 && length(days) < limit)
@@ -93,6 +96,8 @@ getSimilarDaysIndices = function(index, limit, same_season)
                # TODO: exact computation instead of -364
                # 364 = closest multiple of 7 to 365 - drift along the years... but not so many years so OK
                i = i - 364
+
+
        }
 
        return ( days[1:min(limit,length(days))] )