Statistics

Along with acting as the base classes for the distance metrics, the statistics can also return useful information about single datasets.

Distance Metrics

The distance metrics are computed using certain outputs contained in the related statistics class.

Nearly all of the distance metrics are actually “pseudo” - distance metrics. They must have the following properties:

  1. \(d(A, A) = 0\)
  2. Symmetric \(d(A, B) = d(B, A)\)
  3. Triangle Inequality \(d(A, B) \leq d(A, C) + d(B, C)\)

Here \(A\) and \(B\) represent two datasets (either a PPV datacube, column density map, or an associated moment of the cube).

For two datasets with different physical properties, a good statistic will return a large value \(d(A, B) \gg 0\). If the datasets have similar physical properties, the distance should be small \(d(A, B) \approx 0\).

Additionally, the statistics should ideally be insensitive to spatial shifts \(d\left( A\left[ x,y,v \right], A\left[ x+\delta x,y,v \right] \right)=0\) and independent of the noise level (for observational data) \(d\left( A + \mathcal{N}\left(0, \sigma_1^2 \right), A + \mathcal{N}\left(0, \sigma_2^2 \right) \right) \approx 0\).

Source Code

Functions

WidthEstimate1D(inList[, method]) Find widths from spectral eigenvectors.
WidthEstimate2D(inList[, method, noise_ACF, ...]) Estimate spatial widths from a set of autocorrelation images.
stats_wrapper(dataset1, dataset2[, ...]) Function to run all of the statistics on two datasets.

Classes

BiSpectrum(img) Computes the bispectrum (three-point correlation function) of the given image (Burkhart et al., 2010).
BiSpectrum_Distance(data1, data2[, ...]) Calculate the distance between two images based on their bicoherence.
Cramer_Distance(cube1, cube2[, ...]) Compute the Cramer distance between two data cubes.
DeltaVariance(img[, header, weights, ...]) The delta-variance technique as described in Ossenkopf et al.
DeltaVariance_Distance(dataset1, dataset2[, ...]) Compares 2 datasets using delta-variance.
DendroDistance(cube1, cube2[, min_deltas, ...]) Calculate the distance between 2 cubes using dendrograms.
Dendrogram_Stats(data[, header, min_deltas, ...]) Dendrogram statistics as described in Burkhart et al.
Genus(img[, lowdens_percent, ...]) Genus Statistics based off of Chepurnov et al.
GenusDistance(img1, img2[, smoothing_radii, ...]) Distance Metric for the Genus Statistic.
MVC(centroid, moment0, linewidth[, header]) Implementation of Modified Velocity Centroids (Lazarian & Esquivel, 03)
MVC_Distance(data1, data2[, fiducial_model, ...]) Distance metric for MVC.
Mahalanobis(cube) docstring for Mahalanobis
Mahalanobis_Distance(cube1, cube2) docstring for Mahalanobis_Distance
PCA(cube[, n_eigs, distance]) Implementation of Principal Component Analysis (Heyer & Brunt, 2002)
PCA_Distance(cube1, cube2[, n_eigs, ...]) Compare two data cubes based on the eigenvalues of the PCA decomposition.
PDF(img[, min_val, bins, weights, ...]) Create the PDF of a given array.
PDF_Distance(img1, img2[, min_val1, ...]) Calculate the distance between two arrays using their PDFs.
PSpec_Distance(data1, data2[, weights1, ...]) Distance metric for the spatial power spectrum.
PowerSpectrum(img[, header, weights]) Compute the power spectrum of a given image.
SCF(cube[, header, size, roll_lags]) Computes the Spectral Correlation Function of a data cube (Rosolowsky et al, 1999).
SCF_Distance(cube1, cube2[, size, boundary, ...]) Calculates the distance between two data cubes based on their SCF surfaces.
StatMoments(img[, weights, radius, ...]) Statistical Moments of a given image are returned.
StatMoments_Distance(image1, image2[, ...]) Compute the distance between two images based on their moments.
Tsallis(img[, lags, num_bins, periodic]) The Tsallis Distribution (see Tofflemire et al., 2011)
Tsallis_Distance(array1, array2[, lags, ...]) Distance Metric for the Tsallis Distribution.
VCA(cube[, header, slice_size]) The VCA technique (Lazarian & Pogosyan, 2004).
VCA_Distance(cube1, cube2[, slice_size, ...]) Calculate the distance between two cubes using VCA.
VCS(cube[, header, vel_units]) The VCS technique (Lazarian & Pogosyan, 2004).
VCS_Distance(cube1, cube2[, breaks, ...]) Calculate the distance between two cubes using VCS.
Wavelet(data[, header, scales, num, ...]) Compute the wavelet transform of a 2D array.
Wavelet_Distance(dataset1, dataset2[, ...]) Compute the distance between the two cubes using the Wavelet transform.
class turbustat.statistics.base_pspec2.StatisticBase_PSpec2D[source] [edit on github]

Common features shared by 2D power spectrum methods.