An L1-Analysis of the Iterative Proportional Fitting Procedure

  • Convergence of the Iterative Proportional Fitting procedure is analyzed. The input comprises a nonnegative weight matrix, and positive target marginals for rows and columns. The output sought is what is called the biproportional fit, a scaling of the input weight matrix through row and column divisors so as to equate row and column sums to target marginals. The procedure alternates between the fitting of rows, and the fitting of columns. We monitor progress with an L1-error function measuring the distance between current row and column sums and target row and column marginals. The procedure converges to the biproportional fit if and only if the L1-error tends zero. In case of non-convergence the procedure appears to oscillate between two accumulation points. The oscillation result is contingent on the conjecture that the even-step subsequence of the procedure is always convergent. The conjecture is established in the specific setting when the even-step subsequence has a connectedConvergence of the Iterative Proportional Fitting procedure is analyzed. The input comprises a nonnegative weight matrix, and positive target marginals for rows and columns. The output sought is what is called the biproportional fit, a scaling of the input weight matrix through row and column divisors so as to equate row and column sums to target marginals. The procedure alternates between the fitting of rows, and the fitting of columns. We monitor progress with an L1-error function measuring the distance between current row and column sums and target row and column marginals. The procedure converges to the biproportional fit if and only if the L1-error tends zero. In case of non-convergence the procedure appears to oscillate between two accumulation points. The oscillation result is contingent on the conjecture that the even-step subsequence of the procedure is always convergent. The conjecture is established in the specific setting when the even-step subsequence has a connected accumulation point, but remains open in general.show moreshow less

Download full text files

Export metadata

Statistics

Number of document requests

Additional Services

Share in Twitter Search Google Scholar
Metadaten
Author:Friedrich PukelsheimGND
URN:urn:nbn:de:bvb:384-opus4-12739
Frontdoor URLhttps://opus.bibliothek.uni-augsburg.de/opus4/1604
Series (Serial Number):Preprints des Instituts für Mathematik der Universität Augsburg (2012-02)
Type:Preprint
Language:English
Publishing Institution:Universität Augsburg
Release Date:2012/01/26
Tag:alternating scaling algorithm; biproportional fitting; entropy; matrix scaling; RAS procedure
GND-Keyword:Statistik; Kontingenztafelanalyse; L1-Approximation
Institutes:Mathematisch-Naturwissenschaftlich-Technische Fakultät
Mathematisch-Naturwissenschaftlich-Technische Fakultät / Institut für Mathematik
Mathematisch-Naturwissenschaftlich-Technische Fakultät / Institut für Mathematik / Lehrstuhl für Stochastik und ihre Anwendungen
Dewey Decimal Classification:5 Naturwissenschaften und Mathematik / 51 Mathematik / 510 Mathematik
Licence (German):Deutsches Urheberrecht mit Print on Demand