如何按因子计算多列的中位数？

我必须根据其因子计算数据集中各列的中位数。

这是我用来获取除“类型”列以外的各列中位数的代码。类型列是数据框中的第一列，正是该列的值是我的考虑因素。

quant0 = c(0.5)
Median = apply(mydata[2:1051], 2, median, probs = quant0, na.rm = TRUE )

我的数据框看起来像这样：

        Type    x1  x2  x3  ...
1:  Fresh   1.54    1.48    1.88    
2:  Dated   1.46    1.99    1.48
3:  Fresh   2.01    1.02    1.03
...

我想要因子新鲜和过时的x1，x2，...的中值。