Source author record

Jorge F. Silva

Jorge F. Silva appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.IM astro-ph.SR Information Theory Machine Learning math.IT physics.data-an physics.comp-ph physics.ins-det

Catalog footprint

What is connected

6works

8topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2022arXiv

Bayesian inference in single-line spectroscopic binaries with a visual orbit

We present a Bayesian inference methodology for the estimation of orbital parameters on single-line spectroscopic binaries with astrometric data, based on the No-U-Turn sampler Markov chain Monte Carlo algorithm. Our approach is designed to provide a precise and efficient estimation of the joint posterior distribution of the orbital parameters in the presence of partial and heterogeneous observations. This scheme allows us to directly incorporate prior information about the system - in the form of a trigonometric parallax, and an estimation of the mass of the primary component from its spectral type - to constrain the range of solutions, and to estimate orbital parameters that cannot be usually determined (e.g. the individual component masses), due to the lack of observations or imprecise measurements. Our methodology is tested by analyzing the posterior distributions of well-studied double-line spectroscopic binaries treated as single-line binaries by omitting the radial velocity data of the secondary object. Our results show that the system's mass ratio can be estimated with an uncertainty smaller than 10% using our approach. As a proof of concept, the proposed methodology is applied to twelve single-line spectroscopic binaries with astrometric data that lacked a joint astrometric-spectroscopic solution, for which we provide full orbital elements. Our sample-based methodology allows us also to study the impact of different posterior distributions in the corresponding observations space. This novel analysis provides a better understanding of the effect of the different sources of information on the shape and uncertainty of the orbit and radial velocity curve.

preprint2021arXiv

Data-Driven Representations for Testing Independence: Modeling, Analysis and Connection with Mutual Information Estimation

This work addresses testing the independence of two continuous and finite-dimensional random variables from the design of a data-driven partition. The empirical log-likelihood statistic is adopted to approximate the sufficient statistics of an oracle test against independence (that knows the two hypotheses). It is shown that approximating the sufficient statistics of the oracle test offers a learning criterion for designing a data-driven partition that connects with the problem of mutual information estimation. Applying these ideas in the context of a data-dependent tree-structured partition (TSP), we derive conditions on the TSP's parameters to achieve a strongly consistent distribution-free test of independence over the family of probabilities equipped with a density. Complementing this result, we present finite-length results that show our TSP scheme's capacity to detect the scenario of independence structurally with the data-driven partition as well as new sampling complexity bounds for this detection. Finally, some experimental analyses provide evidence regarding our scheme's advantage for testing independence compared with some strategies that do not use data-driven representations.

preprint2021arXiv

Studying the Interplay between Information Loss and Operation Loss in Representations for Classification

Information-theoretic measures have been widely adopted in the design of features for learning and decision problems. Inspired by this, we look at the relationship between i) a weak form of information loss in the Shannon sense and ii) the operation loss in the minimum probability of error (MPE) sense when considering a family of lossy continuous representations (features) of a continuous observation. We present several results that shed light on this interplay. Our first result offers a lower bound on a weak form of information loss as a function of its respective operation loss when adopting a discrete lossy representation (quantization) instead of the original raw observation. From this, our main result shows that a specific form of vanishing information loss (a weak notion of asymptotic informational sufficiency) implies a vanishing MPE loss (or asymptotic operational sufficiency) when considering a general family of lossy continuous representations. Our theoretical findings support the observation that the selection of feature representations that attempt to capture informational sufficiency is appropriate for learning, but this selection is a rather conservative design principle if the intended goal is achieving MPE in classification. Supporting this last point, and under some structural conditions, we show that it is possible to adopt an alternative notion of informational sufficiency (strictly weaker than pure sufficiency in the mutual information sense) to achieve operational sufficiency in learning.

preprint2016arXiv

Analysis of the Bayesian Cramer-Rao lower bound in astrometry: Studying the impact of prior information in the location of an object

Context. The best precision that can be achieved to estimate the location of a stellar-like object is a topic of permanent interest in the astrometric community. Aims. We analyse bounds for the best position estimation of a stellar-like object on a CCD detector array in a Bayesian setting where the position is unknown, but where we have access to a prior distribution. In contrast to a parametric setting where we estimate a parameter from observations, the Bayesian approach estimates a random object (i.e., the position is a random variable) from observations that are statistically dependent on the position. Methods. We characterize the Bayesian Cramer-Rao (CR) that bounds the minimum mean square error (MMSE) of the best estimator of the position of a point source on a linear CCD-like detector, as a function of the properties of detector, the source, and the background. Results. We quantify and analyse the increase in astrometric performance from the use of a prior distribution of the object position, which is not available in the classical parametric setting. This gain is shown to be significant for various observational regimes, in particular in the case of faint objects or when the observations are taken under poor conditions. Furthermore, we present numerical evidence that the MMSE estimator of this problem tightly achieves the Bayesian CR bound. This is a remarkable result, demonstrating that all the performance gains presented in our analysis can be achieved with the MMSE estimator. Conclusions The Bayesian CR bound can be used as a benchmark indicator of the expected maximum positional precision of a set of astrometric measurements in which prior information can be incorporated. This bound can be achieved through the conditional mean estimator, in contrast to the parametric case where no unbiased estimator precisely reaches the CR bound.

preprint2015arXiv

Performance analysis of the Least-Squares estimator in Astrometry

We characterize the performance of the widely-used least-squares estimator in astrometry in terms of a comparison with the Cramer-Rao lower variance bound. In this inference context the performance of the least-squares estimator does not offer a closed-form expression, but a new result is presented (Theorem 1) where both the bias and the mean-square-error of the least-squares estimator are bounded and approximated analytically, in the latter case in terms of a nominal value and an interval around it. From the predicted nominal value we analyze how efficient is the least-squares estimator in comparison with the minimum variance Cramer-Rao bound. Based on our results, we show that, for the high signal-to-noise ratio regime, the performance of the least-squares estimator is significantly poorer than the Cramer-Rao bound, and we characterize this gap analytically. On the positive side, we show that for the challenging low signal-to-noise regime (attributed to either a weak astronomical signal or a noise-dominated condition) the least-squares estimator is near optimal, as its performance asymptotically approaches the Cramer-Rao bound. However, we also demonstrate that, in general, there is no unbiased estimator for the astrometric position that can precisely reach the Cramer-Rao bound. We validate our theoretical analysis through simulated digital-detector observations under typical observing conditions. We show that the nominal value for the mean-square-error of the least-squares estimator (obtained from our theorem) can be used as a benchmark indicator of the expected statistical performance of the least-squares method under a wide range of conditions. Our results are valid for an idealized linear (one-dimensional) array detector where intra-pixel response changes are neglected, and where flat-fielding is achieved with very high accuracy.

preprint2014arXiv

Analysis of the Cramer-Rao lower uncertainty bound in the joint estimation of astrometry and photometry

In this paper we use the Cramer-Rao lower uncertainty bound to estimate the maximum precision that could be achieved on the joint simultaneous (or 2D) estimation of photometry and astrometry of a point source measured by a linear CCD detector array. We develop exact expressions for the Fisher matrix elements required to compute the Cramer-Rao bound in the case of a source with a Gaussian light profile. From these expressions we predict the behavior of the Cramer-Rao astrometric and photometric precision as a function of the signal and the noise of the observations, and compare them to actual observations - finding a good correspondence between them. We show that the astrometric Cramer-Rao bound goes as $(S/N)^{-1}$ (similar to the photometric bound) but, additionally, we find that this bound is quite sensitive to the value of the background - suppressing the background can greatly enhance the astrometric accuracy. We present a systematic analysis of the elements of the Fisher matrix in the case when the detector adequately samples the source (oversampling regime), leading to closed-form analytical expressions for the Cramer-Rao bound. We show that, in this regime, the joint parametric determination of photometry and astrometry for the source become decoupled from each other, and furthermore, it is possible to write down expressions (approximate to first order in the small quantities F/B or B/F) for the expected minimum uncertainty in flux and position. These expressions are shown to be quite resilient to the oversampling condition, and become thus very valuable benchmark tools to estimate the approximate behavior of the maximum photometric and astrometric precision attainable under pre-specified observing conditions and detector properties.