median of otus per category

Since the relative abundance of otus does not follow a normal distribution, the median is a better description of the average value in each category (treatment, e.g.) than the mean.
However, in mothur, I have not seen median, only mean values. For example in metastats.
Is there a function that could output the median values of each otu per category?
In advance, thanks for your help
Bernard

You’re probably stuck with learning R. Which is likely a good thing :slight_smile: