> data.frame(n = tapply(v.move, Nn, length), + med = tapply(v.move, Nn, median, na.rm = TRUE)) n med 0 86 150 1 70 200 2 40 100 3 13 220 4 5 265 5 7 285 6 5 495 7 3 20 8 1 NA 10 1 420 11 1 90