(Postal Mail)
Sumio Watanabe, Ph.D. Professor of Department of Mathematical and Computing Science Tokyo Institute of Technology, Mail-Box W8-42, 2-12-1, Oookayama, Meguro-ku, Tokyo, 152-8552, Japan (E-mail) swatanab (AT) c . titech . ac . jp Japanese Homepage Google Scholar DBLP: Computer Science Bibliography Paper Information Return to Watanabe Lab. |
In 1998, we found a bridge between algebraic geometry and learning theory. |

Please search WAIC in studies of covid-19.

WAIC in Practical Problems

Seasonal Institute of Mathematics : Deepening and Evolution of Applied Singularity Theory, 2022, Nov, 20-25,Yokohama Japan.

Statistical problems caused by singularities can be mathematically resolved based on algebraic geometry.

Sumio Watanabe, Mathematical Theory of Bayesian Statistics for Unknown Information Source

In statistical inference, uncertainty is unknown and all models are wrong. A person who makes a statistical model and a prior distribution is simultaneously aware that they are fictional and virtual candidates. In order to study such cases, several statistical measures have been constructed, such as cross validation, information criteria, and marginal likelihood, however, their mathematical properties have not yet been completely clarified when statistical models are under- and over- parametrized. In this paper, we introduce a place of mathematical theory of Bayesian statistics for unknown uncertainty, on which we show general properties of cross validation, information criteria, and marginal likelihood. The derived theory holds even if an unknown uncertainty is unrealizable by a statistical model or even if the posterior distribution cannot be approximated by any normal distribution, hence it gives a helpful standpoint for a person who cannot believe in any specific model and prior. The results are followings. (1) There exists a more precise statistical measure of the generalization loss than leave-one-out cross validation and information criterion based on the mathematical properties of them. (2) There exists a more efficient approximation method of the free energy, which is the minus log marginal likelihood, even if the posterior distribution cannot be approximated by any normal distribution. (3) And the prior distributions optimized by the cross validation and the widely applicable information criterion are asymptotically equivalent to each other, which are different from that by the marginal likelihood.

https://arxiv.org/abs/2206.05630 , github

Sumio Watanabe, Mathematical theory of Bayesian statistics where all models are wrong. Advancements in Bayesian Methods and Implementations, Handbook of statistics, vol.47, Elsevier, 2022, September.

Nowadays, we know all models are wrong, both subjectively and objectively. We are aware that, if a model is wrong, coherent inference takes us to a wrong prediction and wrong science. Statisticians should construct new theory and method which can be employed even if a person is aware that all models are wrong.

In this chaper of the book, we show that there are mathematical theorems in Bayesian statistics, which hold even if statistical models are wrong or overparametrized. This article introduces such mathematical theorems, which enable us to find useful models and prior distributions in unknown large world.

Lecture Notes on Statstical Leaning Theory

Singular learning theory is introduced for students in mathematics and computer science. Mathematical learning theory was establised based on algebraic geometry, which is truly useful in unknown real world.

Overparametrized statistical models and learning machines can be studied by algebraic geometry. The real world is connected to algebraic geometry via zeta function.

Sumio Watanabe (2021) Information criteria and cross validation for Bayesian inference in regular and singular cases.

Japanese Journal of Statistics and Data Science.

https://doi.org/10.1007/s42081-021-00121-3

In data science, an unknown information source is estimated by a predictive distribution defined from a statistical model and a prior. In an older Bayesian framework, it was explained that the Bayesian predictive distribution should be the best on the assumption that a statistical model is convinced to be correct and a prior is given by a subjective belief in a small world. However, such a restricted treatment of Bayesian inference cannot be applied to highly complicated statistical models and learning machines in a large world. In 1980, a new scientific paradigm of Bayesian inference was proposed by Akaike, in which both a model and a prior are candidate systems and they had better be designed by mathematical procedures so that the predictive distribution is the better approximation of unknown information source. Nowadays, Akaikefs proposal is widely accepted in statistics, data science, and machine learning. In this paper, in order to establish a mathematical foundation for developing a measure of a statistical model and a prior, we show the relation among the generalization loss, the information criteria, and the cross-validation loss, then compare them from three different points of view. First, their performances are compared in singular problems where the posterior distribution is far from any normal distribution. Second, they are studied in the case when a leverage sample point is contained in data. And last, their stochastic properties are clarified when they are used for the prior optimization problem. The mathematical and experimental comparison shows the equivalence and the difference among them, which we expect useful in practical applications.

(Note) The difference between cross validation and information criteria is shown both theoretically and experimentally. Cross validation needs simultaneous independency of {(Xi,Yi)}, whereas information criteria does only conditional independency of {(Yi|Xi)}. That is to say, information criteria can be used even if {Xi} is not independent.

For example, cross validation cannot be used in the autoregressive models, whereas WAIC can be.

Sumio Watanabe (2021) WAIC and WBIC for mixture models. Bihaviormetrika,

https://doi.org/10.1007/s41237-021-00133-z

In this paper, we introduce mathematical foundation and computing methods of WAIC and WBIC in a normal mixture which is a typically singular statistical model, and discuss their properties in statistical inference. Also, we study the case that samples are not independently and identically distributed, for example, they are conditional independent or exchangeable. Furthermore, a simple calculation method of WBIC in mixture models is proposed.

Note: In a Gibbs sampler, it was difficult to calculate WBIC. However, we made a simple calculation method of WBIC in the above paper. The following is a Matlab file.

Matlab file : WAIC and WBIC for a normal mixture

Mathematical Theory of Bayesian Statistics

This book may be useful those who want the scientific Bayesian framework.

Nowadays, we know that all models are wrong, both subjectively and objectively.

In statistics and machine learning, a data-generating distributuon (DGP) is estimated by a predictive distribution defined from a statistical model and a prior. In an older Bayesian framework in the 20th century, it was explained that the Bayesian predictive distribution should be the best decision, on the assumption that a statistical model is convinced to be correct and a prior is given by a subjective belief in a small world. However, such a restricted formalism of Bayesian inference cannot be applied to highly complicated statistical models and learning machines in a large world.

In this book, in order to establish the mathematical foundation for Bayesian stistsics in the large world, it is shown that mathematical theorems universally hold for an arbitrary triple (a data-generating distribution, a statistical model, a prior).

This book may be useful for the readers who are interested in the followings.

(1) In the real world, all models are wrong. If a model is wrong, coherent inference takes us to a wrong prediction.

(2) If we assume that 'uncertainty' is represented by a model and 'subjective belief' is represented by a prior in the older Bayesian framework, we automatically accept the wrong assumption that the model and the prior are DGP.

(3) If we think that our own model and prior are DGP, our decision is not scientific. Because we know that our model is not DGP in the real world.

(4) Thus, we need to check or test both a model and a prior, even in Bayesian data analysis. A pair of a model and a prior is only a candiate system which is not DGP. We should improve both a model and a prior.

(4.1). The older Bayesian formalism prohibited to improve a statistical model and a prior, which is not appropriate in science.

(4.2). The old Bayesian formalism said that a probability distribution can be assigned for even non-probabilistic phenomenon by the belief of subjective decision. In science, such a claim may take us to misconducted researches. In science, we should verify both a model and a prior based on experiments.

(5) There are new mathematical theorems in Bayesian statistics. They hold even if a model is too small or too redundant. We had better know them in the unknown real world. All models are wrong, but we may find some useful models, based on mathematical theory. They depend on the sample size in general.

(6) Even if the posterior distribution cannot be approximated by any normal distribution, its generalization performance is determined by the birational invariant, the real log canonical threshold.

You can learn the new mathematical theory in Bayesian statistics: Statistical Leaning Theory

You can find phase transition in neural network learning. Sumio Watanabe, Algebraic geometrical methods for hierarchical learning machines. Neural Networks, 14(8), 2001, pp.1049-1060

Neural networks and normal mixtures are singular. You can find the phase transition phenomenon in Bayesian posterior distribution, as the sample size increases. For details, see Phase transition .

The reason why singular learning theory is necessary in deep learning.

Singular learning theory can be applied to neural networks, because it holds even if Fisher information matrix is highly singular. The generalization loss and log marginal likelihood are given by the real log canonical threshold, which is the volume dimension of singularities. WAIC and WBIC can be used in neural networks. AIC is far larger than the Generalization loss.

WAIC, LOO, AIC, and Generalization Loss in Neural Network (mp4) .

LOO requires independency of (Xi,Yi), whereas information criteria do only conditional independency of (Yi|Xi). Thus information criteria can be used more widely than LOO, for example AIC and WAIC can be applied to AR model in time series analysis.

(See: S. Watanabe, Mathematical theory of Bayesian statistics, 2018, CRC Press).

Sumio Watanabe, homepage continued