Related papers: Multi-scale time- and frequency-domain likelihood …
The field of gamma ray astronomy relies heavily on the statistical analysis of data. Because of the paucity of data, and the often large errors associated with detecting gamma rays, analysis and interpretation of the data require…
Due to instrumental limitations, the nature of which vary from case to case, spectroscopic galaxy redshift surveys usually do not collect redshifts for all galaxies in the population of potential targets. Especially problematic is the…
Multivariate analyses play an important role in high energy physics. Such analyses often involve performing an unbinned maximum likelihood fit of a probability density function (p.d.f.) to the data. This paper explores a variety of unbinned…
Detection of point sources in images is a fundamental operation in astrophysics, and is crucial for constraining population models of the underlying point sources or characterizing the background emission. Standard techniques fall short in…
Identifying as many gamma-ray pulsars as possible in the Fermi Large Area Telescope (LAT) data helps test pulsar emission models by comparing predicted and observed properties for a large, varied sample with as little selection bias as…
I describe ongoing work developing Bayesian methods for flexible modeling of arrival time series data without binning, aiming to improve detection and measurement of X-ray and gamma-ray pulsars, and of pulses in gamma-ray bursts. The…
Several rapid parameter estimation methods have recently been advanced to deal with the computational challenges of the problem of Bayesian inference of the properties of compact binary sources detected in the upcoming science runs of the…
I describe a new time-domain algorithm for detecting localized structures (bursts), revealing pulse shapes, and generally characterizing intensity variations. The input is raw counting data, in any of three forms: time-tagged photon events…
Frequency-filtered photon correlations have been proven to be extremely useful in grasping how the detection process alters photon statistics. Harnessing the spectral correlations also permits refinement of the emission and unraveling of…
Future ground-based and space-borne interferometric gravitational-wave detectors may capture between tens and thousands of binary coalescence events per year. There is a significant and growing body of work on the estimation of…
Wide-field gamma-ray telescopes typically have highly variable event-by-event resolution which leads to a number of unique and challenging analysis requirements -- particularly when conducting transient searches over multiple time scales.…
We investigate the use of optical photometric variability to select and identify blazars in large-scale time-domain surveys, in part to aid in the identification of blazar counterparts to the ~30% of gamma-ray sources in the Fermi 2FGL…
We present maximum-likelihood search methods for time-dependent fluxes from point sources, such as flares or periodic emissions. We describe a method for the case when the time dependence of the flux can be assumed a priori from other…
We develop two algorithms, based on maximum likelihood (ML) inference, for estimating the parameters of polarized radio sources which emit at a single rotation measure (RM), e.g., pulsars. These algorithms incorporate the flux density…
Parameter estimation of gravitational wave signals is computationally intensive and typically requires millions of likelihood evaluations to construct posterior probability distributions. This computational cost increases significantly in…
Linear time series analysis, mainly the Fourier transform based methods, has been quite successful in extracting information contained in the ever-modulating light curves (Lcs) of active galactic nuclei, and thereby contribute in…
We provide a new approach to measure power spectra and reconstruct time series in active galactic nuclei (AGNs) based on the fact that the Fourier transform of AGN stochastic variations is a series of complex Gaussian random variables. The…
A key goal of radio and $\gamma-$ray observations of active galactic nuclei is to characterize their time variability in order to elucidate physical processes responsible for the radiation. I describe algorithms for relevant time series…
Pulsar Timing Array (PTA) searches for nHz gravitational-wave backgrounds (GWBs) typically model time-correlated noise by assuming a diagonal covariance in Fourier space, neglecting inter-frequency correlations introduced by the finite…
Inferring the source properties of a gravitational wave signal has traditionally been very computationally intensive and time consuming. In recent years, several techniques have been developed that can significantly reduce the computational…