SAPP (version 1.0.4)

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 )

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

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

Y.Ogata, K.Katsura and J.Zhuang (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 Brastings
  ptspec( Brastings, 1000, 0.5, 1.0,, c(2.0, 1.0, 0.5), 4600 )

  data(PoissonData)  # to see the contrasting difference
  ptspec( PoissonData, 1000, 0.5, 1.0,, c(2.0, 1.0, 0.5), 5000 )