ptspec: The Periodogram of Point Process Data

In SAPP: Statistical Analysis of Point Processes

The Periodogram of Point Process Data

Description

Provide the periodogram of point process data with the significant band (0.90, 0.95 and 0.99) of the maximum power in searching a cyclic component, for stationary Poisson Process.

Usage

ptspec(data, nfre, prdmin, prd, nsmooth = 1, pprd, interval, plot = TRUE)

Arguments

data	data of events.
nfre	number of sampling frequencies of spectra.
prdmin	the minimum periodicity of the sampling.
prd	a periodicity for calculating the Rayleigh probability.
nsmooth	number for smoothing of periodogram.
pprd	particular periodicities to be investigated among others.
interval	length of observed time interval of events.
plot	logical. If TRUE (default) the periodogram is plotted.

Value

f	frequency.
db	D.B.
power	power.
rayleigh.prob	the probability of Rayleigh.
distance	= \sqrt(rwx^2+rwy^2).
phase	phase.

References

Ogata, Y., Katsura, K. and Zhuang, J. (2006) Computer Science Monographs, No.32, TIMSAC84: STATISTICAL ANALYSIS OF SERIES OF EVENTS (TIMSAC84-SASE) VERSION 2. The Institute of Statistical Mathematics.

Examples

data(Brastings)  # The Occurrence Times Data of 627 Blastings
ptspec(Brastings, nfre = 1000, prdmin = 0.5, prd = 1.0, pprd = c(2.0, 1.0, 0.5),
       interval = 4600)

data(PoissonData)  # to see the contrasting difference
ptspec(PoissonData, nfre = 1000, prdmin = 0.5, prd =  1.0, pprd = c(2.0, 1.0, 0.5),
       interval = 5000)

SAPP documentation built on June 7, 2023, 5:45 p.m.