The great findInterval() function in R uses left-closed sub-intervals in its vec argument, as shown in its docs:
findInterval()
vec
if i
i
If your limits are intervals you simply can grow the right interval a bit: interval+c(0,0.1) would do: findinterval(value, interval+c(0,0.1))