Используя только функции base
R, вы можете использовать функцию by
для выполнения функций по группам, а затем cbind
список результатов в data.frame следующим образом:
do.call(cbind,
by(proj$Gross.Salary, proj$Year.of.Inflow,
function(x) quantile(x, c(0.25,0.5,0.75)))
)
пример данных:
set.seed(0L)
yrs <- rep(2009:2017, each=10)
proj <- data.frame(Year.of.Inflow=yrs, Gross.Salary=runif(length(yrs))*1e6)
выход: * +1010 *
2009 2010 2011 2012 2013 2014 2015 2016
25% 422306.3 249489.7 225912.1 375781.8 490352.4 262665.9 303301.2 357510.0
50% 644955.9 440901.5 384251.0 546553.6 600048.3 457858.4 432948.0 626836.0
75% 897966.6 709969.6 746002.4 709899.9 770278.0 626162.4 607789.2 858114.7
2017
25% 265454.4
50% 417326.9
75% 712166.3