ptspec: The Periodogram of Point Process Data In SAPP: Statistical Analysis of Point Processes

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

 `1` ```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

 ```1 2 3 4 5 6 7``` ```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 July 2, 2020, 2:59 a.m.