Kalpana Kalpana (Editor)

SigSpec

Updated on
Edit
Like
Comment
Share on FacebookTweet on TwitterShare on LinkedInShare on Reddit

SigSpec is an acronym of "SIGnificance SPECtrum" and addresses a statistical technique to provide the reliability of periodicities in a measured (noisy and not necessarily equidistant) time series. It relies on the amplitude spectrum obtained by the Discrete Fourier transform (DFT) and assigns a quantity called the spectral significance (frequently abbreviated by “sig”) to each amplitude. This quantity is a logarithmic measure of the probability that the given amplitude level is due to white noise, in the sense of a type I error. It represents the answer to the question, “What would be the chance to obtain an amplitude like the measured one or higher, if the analysed time series were random?”

Contents

SigSpec may be considered a formal extension to the Lomb-Scargle periodogram, appropriately incorporating a time series to be averaged to zero before applying the DFT, which is done in many practical applications. When a zero-mean corrected dataset has to be statistically compared to a random sample, the sample mean (rather than the population mean only) has to be zero.

Probability density function (pdf) of white noise in Fourier space

Considering a time series to be represented by a set of K pairs ( t k , x k ) , the amplitude pdf of white noise in Fourier space, depending on frequency and phase angle may be described in terms of three parameters, α 0 , β 0 , θ 0 , defining the “sampling profile”, according to

tan ⁡ 2 θ 0 = K ∑ k = 0 K − 1 sin ⁡ 2 ω t k − 2 ( ∑ k = 0 K − 1 cos ⁡ ω t k ) ( ∑ k = 0 K − 1 sin ⁡ ω t k ) K ∑ k = 0 K − 1 cos ⁡ 2 ω t k − ( ∑ k = 0 K − 1 cos ⁡ ω t k ) 2 + ( ∑ k = 0 K − 1 sin ⁡ ω t k ) 2 , α 0 = 2 K 2 ( K ∑ k = 0 K − 1 cos 2 ⁡ ( ω t k − θ 0 ) − [ ∑ l = 0 K − 1 cos ⁡ ( ω t k − θ 0 ) ] 2 ) , β 0 = 2 K 2 ( K ∑ k = 0 K − 1 sin 2 ⁡ ( ω t k − θ 0 ) − [ ∑ l = 0 K − 1 sin ⁡ ( ω t k − θ 0 ) ] 2 ) .

In terms of the phase angle in Fourier space, θ , with

tan ⁡ θ = ∑ k = 0 K − 1 sin ⁡ ω t k ∑ k = 0 K − 1 cos ⁡ ω t k ,

the probability density of amplitudes is given by

ϕ ( A ) = K A ⋅ sock 2 < x 2 > exp ⁡ ( − K A 2 4 < x 2 > ⋅ sock ) ,

where the sock function is defined by

sock ⁡ ( ω , θ ) = [ cos 2 ⁡ ( θ − θ 0 ) α 0 2 + sin 2 ⁡ ( θ − θ 0 ) β 0 2 ]

and < x 2 > denotes the variance of the dependent variable x k .

False-alarm probability and spectral significance

Integration of the pdf yields the false-alarm probability that white noise in the time domain produces an amplitude of at least A ,

Φ FA ( A ) = exp ⁡ ( − K A 2 4 < x 2 > ⋅ sock ) .

The sig is defined as the negative logarithm of the false-alarm probability and evaluates to

sig ⁡ ( A ) = K A 2 log ⁡ e 4 < x 2 > ⋅ sock .

It returns the number of random time series one would have to examine to obtain one amplitude exceeding A at the given frequency and phase.

Applications

SigSpec is primarily used in asteroseismology to identify variable stars and to classify stellar pulsation (see references below). The fact that this method incorporates the properties of the time-domain sampling appropriately makes it a valuable tool for typical astronomical measurements containing data gaps.

References

SigSpec Wikipedia
(Text) CC BY-SA