SOPIE (version 1.4)

SOPIE: Sequential Off-Pulse Interval Estimation of a Pulsar Light Curve

Description

SOPIE is a wrapper-function that utilises findh, circ.kernel, a.estimate and b.estimate to produce the estimated off-pulse intervals in an easy readable matrix format, together with a graph.

Usage

SOPIE(data, h = 1, to = 1, alpha = 0.05, g = 20, r = 10, m = 1, grid = 512)

Arguments

data

the data vector within which to find the estimated smoothing parameter $\hat h$ that will be used in the circular kernel density estimator. After obtaining the minimum point(s) from the circular kernel density estimator, the estimate off-pulse interva

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$.

alpha

significance level ($\alpha$) that will be used during the sequential application of the goodness-of-fit tests for uniformity when estimating the off-pulse interval.

g

the value of the incremental growth of each subsequent interval over which uniformity is tested. In the suggested procedure, uniformity is sequentially tested, with the interval used in the test growing by g observations after every iterat

r

the number of subsequent intervals that must result in the rejection of uniformity before the function will stop. The choice of r must therefore be linked to the choice of g as explained above. For smaller values of g

m

the number of local minimum points included in the output.

grid

the number of equally spaced grid points at which the density is to be estimated.

Value

• The output produced by the function is a list containing the following:
• summary
  is a matrix that contains the estimated value of $a$ and $b$ for each of the four goodness-of-fit tests, namely the Anderson-Darling, Kolmogorov-Smirnov, Cramer-von Mises and the Rayleigh goodness-of-fit test. Based on the four estimated values of $a$ and $b$, the median values of $a$ and $b$ are also calculated. This median off-pulse interval is the recommended interval and also the interval that is depicted on the graph.
• general
  is a list containing the function call, the minimum value(s) used in the estimation, the level of significance ($\alpha$), the value of g and the value of r.
• A histogram estimate of the data is produced with the circular kernel density estimate overlaid. An indication of the estimated median off-pulse interval derived from the four goodness-of-fit tests is illustrated with two solid vertical lines.

Details

SOPIE is a wrapper-function in the sense that it utilises the function findh, circ.kernel, a.estimate and b.estimate to produce the estimated off-pulse intervals of a pulsar light curve in an easy readable matrix format, together with a graph consisting of the histogram estimate of the sample data, the kernel density estimator, and a visual representation of the estimated off-pulse intervals.

References

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

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)
data(crab)
SOPIE(crab)
data(J1709)
SOPIE(J1709)