next up previous contents
Next: Correlograms Up: SigSpec User's Manual by Previous: MultiSine profiles   Contents


Preview

Since the prewhitening cascade performed by SIGSPEC may be extremely time consuming, the program can compute a preview. This add-on is activated by the keyword preview in the .ini file.

Whereas the significance spectra rely on the False-Alarm Probability compared to a noise dataset with the same rms error as the given time series (or series of residuals, respectively), the significance spectrum provided in the file preview.dat represents a set of identified maxima in the significance spectrum of the original time series, but based on the point-to-point scatter in the time domain rather than on the standard deviation of observables. The lower sig limit for writing a local maximum to the file preview.dat is specified as the argument to the keyword preview in the .ini file.

The calculation of the sig is based on the assumption that only the point-to-point scatter is random, and everything else contributing to the rms error represents signal that will be prewhitened in the course of the subsequent loop. The preview output is to be considered as a rough estimate for the final result obtained by step-by-step prewhitening and contains not only the intrinsic variations but also all aliases, which will not occur in the following analysis. The file preview.dat consists of four columns referring to

  1. frequency [inverse time units],
  2. sig,
  3. DFT amplitude [units of observable],
  4. phase [rad].



Example. The sample project preview contains a preview file for the V photometry of IC4996#89. In the file preview.ini, the line

preview:siglimit 5

sets the sig threshold to 5. The output file preview/preview.dat contains 11 components, sorted by frequency. The frequencies and corresponding sigs in the first two columns are

 0.9945158494303480	  6.1674140356166323
 1.9917563998155081	  6.9302735632389876
 2.1388902515119841	  6.7175642729893710
 2.9835475482878717	  8.6773027802854656
 3.1361308018982843	  9.4899187898938777
 3.9862375005883859	  8.9589776551282210
 4.1333713522847031	  8.7492615592402885
 4.9780286490607102	  5.0523760377159039
 5.1360613045861099	  5.3438911274207790
11.0268647743572874	  5.5214237212500406
12.0241053247411784	  5.6674302270769710

Fig.17 displays the significance and amplitude spectrum of the original time series. Since the preview does not employ any prewhitenings, aliases are present in the file.

  • The signal at 3.132 cycles per day corresponds to components # 3, 5, 7, and 9.
  • The signal at 3.986 cycles per day corresponds to components # 1, 2, 4, 6, and 8.
  • The signal at 5.409 cycles per day is not found in the preview. In the result of the prewhitening sequence, its sig is 5.02. Since the sig in the preview relies on the rms deviation of the original time series, whereas the final sig is based on the rms deviation after the previous prewhitening step, the sig associated to this frequency falls below the pre-selected threshold of 5 in the preview. The significance spectrum (grey line in the left panel of Fig.17) shows a peak at the frequency under consideration the sig of which is $\approx$ 4.8.
  • Components # 10 and 11 are 1-cycle-per-day aliases of each other, but do not show up in the final result, preview/result.dat.

Figure 17: Grey: Fourier spectra for the sample project preview. Left: significance spectrum. Right: DFT amplitudes. The significant components in the preview are indicated by dots with dashed drop lines (file preview/preview.dat). The default sig threshold of 5 is represented by a horizontal dashed line in the left panel.
\includegraphics[clip,angle=0,width=109mm, clip]{eps/preview.eps}


next up previous contents
Next: Correlograms Up: SigSpec User's Manual by Previous: MultiSine profiles   Contents
Piet Reegen 2009-09-23