SOPIE. Users should start by looking at the documentation of the function findh, circ.kernel and SOPIE.findhis the function used to obtain the estimated smoothing parameter$\hat h$that will be
used in the circular kernel density estimator.circ.kernelis the function used to perform circular kernel density estimation on the
sample data set in order to obtain the minimum points of the kernel density estimator.
This is essentially the first step of the suggested procedure, as described in the second reference listed below. The output can also
be used to draw a graph of the circular kernel density estimator.a.estimateandb.estimateis almost identical function.a.estimateis the function
used to obtain the estimated values of$a$, i.e.$\hat a$, for the off-pulse interval of a pulsar light curve.b.estimateis the function
used to obtain the estimated values of$b$, i.e.$\hat b$, for the off-pulse interval of a pulsar light curve.SOPIEis a wrapper-function in the sense that it utilises all of the above function to produce
the estimated off-pulse intervals in an easy readable matrix format, together with
a graph consisting of the histogram estimate of the sample data, the kernel density
estimator and an indication of the estimated median off-pulse interval.set.seed(777)
simdata<-von_mises_sim(n=5000,k=1,c=0.3,noise=0.2)
SOPIE(simdata,h=1,to=1,alpha=0.05,g=5,r=10,m=1,grid=100)Run the code above in your browser using DataLab