Dear list, Is it possible to replace the median in the following code so as to calculate 5 and 95% tiles? with( raunheim, aggregate(no, by = list(wd,hr), FUN = "median", na.rm=TRUE) ) Regards, Philip Sinclair