This page was made by George Rybicki and Bill Press, as a location for our links to the subject, and as a distribution point for our source code that implements so-called "fast" statistical methods.
We define "fast" methods as being methods that
G.B. Rybicki and W.H. Press, ``A Class of Fast Methods for Processing Irregularly Sampled or Otherwise Inhomogeneous One-Dimensional Data,'' Physical Review Letters, 74, 1060 (1995).Here is a downloadable PDF file.
The above paper is rather terse, due to the requirements of the journal. We have written related papers on statistical reconstruction (though not via fast methods) which may help you to understand the above fast-methods paper:
G.B. Rybicki and W.H. Press, ``Interpolation, Realization, and Reconstruction of Noisy, Irregularly Sampled Data,'' Astrophysical Journal, 398, 169 (1992). (Reprint in ADS.)
W.H. Press and G.B. Rybicki, ``Large-Scale Linear Methods for Interpolation, Realization, and Reconstruction of Noisy, Irregularly Sampled Data,'' in Predicting the Future and Understanding the Past: a Comparison of Approaches, A.S. Weigend and N.A. Gershenfeld, eds. (Reading, MA: Addison-Wesley, 1993).
Some Notes on Wiener Reconstruction (10 October 1997). PDF file. (Shows how to rewrite the unbiased estimator formulas in terms of the structure function alone, instead of the correlation function.)
Notes on the Ornstein-Uhlenbeck Process (2 December 1994).PDF file. (This is the process with the exponential correlation function that, when used as an approximation, enables the fast methods. Here it is looked at in its own right.)
Fast Operations Involving the Matrix |ti-tj| (4 March 1998). PDF file. (Shows how to take the limit to pure random walk processes, where the structure function is linear in time difference. Sample code is in the file lininv.f, previously listed.)
Reduction Method for Duplicated Times (3 May 1995). PDF file. (Shows how to prevent duplicate times from causing the fast-method to become ill-posed. This technique is used in some of the sample codes above.)
Estimation of Correlation Properties (10 Jan 1995). PDF file. (Shows how to estimate the parametrized correlation matrix of a process from data by a method that can be made fast using the same sorts of techniques.)
W.H. Press and G.B. Rybicki, ``Fast Algorithm for Spectral Analysis of Unevenly Sampled Data," Astrophysical Journal, 338, 277 (1989). (Reprint in ADS.)
G.B. Rybicki and D.G. Hummer, ``Fast Solution for the Diagonal Elements of the Inverse of a Tridiagonal Matrix," Appendix B of ``An Accelerated Lambda Iteration Method for Multilevel Radiative Transfer," Astronomy and Astrophysics 245, 171 (1991). Downloadable as PDF file.