I\'ve attached two example datasets
resind2 <- structure(list(Station = structure(c(2L, 2L, 3L, 3L, 4L, 4L, 5L, 5L, 6L, 6L, 26L, 27L, 27L, 28L, 28L, 29L,