SAPP (version 1.0.9-1)

ptspec: 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)

Value

f

frequency.

db

D.B.

power

power.

rayleigh.prob

the probability of Rayleigh.

distance

= \(\sqrt(rwx^2+rwy^2)\).

phase

phase.

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.

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)