Of course, you can always simply ask R to output summary statistics:

12345678910111213

> n <- 10000
> d <- data.frame(unif=runif(n=n), norm=rnorm(n=n), exp=rexp(n=n),
+ pois=rpois(n, lambda=1))
>
> summary(d)
> summary(d)
unif norm exp pois
Min. :2.449e-05 Min. :-3.603527 Min. :0.0001143 Min. :0.000
1st Qu.:2.489e-01 1st Qu.:-0.657498 1st Qu.:0.2963731 1st Qu.:0.000
Median :4.991e-01 Median :-0.003342 Median :0.6927739 Median :1.000
Mean :5.001e-01 Mean : 0.008527 Mean :0.9958194 Mean :1.006
3rd Qu.:7.498e-01 3rd Qu.: 0.662084 3rd Qu.:1.3718103 3rd Qu.:2.000
Max. :1.000e+00 Max. : 3.512077 Max. :9.4726291 Max. :6.000

Or perhaps, we could be even more detailed and ask for more finely detailed order statistics. The quantile function will calculate arbitrary order statistics:

Quantile quantile plots are typically used when you want to see if a sample is a particular distribution. You plot the quantiles of the sample versus the assumed distribution and compare; if they are the same distribution you should get a straight line at roughly a forty five degree angle.

12345678

> par(mfrow=c(1,2), oma=c(0,0,2,0))
>
> qqplot(x=d$exp, y=d$pois,
+ xlab='exponential', ylab='poisson', main = 'qq plot: different distributions')
> qqplot(x=d$norm, y=rnorm(n=n),
+ xlab='our normal sample in d', ylab='new normal sample', main = 'qq plot: same distribution, resampled')
>
> mtext('QQ plots -- different vs same distributions', side=3, outer=T)

Finally, for univariate distributions, I prefer box-percentile plots. These are similar to boxplots, but the width of the distribution graphs are proportional to the percent of observations more extreme in that direction. They are also marked at the 25th, 50th, and 75th percentiles.

And if the distributions are genuinely different, let’s examine 5 univariate Normal Distributions: our already sampled N(0,1), a resampled N(0,1), N(1,1), N(0,3), N(-1,0.5). Side by side, the box percentile plot really draws out the differences in the distributions.