We present a multilinear analysis to determine the significant predictors of star formation in galaxies using the combined EDGE-CALIFA sample of galaxies. We analyze 1845 kpc-scale lines of sight across 39 galaxies with molecular line emission measurements from EDGE combined with optical IFU data drawn from CALIFA. We use the Least Absolute Shrinkage and Selection Operator (LASSO) to identify significant factors in predicting star formation rates. We find that the local star formation rate surface density is increased by higher molecular gas surface densities and stellar surface densities. In contrast, we see lower star formation rates in systems with older stellar populations, higher gas- and stellar-phase metallicities and larger galaxy masses. We also find a significant increase in star formation rate with galactocentric radius normalized by the disk scale length, which suggests additional parameters regulating star formation rate not explored in this study.

We present TurbuStat (v1.0): a Python package for computing turbulence statistics in spectral-line data cubes. TurbuStat includes implementations of fourteen methods for recovering turbulent properties from observational data. Additional features of the software include: distance metrics for comparing two data sets; a segmented linear model for fitting lines with a break-point; a two-dimensional elliptical power-law model; multi-core fast-fourier-transform support; a suite for producing simulated observations of fractional Brownian Motion fields, including two-dimensional images and optically-thin HI data cubes; and functions for creating realistic world coordinate system information for synthetic observations. This paper summarizes the TurbuStat package and provides representative examples using several different methods. TurbuStat is an open-source package and we welcome community feedback and contributions.

We present a statistical framework to compare spectral-line data cubes of molecular clouds and use the framework to perform an analysis of various statistical tools developed from methods proposed in the literature. We test whether our methods are sensitive to changes in the underlying physical properties of the clouds or whether their behaviour is governed by random fluctuations. We perform a set of 32 self-gravitating magnetohydrodynamic simulations that test all combinations of five physical parameters -- Mach number, plasma parameter, virial parameter, driving scales, and solenoidal driving fraction -- each of which can be set to a low or high value. We create mock observational data sets of ${\rm ^{13}CO}$(1-0) emission from each simulation. We compare these mock data to a those generated from a set of baseline simulations using pseudo-distance metrics based on 18 different statistical techniques that have previously been used to study molecular clouds. We analyze these results using methods from the statistical field of experimental design and find that several of the statistics can reliably track changes in the underlying physics. Our analysis shows that the interactions between parameters are often among the most significant effects. A small fraction of statistics are also sensitive to changes in magnetic field properties. We use this framework to compare the set of simulations to observations of three nearby star-forming regions: NGC 1333, Oph A, and IC 348. We find that no one simulation agrees significantly better with the observations, although there is evidence that the high Mach number simulations are more consistent with the observations.

We propose a new approach to comparing simulated observations that enables us to determine the significance of the underlying physical effects. We utilize the methodology of experimental design, a subfield of statistical analysis, to establish a framework for comparing simulated position-position-velocity data cubes to each other. We propose three similarity metrics based on methods described in the literature: principal component analysis, the spectral correlation function, and the Cramer multi-variate two sample similarity statistic. Using these metrics, we intercompare a suite of mock observational data of molecular clouds generated from magnetohydrodynamic simulations with varying physical conditions. Using this framework, we show that all three metrics are sensitive to changing Mach number and temperature in the simulation sets, but cannot detect changes in magnetic field strength and initial velocity spectrum. We highlight the shortcomings of one-factor-at-a-time designs commonly used in astrophysics and propose fractional factorial designs as a means to rigorously examine the effects of changing physical properties while minimizing the investment of computational resources.