kth.sePublications
Change search
CiteExportLink to record
Permanent link

Direct link
Cite
Citation style
  • apa
  • ieee
  • modern-language-association-8th-edition
  • vancouver
  • Other style
More styles
Language
  • de-DE
  • en-GB
  • en-US
  • fi-FI
  • nn-NO
  • nn-NB
  • sv-SE
  • Other locale
More languages
Output format
  • html
  • text
  • asciidoc
  • rtf
Data Assimilation for Sign-indefinite Priors: A generalization of Sinkhorn's algorithm
KTH, School of Electrical Engineering and Computer Science (EECS), Intelligent systems, Decision and Control Systems (Automatic Control). KTH, School of Engineering Sciences (SCI), Mathematics (Dept.).
Mechanical and Aerospace Engineering, University of California, Irvine, CA 92697, USA.
Formerly with Computer Science and Applied Mathematics & Statistics, State University of New York, Stony Brook, NY 11794, USA.
2025 (English)In: Automatica, ISSN 0005-1098, E-ISSN 1873-2836, Vol. 177, article id 112283Article in journal (Refereed) Published
Abstract [en]

We deal with a hitherto unstudied problem that is motivated by the need to model promotion/inhibition contributions of relations, as in gene regulatory networks, to explain observed empirical distributions. There is no known pre-existing work focusing on this problem which, in essence, consists in calibrating a sign-indefinite data set representing prior information on the joint strength of promotion/inhibition, to restore consistency with specified observed marginals. We address the problem in the general setting of relations between n objects, and then specialize to binary relations to compare with a possible adhoc alternative following a Metropolis–Hastings rationale. Thus, our paper begins by addressing the general problem to find a multi-dimensional array (representing relative contribution of the n factors) to match specified marginals; the salient feature is that the sign (promotion/inhibition) of entries should match that of the prior. We develop an approach that is based on minimizing a suitable entropy functional subject to the constraints. Our approach follows verbatim Schrödinger's rationale (Schrödinger, 1931, 1932) that underlies the celebrated Sinkhorn–Knopp algorithm in statistics used to reconcile contingency tables (of positive/frequency data sets) with marginal data (Knight, 2008). The resulting algorithm for our case of sign-indefinite prior generalizes the Sinkhorn–Knopp algorithm in that it amounts to iterative scaling of the entries of the array along different coordinate directions. The scaling is multiplicative but also, in contrast to Sinkhorn–Knopp, inverse-multiplicative depending on the sign of entries. When the entries of the prior are non-negative, our algorithm reduces to the classical Sinkhorn–Knopp algorithm which, in recent years, has become the computational workhorse of optimal transport (Peyré and Cuturi, 2019). Interestingly, the extension of the Sinkhorn–Knopp algorithm herein is perhaps the only example of a coordinate gradient ascent algorithm, other than the classical Sinkhorn, that has a non-trivial closed-form solution.

Place, publisher, year, edition, pages
Elsevier BV , 2025. Vol. 177, article id 112283
Keywords [en]
Multi-marginal Schrödinger bridges, Negative probabilities, Sinkhorn algorithm
National Category
Computer Sciences
Identifiers
URN: urn:nbn:se:kth:diva-362508DOI: 10.1016/j.automatica.2025.112283Scopus ID: 2-s2.0-105001983143OAI: oai:DiVA.org:kth-362508DiVA, id: diva2:1952956
Note

QC 20250422

Available from: 2025-04-16 Created: 2025-04-16 Last updated: 2025-04-22Bibliographically approved

Open Access in DiVA

No full text in DiVA

Other links

Publisher's full textScopus

Authority records

Dong, Anqi

Search in DiVA

By author/editor
Dong, Anqi
By organisation
Decision and Control Systems (Automatic Control)Mathematics (Dept.)
In the same journal
Automatica
Computer Sciences

Search outside of DiVA

GoogleGoogle Scholar

doi
urn-nbn

Altmetric score

doi
urn-nbn
Total: 14 hits
CiteExportLink to record
Permanent link

Direct link
Cite
Citation style
  • apa
  • ieee
  • modern-language-association-8th-edition
  • vancouver
  • Other style
More styles
Language
  • de-DE
  • en-GB
  • en-US
  • fi-FI
  • nn-NO
  • nn-NB
  • sv-SE
  • Other locale
More languages
Output format
  • html
  • text
  • asciidoc
  • rtf