Phase dispersion minimization

Last updated
PDM2 analysis for the double mode Cepheid variable TU Cas, primary period Tu cas.jpg
PDM2 analysis for the double mode Cepheid variable TU Cas, primary period
Raw data set for the analysis shown above. The large gaps are often found in ground-based observations. TUCasData.jpg
Raw data set for the analysis shown above. The large gaps are often found in ground-based observations.

Phase dispersion minimization (PDM) is a data analysis technique that searches for periodic components of a time series data set. It is useful for data sets with gaps, non-sinusoidal variations, poor time coverage or other problems that would make Fourier techniques unusable. It was first developed by Stellingwerf in 1978 [1] and has been widely used for astronomical and other types of periodic data analyses. Source code is available for PDM analysis. The current version of this application is available for download. [2]

Contents

Background

PDM is a variant of a standard astronomical technique called data folding. This involves guessing a trial period for the data, and cutting, or "folding" the data into multiple sub-series with a time duration equal to the trial period. The data are now plotted versus "phase", or a scale of 0->1, relative to the trial period. If the data is truly periodic with this period a clean functional variation, or "light curve", will emerge. If not the points will be randomly distributed in amplitude.

As early as 1926 Whittiker and Robinson [3] proposed an analysis technique of this type based on maximizing the amplitude of the mean curve. Another technique focusing on the variation of data at adjacent phases was proposed in 1964 by Lafler and Kinman. [4] Both techniques had difficulties, particularly in estimating the significance of a possible solution.

PDM analysis

PDM divides the folded data into a series of bins and computes the variance of the amplitude within each bin. The bins can overlap to improve phase coverage, if needed. The bin variances are combined and compared to the overall variance of the data set. For a true period the ratio of the bin to the total variances will be small. For a false period the ratio will be approximately unity. A plot of this ratio versus trial period will usually indicate the best candidates for periodic components. Analyses of the statistical properties of this approach have been given by Nemec & Nemec [5] and Schwarzenberg-Czerny. [6]

PDM2 updates

The original PDM technique has been updated (PDM2) in several areas::

PDM2b analysis of RR Lyrae rich data set. Mean curve is in red, using 100 bins and a Spline fit. RR Lyrae, rich data example.jpg
PDM2b analysis of RR Lyrae rich data set. Mean curve is in red, using 100 bins and a Spline fit.

See reference (2) for a detailed technical discussion, test cases, C source code, and a Windows application package.

Binless PDM

In Plavchan et al. 2008, [7] Plavchan introduced a binless version of the phase dispersion minimization algorithm. The algorithm was further revised in 2014 in Parks, Plavchan et al. 2014, [8] and is available for highly-parallel use online at the NASA Exoplanet Archive. [9] The binned PDM approach is susceptible to period aliases when the cadence is semi-regular (e.g., nightly observations of a star brightness). Plavchan and colleagues avoided this aliasing by computing a box-car smoothed phased time-series, where the box-car width can be thought of as the old bin size. The original folded time-series is compared to the smoothed time-series, and the best period is found when the time-series are most similar. See the NASA Exoplanet Archive for more information on statistical significance and approaches.

Related Research Articles

Analysis of variance (ANOVA) is a collection of statistical models and their associated estimation procedures used to analyze the differences among means. ANOVA was developed by the statistician Ronald Fisher. ANOVA is based on the law of total variance, where the observed variance in a particular variable is partitioned into components attributable to different sources of variation. In its simplest form, ANOVA provides a statistical test of whether two or more population means are equal, and therefore generalizes the t-test beyond two means. In other words, the ANOVA is used to test the difference between two or more means.

<span class="mw-page-title-main">Lyra</span> Constellation in the northern celestial hemisphere

Lyra is a small constellation. It is one of the 48 listed by the 2nd century astronomer Ptolemy, and is one of the modern 88 constellations recognized by the International Astronomical Union. Lyra was often represented on star maps as a vulture or an eagle carrying a lyre, and hence is sometimes referred to as Vultur Cadens or Aquila Cadens, respectively. Beginning at the north, Lyra is bordered by Draco, Hercules, Vulpecula, and Cygnus. Lyra is nearly overhead in temperate northern latitudes shortly after midnight at the start of summer. From the equator to about the 40th parallel south it is visible low in the northern sky during the same months.

<span class="mw-page-title-main">Variable star</span> Star whose brightness as seen from Earth fluctuates

A variable star is a star whose brightness as seen from Earth changes with time. This variation may be caused by a change in emitted light or by something partly blocking the light, so variable stars are classified as either:

<span class="mw-page-title-main">Light curve</span> Graph of light intensity of a celestial object or region, as a function of time

In astronomy, a light curve is a graph of the light intensity of a celestial object or region as a function of time, typically with the magnitude of light received on the y-axis and with time on the x-axis. The light is usually in a particular frequency interval or band.

<span class="mw-page-title-main">Cepheid variable</span> Type of variable star that pulsates radially

A Cepheid variable is a type of variable star that pulsates radially, varying in both diameter and temperature. It changes in brightness, with a well-defined stable period and amplitude.

<span class="mw-page-title-main">RR Lyrae variable</span> Type of variable star

RR Lyrae variables are periodic variable stars, commonly found in globular clusters. They are used as standard candles to measure (extra) galactic distances, assisting with the cosmic distance ladder. This class is named after the prototype and brightest example, RR Lyrae.

PDM may stand for:

Linear trend estimation is a statistical technique to aid interpretation of data. When a series of measurements of a process are treated as, for example, a sequences or time series, trend estimation can be used to make and justify statements about tendencies in the data, by relating the measurements to the times at which they occurred. This model can then be used to describe the behaviour of the observed data, without explaining it.

A test statistic is a statistic used in statistical hypothesis testing. A hypothesis test is typically specified in terms of a test statistic, considered as a numerical summary of a data-set that reduces the data to one value that can be used to perform the hypothesis test. In general, a test statistic is selected or defined in such a way as to quantify, within observed data, behaviours that would distinguish the null from the alternative hypothesis, where such an alternative is prescribed, or that would characterize the null hypothesis if there is no explicitly stated alternative hypothesis.

<span class="mw-page-title-main">V399 Carinae</span> Star in the constellation Carina

V399 Carinae is a variable star in the constellation Carina.

In neurophysiology, peristimulus time histogram and poststimulus time histogram, both abbreviated PSTH or PST histogram, are histograms of the times at which neurons fire. It is also sometimes called pre event time histogram or PETH. These histograms are used to visualize the rate and timing of neuronal spike discharges in relation to an external stimulus or event. The peristimulus time histogram is sometimes called perievent time histogram, and post-stimulus and peri-stimulus are often hyphenated.

<span class="mw-page-title-main">Least-squares spectral analysis</span> Periodicity computation method

Least-squares spectral analysis (LSSA) is a method of estimating a frequency spectrum based on a least-squares fit of sinusoids to data samples, similar to Fourier analysis. Fourier analysis, the most used spectral method in science, generally boosts long-periodic noise in the long and gapped records; LSSA mitigates such problems. Unlike in Fourier analysis, data need not be equally spaced to use LSSA.

In probability theory and statistics, the index of dispersion, dispersion index,coefficient of dispersion,relative variance, or variance-to-mean ratio (VMR), like the coefficient of variation, is a normalized measure of the dispersion of a probability distribution: it is a measure used to quantify whether a set of observed occurrences are clustered or dispersed compared to a standard statistical model.

<span class="mw-page-title-main">RR Lyrae</span> Star in the constellation Lyra

RR Lyrae is a variable star in the Lyra constellation, figuring in its west near to Cygnus. As the brightest star in its class, it became the eponym for the RR Lyrae variable class of stars and it has been extensively studied by astronomers. RR Lyrae variables serve as important standard candles that are used to measure astronomical distances. The period of pulsation of an RR Lyrae variable depends on its mass, luminosity and temperature, while the difference between the measured luminosity and the actual luminosity allows its distance to be determined via the inverse-square law. Hence, understanding the period-luminosity relation for a local set of such stars allows the distance of more distant stars of this type to be determined.

SigSpec is a statistical technique to provide the reliability of periodicities in a measured time series. It relies on the amplitude spectrum obtained by the Discrete Fourier transform (DFT) and assigns a quantity called the spectral significance to each amplitude. This quantity is a logarithmic measure of the probability that the given amplitude level would be seen in 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?”

<span class="mw-page-title-main">Stellar pulsation</span>

Stellar pulsations are caused by expansions and contractions in the outer layers as a star seeks to maintain equilibrium. These fluctuations in stellar radius cause corresponding changes in the luminosity of the star. Astronomers are able to deduce this mechanism by measuring the spectrum and observing the Doppler effect. Many intrinsic variable stars that pulsate with large amplitudes, such as the classical Cepheids, RR Lyrae stars and large-amplitude Delta Scuti stars show regular light curves.

<span class="mw-page-title-main">Surface wave inversion</span>

Seismic inversion involves the set of methods which seismologists use to infer properties through physical measurements. Surface-wave inversion is the method by which elastic properties, density, and thickness of layers in the subsurface are obtained through analysis of surface-wave dispersion. The entire inversion process requires the gathering of seismic data, the creation of dispersion curves, and finally the inference of subsurface properties.

<span class="mw-page-title-main">V473 Lyrae</span> Star in the constellation Lyra

V473 Lyrae is a variable star in the constellation Lyra. It is an unusual Classical Cepheid variable with a visual range of 5.99 to 6.35.

<span class="mw-page-title-main">DY Pegasi</span> Star in the constellation Pegasus

DY Pegasi, abbreviated DY Peg, is a binary star system in the northern constellation of Pegasus. It is a well-studied SX Phoenicis variable star with a brightness that ranges from an apparent visual magnitude of 9.95 down to 10.62 with a period of 1.75 hours. This system is much too faint to be seen with the naked eye, but can be viewed with large binoculars or a telescope. Based on its high space motion and low abundances of heavier elements, it is a population II star system.

<span class="mw-page-title-main">TU Ursae Majoris</span>

TU Ursae Majoris is a variable star in the northern circumpolar constellation of Ursa Major. It is classified as a Bailey-type 'ab' RR Lyrae variable with a period of 0.557648 days that ranges in brightness from apparent visual magnitude of 9.26 down to 10.24. The distance to this star is approximately 2,090 light years based on parallax measurements. It is located near the north galactic pole at a distance that indicates this is a member of the galactic halo.

References

  1. "Period Determination Using Phase Dispersion Minimization", Stellingwerf, R.F., Astrophysical.J. v224, p953, 1978.
  2. "PDM2 Application, Technical Manual, and test data sets", Stellingwerf, R. F., 2006.
  3. "The Calculus of Observations", Whittiker, E. T., Robinson, G. (London: Blackie and Son) 1926.
  4. "An RR Lyrae Star Survey with Ihe Lick 20-INCH Astrograph II. The Calculation of RR Lyrae Periods by Electronic Computer", Lafler, J., Kinman, T. D. Astrophysical J., v11, p216, 1965.
  5. "A test of significance for periods derived using phase-dispersion-minimization techniques," Nemec & Nemec, Astronomical.J. v90, p2317, 1985.
  6. "The Correct Probability Distribution for the Phase Dispersion Minimization Periodogram", Schwarzenberg-Czerny, A., Astrophysical J. v489, p941, 1997.
  7. Plavchan, Peter; Jura, M.; Kirkpatrick, J. Davy; Cutri, Roc M.; Gallagher, S. C. (2008). "Near-Infrared Variability in the 2MASS Calibration Fields: A Search for Planetary Transit Candidates". The Astrophysical Journal Supplement Series. 175 (1): 191–228. arXiv: 0709.1182 . Bibcode:2008ApJS..175..191P. doi:10.1086/523644. S2CID   30540669.
  8. Parks, J. Robert; Plavchan, Peter; White, Russel J.; Gee, Alan H. (2014). "Periodic and Aperiodic Variability in the Molecular Cloud rho Ophiuchus". The Astrophysical Journal Supplement Series. 211 (1): 3. arXiv: 1309.5300 . Bibcode:2014ApJS..211....3P. doi:10.1088/0067-0049/211/1/3. S2CID   51438707.
  9. http://exoplanetarchive.ipac.caltech.edu/