circ_ksdensity FUNCTION - Compute a kernel density estimate over a periodic domain
Usage: [vfEstimate] = circ_ksdensity(vfObservations, vfPDFSamples, <vfDomain, fSigma, vfWeights>)
This function calculates a kernel density estimate of an (optionally weighted) data sample, over a periodic domain.
'vfObservations' is a set of observations made over a periodic domain, optionally defined by 'vfDomain': [fMin fMax]. The default domain is [0..2*pi]. 'vfPDFSamples' defines the sample points over which to perform the kernel density estimate, over the same domain as 'vfObservations'.
Weighted estimations can be performed by providing the optional argument 'vfWeights', where each element in 'vfWeights' corresponds to the matching element in 'vfObservations'.
The kernel density estimate will be performed using a wrapped Gaussian kernel, with a width estimated as
(4/3)^0.2 * circ_std(vfObservations, vfWeights) *(length(vfObservations^-0.2)
The optional argument 'fSigma' can be provided to set the width of the kernel.
'vfEstimate' will be a vector with a (weighted) estimate of the underlying distribution, with an entry for each element of 'vfPDFSamples'.