fdarvas

fdarvas / GregoryLoredoAlgorithm / 0.2.0

README.md

Computes the Gregory-Loredo algorithm on a list of arrival times

This function computes the likelihood of a set of arrival times originating from a periodic process rather than constant rate (Poisson) process (e.g. background noise).

based on Gregory, P. C. and Thomas. J. Loredo, 1992, "A New Method For The Detection Of A Periodic Signal Of Unknown Shape And Period" in The Astrophysical Journal, Astrophysical J., 398, p.146

Inputs

  • Tlist (required)- a list of arrival times
  • m_max  (optional, default m_max=12) - number of phase bins
  • w_range (optional, default min(20,N/10)*pi/T, N= number of arrival times, T= max(Tlist)) - frequency range to scan for periodicity
  • ni (optional, default ni=10) - number of bins for numerical integration of the phase from 0:2pi/m
  • parallel (optional, default parallel=False) - parallel execution flag

Outputs

  • O_period - odds ratio for periodic process
  • p_period - probability of periodic process (0<=p_period<=1)
  • m_opt - optimal bin size
  • S - spectrum for m_opt
  • w - frequency range of S
  • w_peak - most likely frequency
  • w_mean - mean frequency
  • w_conf - 95% confidence interval of frequency  

Other sample inputs: 

{"Tlist" : [80,921,2281,5831,6095,7703,8205,8761,11010,11095,14201,22531], "m_max" : 5}

{"Tlist" : [80,921,2281,5831,6095,7703,8205,8761,11010,11095,14201,22531],"m_max" :5, "ni" : 20}