SOPIE-package | R Documentation |
The package 'SOPIE' provides functions to non-parametrically estimate the off-pulse interval of a source function originating from a pulsar.
This technique is based on a sequential application of P-values
obtained from goodness-of-fit tests for the uniform distribution. The well-known Kolmogorov-Smirnov, Cramer-von Mises,
Anderson-Darling and Rayleigh test statistics are applied sequentially on subintervals of [0 ; 1].
The most important functions in the package
are combined in a wrapper function called SOPIE
. Users should start by looking at the documentation of the function findh
, circ.kernel
and SOPIE
.
Package: | SOPIE |
Type: | Package |
Version: | 1.6 |
Date: | 2022-02-23 |
License: | GPL-3 |
LazyLoad: | yes |
The SOPIE package consists of 4 main functions. Each of these functions are
discussed in terms of its functioning, structure, arguments and output in the help documentation of each function.
findh
is the function used to obtain the estimated smoothing parameter \hat h that will be
used in the circular kernel density estimator.
circ.kernel
is 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.estimate
and b.estimate
is almost identical function. a.estimate
is 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.estimate
is the function
used to obtain the estimated values of b, i.e. \hat b, for the off-pulse interval of a pulsar light curve.
SOPIE
is 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.
Willem Daniel Schutte
Maintainer: Willem Daniel Schutte
Jammalamadaka, S. Rao and SenGupta, A. (2001). Topics in Circular Statistics, World Scientific Publishing Co. Pte. Ltd.
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
Schutte WD, Swanepoel JWH (2016). SOPIE: an R package for the non-parametric estimation of the off-pulse interval of a pulsar light curve. Monthly Notices of the Royal Astronomical Society, 461, 627-640.
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)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.