基于表格格式的R中另一个变量的许多变量的汇总统计数据？

Question

我有一个由10个变量组成的数据框。一个变量是站ID。我有10个站（A，B，C和J），每个站每天观察五年。我总共有18260个观测值。我的观察结果包含了很多天的NA。我的数据框看起来像这样

stationID      x1      x2     x3      x4       x10
A
A
A
A
B
B
B
B
C
C
C
C

J
J
J

我想根据电台名称获取每个变量的摘要统计信息（主要是n，mean，median和sd）。

我虽然关于

library(psych) describeBy(mydf, mydf$stationID)

我得到了所有变量的摘要统计信息，作为每个站的表。 但我希望得到类似的东西

                   station A        station B        station C    station J
var.2         n
           mean
             sd
         median
Var.3         n
           mean
             sd
         median

Var.10        n
           mean
             sd
         median

当变量有NA时，我没有得到摘要统计信息。如何显示所有变量的摘要统计信息输出，包括那些有NA？

的变量

我认为这个问题与其他建议的问题/答案不同，因为它有多个统计参数，我的观察结果包含NA