\begin{Schunk}
\begin{Sinput}
 sumdat <- apply(dat, 2, quantile, na.rm = TRUE)
 sumdat <- rbind(sumdat, mean= apply(dat, 2, mean, na.rm = TRUE))
 sumdat <- rbind(sumdat, sd= apply(dat, 2, sd, na.rm = TRUE))
 sumdat <- rbind(sumdat, var= apply(dat, 2, var, na.rm = TRUE))
 sumdat
\end{Sinput}
\begin{Soutput}
          GRIP       ARM   RATINGS       SIMS
0%    29.00000  19.00000 21.600000 -4.1700000
25%   94.00000  64.50000 34.800000 -0.9650000
50%  111.00000  81.50000 41.300000  0.1600000
75%  124.50000  94.00000 47.700000  1.0700000
100% 189.00000 132.00000 57.200000  5.1700000
mean 110.23129  78.75170 41.009878  0.2017687
sd    23.62987  21.10933  8.521865  1.6789742
var  558.37079 445.60402 72.622184  2.8189544
\end{Soutput}
\end{Schunk}
