+ if (local)
+ {
+ # Neighbors: days in "same season"; TODO: 60 == magic number...
+ fdays = getSimilarDaysIndices(today, data, limit=60, same_season=TRUE,
+ days_in=fdays_cut)
+ if (length(fdays) <= 1)
+ return (NA)
+ levelToday = data$getLevel(today)
+ distances = sapply(fdays, function(i) abs(data$getLevel(i)-levelToday))
+ #TODO: 2, 3, 5, 10 magic numbers here...
+ dist_thresh = 2
+ min_neighbs = min(3,length(fdays))
+ repeat
+ {
+ same_pollution = (distances <= dist_thresh)
+ nb_neighbs = sum(same_pollution)
+ if (nb_neighbs >= min_neighbs) #will eventually happen
+ break
+ dist_thresh = dist_thresh + 3
+ }
+ fdays = fdays[same_pollution]
+ max_neighbs = 10
+ if (nb_neighbs > max_neighbs)
+ {
+ # Keep only max_neighbs closest neighbors
+ fdays = fdays[
+ sort(distances[same_pollution],index.return=TRUE)$ix[1:max_neighbs] ]
+ }
+ if (length(fdays) == 1) #the other extreme...
+ {
+ if (final_call)
+ {
+ private$.params$weights <- 1
+ private$.params$indices <- fdays
+ private$.params$window <- 1
+ }
+ return ( data$getSerie(fdays[1])[1:horizon] ) #what else?!
+ }
+ }
+ else
+ fdays = fdays_cut #no conditioning
+