1 / 9

Comments on John Rice’s paper “On Detecting Periodicity in Astronomical Point Processes”

Comments on John Rice’s paper “On Detecting Periodicity in Astronomical Point Processes” Jeff Scargle. This paper should be quite accessible to astronomers: Clear description of a well-defined, practical problem Does not go off into asymptotia Choices for the analyst not overwhelming

marli
Download Presentation

Comments on John Rice’s paper “On Detecting Periodicity in Astronomical Point Processes”

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Comments on John Rice’s paper “On Detecting Periodicity in Astronomical Point Processes” Jeff Scargle

  2. This paper should be quite accessible to astronomers: • Clear description of a well-defined, practical problem • Does not go off into asymptotia • Choices for the analyst not overwhelming • Nice description of problem of unspecified alternative -- geometrical view: tests good in restricted sets of directions in parameter space • Ideas have direct parallels with Bayesian methods • [Bretthorst periodogram is similar to Rice’s power formula for the case of fundamental only (no harnomics; Rayleigh test)]

  3. Sampling in frequency: In periodogram/spectral analysis with evenly spaced data, one normally evaluates the power spectrum at the Fourier frequencies fn = n / T n = 0, 1, 2, … (N/2) - 1 But often one is tempted to oversample, to avoid “missing a real peak” For uneven sampling or for point data as Rice considers, the Nyquist frequency is not well defined and it is not obvious how many and which frequencies should be sampled. He proposes a novel and potentially extremely useful procedure to integrate over finite frequency bands. One needs to carefully assess the statistical significance of any peak found this way. Rice’s analysis of this is a good beginning, but needs more development … examples, etc., the usual stuff astronomers need to feel good about a method.

  4. Incorporating Frequency Drift • Simple extension of model to include frequency derivative: • Power = P( f, df/dt ) • Perhaps this can be interpreted as a time-frequency plot: • change of variables, such as • ( f , df/dt ) ? (f true , t ) where t = ( f true – f ) / df/dt • This would be an excellent tool for astronomers, for not just pulsars, but for tracking QPOs in accreting black hole/neutron star systems, etc.!: • Photons in • P( f, t ) out • knob tunes time-frequency resolution trade-off

  5. Alice Harding

  6. Alice Harding

  7. Bayesian Estimation of Time Series Lags and Structure , J. Scargle, in Bayesian Inference and Maximum Entropy in Science and Engineering, 2001. AIP Conference Proceedings, ed. Robert L. Fry. http://proceedings.aip.org/ (volume617),

  8. Bayesian Time Series Tools Periodogram sin( ωt + a ) Pr(ω) ~ exp( - P(ω) / σ2 ) Bretthorst period/phase Cross-corr. f1(t), f2( t – τ ) Pr(τ) ~ exp( - ρ(τ) / const ) Scargle Time Lag Scalegram s = scale Pr(s) ~ exp( - W(s) / const ) Scargle, Loredo Time-Freq. Bayesianization of Rice’s frequency drift tracker? Distribution Others ?

More Related