Mam następującą ramkę danych (df) z 29 obserwacji 5 zmiennych:

    age   height_seca1 height_chad1 height_DL weight_alog1
1   19         1800         1797       180           70
2   19         1682         1670       167           69
3   21         1765         1765       178           80
4   21         1829         1833       181           74
5   21         1706         1705       170          103
6   18         1607         1606       160           76
7   19         1578         1576       156           50
8   19         1577         1575       156           61
9   21         1666         1665       166           52
10  17         1710         1716       172           65
11  28         1616         1619       161           66
12  22         1648         1644       165           58
13  19         1569         1570       155           55
14  19         1779         1777       177           55
15  18         1773         1772       179           70
16  18         1816         1809       181           81
17  19         1766         1765       178           77
18  19         1745         1741       174           76
19  18         1716         1714       170           71
20  21         1785         1783       179           64
21  19         1850         1854       185           71
22  31         1875         1880       188           95
23  26         1877         1877       186          106
24  19         1836         1837       185          100
25  18         1825         1823       182           85
26  19         1755         1754       174           79
27  26         1658         1658       165           69
28  20         1816         1818       183           84
29  18         1755         1755       175           67

Chcę uzyskać średnią, odchylenie standardowe, medianę, minimalną, maksymalną i wielkość próbki każdej ze zmiennych i uzyskać wynik jako ramkę danych. Próbowałem użyć poniższego kodu, ale wtedy praca z tapply lub agregacją staje się niemożliwa. Wydaje mi się, że nie jestem początkującym programistą R. Moje zlecenie nie wymaga ode mnie żadnych dodatkowych pakietów R.

apply(df, 2, mean)
apply(df, 2, sd)
apply(df, 2, median)
apply(df, 2, min)
apply(df, 2, max)
apply(df, 2, length)

W idealnej sytuacji tak powinna wyglądać ramka danych wyjściowych, w tym nagłówki wierszy dla każdej funkcji statystycznej:

             age height_seca1 height_chad1    height_DL weight_alog1 
mean          20         1737         1736          173           73
sd            3.3         91.9         92.7          9.7         14.5 
median        19         1755         1755          175           71
minimum       17         1569         1570          155           50 
maximum       31         1877         1880          188          106
sample size   29           29           29           29           29 

Każda pomoc byłaby bardzo mile widziana.