Jeff Scargle

Thanks! ... I can point you to a short and simple MatLab script that essentially updates the FORTRAN (do you remember FORTRAN?) in "Studies in astronomical time series analysis. III...

Hi ... the main difference is that it computes the full Fourier transform of arbitrarily spaced data. The absolute square is then equivalent to LS ... but the implementation is...

... oh, and you can also do the co-spectrum if you have the full FT! In this way I have implemented a bin-free co-spectrum estimator for event (time-tagged photon) data.

Hi Matteo ... good to hear from you. First: I do not check my nasa mail account very often and mostly use this gmail account ... I forget the details...

Good comments. In the restricted context of a single sinusoid don't you think Larry Bretthorst's treatment is pretty compelling? Sure, you have to choose what parameters you take as "nuisance"...

… on a different topic. In the discussion of Fig. 3 you state: "Second, notice that as the number of terms is increased, the general \background" level of the periodogram...