Function to compute the kernel density estimation on a grid of the simplex
kdeDirichlet(x, adj = 1, n = 200, kdegrid = NULL, delta = FALSE)A list of two or three elements, depending on the value of `kdegrid`. If `kdegrid` is null, the function is assumed to be used for two-dimensional plotting, and the output is one compatible with the function [image()], i.e. a list of three elements (vector of x-values, vector of y-values and matrix of density values computed). If `kdegrid` is a grid, then the output has two elements: the input grid and a vector of computed densities.
NOTE: no effort is made to check that `kdegrid` has the right class, dimension or content.
data set of (complete) compositional data, i.e. data summing to 1 by columns
accessory scaling factor, for modifying the bandwith in analogy to function [MASS::kde()]
integer, number of grid nodes on each component, where to estimate the density; ignored if `kdegrid` is given
data frame, set of locations where to estimate the density; either specify `n` or `kdegrid`
logical or real controlling if/how zeroes in `x` are treated; logical works only for `kdegrid=NULL` and uses correction by half the grid cell size, otherwise give a real value; this value will be added to the whole composition, including the non-zero values.
This function computes the kde (kernel density estimation) of the probability density function of a random composition on the simplex, by using Dirichlet kernels. The method was proposed by Aitchison and Lauder (1985).
Aitchison J., Lauder I.J. (1985) Kernel density estimation for compositional data; _J. Roy. Statist. Soc. Ser. C_, 34 (2): 129-137.
Ouimet, F. and Tolosana-Delgado, R. (2022) Asymptotic properties of Dirichlet kernel density estimators; _Journal of Multivariate Analysis_ 187: 104832, tools:::Rd_expr_doi("10.1016/j.jmva.2021.104832")