Astron. Astrophys. Suppl. Ser. 132, 107-119

The art of fitting p-mode spectra

I. Maximum likelihood estimation

T. Appourchaux¹ - L. Gizon^1,2 - M.-C. Rabello-Soares^1,3

Send offprint request: thierrya@so.estec.esa.nl

1 - Space Science Department of ESA, ESTEC, NL-2200 AG Noordwijk, The Netherlands
2 - W.W. Hansen Experimental Physics Laboratory, Center for Space Science and Astrophysics, Stanford University, Stanford, CA 94305-4085, U.S.A.
3 - Teoretisk Astrofysik Center, Århus Universitet, DK-8000 Århus C, Denmark

Received October 17, 1997; accepted March 24, 1998

Abstract:

In this article we present our state of the art of fitting helioseismic p-mode spectra. We give a step by step recipe for fitting the spectra: statistics of the spectra both for spatially unresolved and resolved data, the use of Maximum Likelihood estimates, the statistics of the p-mode parameters, the use of Monte-Carlo simulation and the significance of fitted parameters. The recipe is applied to synthetic low-resolution data, similar to those of the LOI, using Monte-Carlo simulations. For such spatially resolved data, the statistics of the Fourier spectrum is assumed to be a multi-normal distribution; the statistics of the power spectrum is not a $\chi^{2}$ with 2 degrees of freedom. Results for l=1 shows that all parameters describing the p modes can be obtained with negligible bias and with minimum variance provided that the leakage matrix is known. Systematic errors due to an imperfect knowledge of the leakage matrix are derived for all the p-mode parameters.

Key words: methods: analytical; data analysis; statistical -- Sun: oscillations