Learn R Programming

SOPIE (version 1.4)

findh: Calculate the Estimated Smoothing Parameter

Description

This function is used to obtain the estimated smoothing parameter $\hat h$ that will be used in the circular kernel density estimator (see circ.kernel).

Usage

findh(data, h = 1, to = 1)

Arguments

data
the data vector from which to calculate the estimated smoothing parameter $\hat h$ that will be used in the circular kernel density estimator.
h
integer value from 1 to 9, specifying the smoothing parameter to calculate according to the following table: l{ $\hat h_1=1.06sn^{-1/5}$ $\hat h_2=1.06s_\circ n^{-1/5}$ $\hat h_3=1.06\bar{D}_\circ n^{-1/5}$ $\hat h_4=1.06|{D}_\circ|n
to
the value of the maximum domain of the data. Values will usually either be 1 or 2$\pi$.

Value

  • The function produces a single real value between 0 and 1, representing the rounded value (to 2 decimal places) of the estimating smoothing parameter.

References

Jammalamadaka S, SenGupta A (2001). "Topics in circular statistics." World Scientific Publishing Co. Pte. Ltd. Hall P, Watson G, Cabrera J (1987). "Kernel density estimation with spherical data." Biometrika, 74 (4), 751-762. Schutte WD (2014). Nonparametric estimation of the off-pulse interval(s) of a pulsar light curve. Ph.D. thesis, North-West University. URL http://hdl.handle.net/10394/12199

Examples

Run this code
simdata<-von_mises_sim(n=5000,k=1,c=0.3,noise=0.2)
findh(simdata,h=9,to=1)

Run the code above in your browser using DataLab