PU
это точно, кстати для моды есть что-нибудь лаконичнее чем
?
d <- density(df$val)
moda_val <- d$x[which.max(d$y)]
?
это плохое решение, так как в случае мультимодальности which.max может отдать случайное или первое значение