R’s base graphics system is notable for the minimal design of its plots. Basic usage is very simple, although more complex customization capabilities are not user friendly. Hence I wrapped the plot and hist functions to improve their default behavior.
Any argument usually passed to plot or hist can also be passed to the two wrapper functions pretty_plot and pretty_hist. A comparison is shown below; “prettified” functions are on the right (obviously!).
Read More »
UQAM’s logo drawn in our bivariate angular density model. It is the type of picture we’d expect from convolution of an image with a gaussian kernel, but here it is obtained using finite mixtures. The difficulty lies in choosing the density basis in such a way as to constrain “bumpyness”, even though a finite number of “bumps” are summed together.
I use these grids to construct smooth approximations to bivariate angular probability densities (of interest in protein bioinformatics) and to specify semiparametric priors on bivariate angular density spaces. The grid on the Torus illustrates the proper angular wrapping behavior of planar grid. I drew them in Mathematica.
Code and variants below.
Read More »