- distances = sapply(indices, function(i) abs(data$getLevel(i)-levelToday))
- # 2 and 5 below == magic numbers (determined by Bruno & Michel)
- same_pollution = (distances <= 2)
- if (sum(same_pollution) == 0)
+ 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)