X-Git-Url: https://git.auder.net/?p=talweg.git;a=blobdiff_plain;f=pkg%2FR%2Futils.R;h=f79c30083ec1c1ddcb4334976f6589efcf9fcc5e;hp=712a4f8d37f461be977f57a0c03f0b48d33d01e4;hb=9db234c56c330bb3f652718c5ee1eb16bc1f6fc7;hpb=2ae382665c51cdebbe3f29d58da0e69c494c4c91 diff --git a/pkg/R/utils.R b/pkg/R/utils.R index 712a4f8..f79c300 100644 --- 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))] )