Coding Projects

As an astronomer, I spend a lot of time programming. While not all of my code is suitable for human consumption (i.e., not all is documented), below are a few of the python scripts I have written for general use. Source code for these projects is publicly available on my GitHub page under the GNU/GPL v3.

OutLines

Determines the emission, absorption, and P Cygni line profiles arising from winds, outflows, and bubbles under the Sobolev approximation following ( Flury 2025). Can be used for modeling observed spectral features or for predictive purposes.

towrecan

Characterizes the orthogonal dispersion in two-dimensional data sets using a combination of minimum Euclidean distance and principal component analysis. Includes error propagation.

heatcool

Computes heating and cooling rates of H II regions and the temperature solution for thermal balance under a variety of conditions. Includes new, state-of-the-art cooling functions for equilibrium and non-equilibrium scenarios.

kendall

Calculates the non-parameteric correlation coefficient Kendall's tau, accounting for censoring (upper limits). Also includes assessment of effects of outliers and uncertainties on the correlation estimator by bootstrapping and Monte Carlo simulation.

LyCsurv

Predicts Lyman continuum (LyC) luminosity or escape fraction using various multivariate survival analysis techniques, accounting for upper limits on LyC in cases of non-detections. Reference data set is the set of z~0.3 galaxies from the combined Low-redshift Lyman Continuum Survey ( LzLCS; Flury et al. 2022a).

KaplanMeier

Calculates and evaluates non-parametric survival functions for data with censoring (lower/upper limits) to determine if a measurement is consistent with a reference population / data set. Includes censored KS-style tests.

FeldCous

Handles robust statistics of low count sources, including detection probabilities, confidence intervals, and sampling of the highly non-Gaussian probability distributions in cases where count signals are exceptionally faint.

histogram

Solves for the "true" histogram or population fraction (ratio of a subset to the total distribution) for a dataset by accounting for uncertainties in the data. Each histogram bin is therefore treated as a Poisson binomial variate.