Information criteria and cross validation for Bayesian inference in regular and singular cases

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, Akaike’s 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.


Introduction
In data science, we estimate an unknown information source by a predictive distribution defined from a statistical model and a prior. In an older framework of Bayesian statistics in the 20th century, it was assumed that a statistical model is convinced to be correct and a prior is given by a subjective belief, resulting that the predictive distribution is believed to be subjectively optimal solution without any check or test. However, such a restricted treatment of Bayesian inference cannot be the foundation of today's data science in the real world, because almost all statistical models and learning machines are different from the data generating distributions, and their priors are chosen not by any a priori knowledge or subjective belief, but by the mathematical laws for the minimum generalization loss, the maximum marginal likelihood, or other optimization methods.
In statistics, a new paradigm was already proposed by Akaike (1980aAkaike ( , 1980b) that both a statistical model and a prior are understood as only candidate systems and they had better be optimized so that the predictive distributions are controlled to be the better approximations of the unknown probability distribution. Also it was proposed by Good (1952); ; Gelman and Shalizi (2013) that the marginal likelihood or the predictive distribution had better be checked and tested because Bayesian confirmation theory does not ensure the predictive distribution is suitable for the unknown distribution when a statistical model is too simple or too redundant. Nowadays, their proposal is widely accepted in statistics, data science, and machine learning, on which many statistical systems and learning machines are being applied to scientific and practical problems.  Wang et al. (2018).
In this paper, to develop the scientific evaluation measures of both a statistical model and a prior, we study the mathematical relation among the generalization loss, the information criteria AIC by Akaike (1974), DIC by Spiegelhalter et al. (2002), WAIC by Watanabe (2010), and the leave-one-out cross-validation loss by Gelfand et al. (1992); Vehtari and Lampinen (2002) in Bayesian inference from the following three points of view.
First, they are compared in the singular condition that the posterior distribution can not be approximated by any normal distribution. In many statistical models which have hierarchical structures or latent variables, the posterior distributions can not be approximated by any normal distribution, even asymptotically. In such cases, the generalization loss can not be estimated by information criteria AIC or DIC, however, it can be estimated by WAIC and the cross validation.
Second, they are studied in the case when a leverage sample point is contained in data. If a sample consists of independently and identically distributed random variables, then information criteria and the leave-one-out cross-validation loss are asymptotically equivalent, however, if otherwise, they are not. If a conditional probability is estimated using not an i.i.d. input sample, then there are cases when information criteria are better estimators of generalization loss than the cross-validation loss.
And last, their properties are clarified when they are used for the prior optimization problem. Information criteria and cross validation loss can be understood as the evaluation measures of the prior design. Minimization of WAIC and leave-oneout cross validation loss makes the average generalization loss, however, it does not minimize the generalization loss as a random variable.
This paper consists of five sections. In the second section, we define the generalization loss, the information criteria, and the cross-validation loss. In the third section, they are compared from the three different points of view. In the fourth section, their statistical difference is discussed, and in the last section, we conclude the paper.

Definitions of information criteria
In this section, we define the generalization loss, information criteria, and cross validation.

Definitions of statistical inference
Let q(x) be a probability density function of x ∈ ℝ N and X 1 , X 2 , ..., X n be an i.i.d. sample, in other words, a sample is a set of independent random variables whose probability density function is q(x). In this paper, we mainly study i.i.d. cases. For the case when a sample is not i.i.d., see the Sect. 3.2. The cross validation procedure needs the i.i.d. condition, whereas information criteria can be used in several not i.i.d. cases as shown in Watanabe (2021).
A statistical inference is defined by a map where p * (x) is an estimated probability density function on ℝ N , which is often called a predictive distribution. There are many procedures in statistical inference, for example, the maximum likelihood method, the maximum a posteriori (MAP) method, the Bayesian method, the sparse estimation method, and so on. In this paper, we study the Bayesian method. Let p(x|w) and (w) be a statistical model and a prior, respectively, where w ∈ ℝ d is a parameter. In the maximum likelihood method, the predictive distribution is defined by p * (x) = p(x|w * ) , where w * is the maximum likelihood estimator. In the MAP method, the maximum likelihood estimator is replaced by the maximum a posteriori estimator. In the Bayesian method, the posterior distribution is defined by where Z n is a normalizing constant referred to as the marginal likelihood. For an arbitrary function of f(w), let w [f (w)] and w [f (w)] denote the posterior average and 1 3 the posterior variance of f(w), respectively. The Bayesian predictive distribution is defined by Note that, in data science, we do not know whether a statistical model and a prior are appropriate or not. For a given predictive distribution p * (x) , its generalization loss G n is defined by Since X n is a random variable, G n is also a random variable. Then minimizing G n is equivalent to minimizing Kullback-Leibler divergence K(q||p * ) , where Hence the generalization loss G n gives the quantitative evaluation of a pair of a model and a prior according to the Kullback-Leibler divergence, however, since q(x) is unknown in general, we need mathematical theory to estimate G n . The training loss T n is defined by The generalization loss G n can not be estimated by T n , because the expectation values of G n and T n are not equal to each other. If a more complicated model is employed, then the training loss decreases, but the generalization loss may increase.
Remark It is shown by the Bayesian decision theory that, if w is generated from some true Φ(w) and if X n is generated from some true ∏ n i=1 P(X i �w) , then the expectation value of the generalization loss is made smallest if and only if (w) = Φ(w) and p(x|w) = P(x|w) . However, in almost all situations of statistical inference, both the true prior and the true model are unknown, hence the decision theory can not be the base of Bayesian inference of unknown information source q(x) in a real world. If q(x) is located outside of a set of statistical model, then min max theory can not be used for designing the prior.

Asymptotic generalization loss
In this subsection, we explain the asymptotic behavior of the generalization loss of Bayesian inference. We need several mathematical definitions. Let W 0 ⊂ ℝ d be the set of all parameters that minimize If there exists w 0 ∈ W 0 such that q(x) = p(x|w 0 ) , then an unknown distribution q(x) is said to be realizable by a statistical model p(x|w). If W 0 consists of a single element w 0 and the Hessian matrix is positive definite, then q(x) is said to be regular for p(x|w). Many statistically asymptotic theories assume the regularity condition, however, if a statistical model has a hierarchical structure or latent variables, then regularity condition is not satisfied. In fact, neural networks, normal mixtures, hidden Markov models, matrix factorizations, latent Dirichlet allocations, and so on do not satisfy the regularity condition, resulting that their posterior distributions are far from any normal distribution. Even in such singular cases, the generalization performance of Bayesian inference is clarified based on algebraic geometry by Watanabe (2009). Let L n (w) be the empirical log loss function, We assume that, for an arbitrary parameter w 0 ∈ W 0 , p(x|w 0 ) represents the same probability density function. Note that, if the unknown distribution is realizable by a statistical model, then L(w 0 ) and L n (w 0 ) are equal to the entropy S and the empirical entropy S n of q(x), respectively. Even if an unknown distribution q(x) is neither realizable by nor regular for a statistical model, it is proved by singular learning theory in Watanabe (2009) that there exists a random variable n such that which satisfies the convergence in distribution n → . We need the zeta function of Bayesian statistics for z ∈ ℂ by It is proved based on the resolution theorem in Hironaka (1964) that (z) is a holomorphic function in ℜ(z) > 0 which can be analytically continued to the unique meromorphic function onto the entire complex plane, whose poles are all real and negative values. Let (− (w 0 )) be the largest pole of the zeta function. Then, by using algebraic geometry, it is proved that Therefore,

3
This equation shows that the generalization loss is the sum of the bias L(w 0 ) and the variance (w 0 )∕n . If an unknown distribution is regular for a statistical model, then (w 0 ) = d∕2 , where d is the dimension of the parameter space, if otherwise, (w 0 ) ≤ d∕2 with an assumption (w 0 ) > 0 for some w 0 ∈ W 0 . The constant (w 0 ) is referred to as the real log canonical threshold, which is a birational invariant, in other words, the generalization loss of Bayes inference is determined by algebrogeometric structure of a statistical model.
Concrete values of them were found in the matrix factorization by Aoyagi and Watanabe (2005), the Poisson mixture by Sato and Watanabe (2019), the latent Dirichlet allocation by Hayashi (2021), and many statistical models and learning machines such as Yamazaki and Watanabe (2003); Yamazaki (2016); Yamazaki and Kaji (2013); Zwiernik (2011);Watanabe (2009). They are useful in the design of Markov chain Monte Calro by Nagata and Watanabe (2008) and the singular Bayesian information criterion by Drton and Plummer (2017).
Remark From the mathematical point of view, the real log canonical threshold (w 0 ) can be understood as a volume dimension of the neighborhood of a set {w ; L(w) − L(w 0 ) = 0} . In fact, it is proved in Watanabe (2009) that The prediction accuracy of Bayes inference is determined by the volume of the set of the almost optimal parameters.

Information criteria and cross validation
In this subsection, we introduce several information criteria to estimate the generalization loss. For the other information criteria to estimate the marginal likelihood, see Sect. 4. The asymptotic behavior of the generalization loss was clarified by Eq.(3), however, the real log canonical threshold depends on the unknown distribution, resulting that w 0 is also unknown, hence it cannot be directly applied to practical problems.

Remark
In this paper, we study the scale of information criteria and cross validation as estimators of the generalization loss [G n ] . In many practical problems, it should be remarked that their scales are defined so as to estimate 2n × [G n ] , according to the Akaike's pioneer work.
The concept of the information criterion was firstly proposed by Akaike (1974). The definition of AIC is given by which is called Akaike information criterion (AIC). At first, AIC was made for estimating the generalization loss of the maximum likelihood method, but it can be employed in Bayes, because if the unknown distribution is realizable by and regular for a statistical model, then the order of the difference between the above AIC and AIC using the maximum likelihood estimator is smaller than 1/n. For Bayesian estimation, the deviance information criterion (DIC) was proposed by Spiegelhalter et al. (2002), by which the Bayesian generalization loss can be estimated, if an unknown distribution is realizable by and regular for a statistical model. The widely applicable information criterion (WAIC) was proposed by The value w [log p(X i |w)] shows the fluctuation of log p(X i |w) according to the posterior distribution. The mathematical relation between the generalization loss and this value is shown by using the functionally partial integration over the limit Gaussian process of the empirical process in Watanabe (2009). This relation has the same mathematical structure as the fluctuation-dissipation theorem which shows the correspondence from the system is equal to the fluctuation of the system. The leave-one-out cross-validation loss in Gelfand et al. (1992); Vehtari and Lampinen (2002); Vehtari et al. (2017) is defined by where X n ⧵ X i is the sample which does not contain X i . If the posterior distribution is precisely realized, then LOO can be numerically approximated by using the importance sampling cross-validation loss, Even if an unknown distribution is unrealizable by or singular for a statistical model, the generalization loss can be estimated by WAIC, LOO, and ISCV. Note that AIC, DIC, and WAIC are calculated numerically using Markov chain Monte Carlo (MCMC) method, if the posterior variance V(i) = w [log p(X i |w)] is finite, The value V(i) shows the leverage of a sample point X i which will be discussed in Sect. 3.2. Note that LOO needs n different posterior distributions for X n ⧵ X i , whereas ISCV can be calculated by one posterior distribution for X n . However, in ISCV, the average w [1∕p(X i |w)] or variance w [1∕p(X i |w)] may be larger or infinite if a leverage sample point is contained, which was proved by Peruggia (1997); Epifani et al. (2008), resulting that the posterior calculation by MCMC method fails, 1 3 and a numerical improvement was developed by Gelman et al. (2014). The difference of information criteria and the leave-one-out cross validation in influential observation cases is explained in Sect. 3.2.

Comparison of information criteria and cross validation
In this section, we compare information criteria and cross validation loss in three different problems, that is to say, (1) in singular posterior cases, (2) in influential observation cases, and (3) in hyperparameter optimization cases.

Regular and singular cases
Firstly, we compare information criteria and leave-one-out cross-validation loss in regular and singular cases.

Regular Cases.
If an unknown distribution q(x) is realizable by and regular for a statistical model p(x|w), then the posterior distribution can be approximated by a normal distribution when the sample size n tends to infinity. On such a regularity condition, AIC, DIC, WAIC, LOO, and ISCV are asymptotically equivalent to each other as random variables. That is to say, They are all asymptotically unbiased estimators of the generalization loss, However, the generalization loss and AIC have asymptotically inverse correlation as is shown by Watanabe (2018), and DIC, WAIC, LOO, and ISCV satisfy the same equation as Eq.(5). This property strongly affects the model selection procedures by AIC, DIC, WAIC, LOO, and ISCV, which is the weak point of both the information criteria and the cross-validation loss. It should be emphasized that, since Eq.(5) is not trivial, many users of the information criteria and the cross validation loss may not be aware of such a risk.

Singular Cases.
In general, the posterior distribution can not be approximated by any normal distribution even if the sample size tends to infinity in . (5) Even if the posterior distribution is far from any normal distribution, it was proved in Watanabe (2010) that and that In such singular cases, neither AIC nor DIC satisfies these equations. In singular cases, the posterior expectation value of the parameter w [w] in DIC is not appropriate for statistical inference. The generalization loss and WAIC have asymptotically inverse correlation as is shown by Watanabe (2010), and both LOO and ISCV satisfy the same equation as Eq. (7). In singular cases, the real log canonical threshold (w 0 ) is not larger than the regular cases. Hence, even if a statistical model p(x|w) is redundant for an information source q(x), the increase of the Bayesian generalization loss is smaller than the regular cases, which is one of the good properties of Bayesian inference. In the model selection problem, we should remark the fact that the increases of WAIC, LOO, and ISCV are also smaller than the regular cases.
Example 1 (Normal mixture in regular and singular cases). Let us compare AIC, DIC, WAIC, and ISCV in regular and singular cases. A normal distribution of x ∈ ℝ N which has a mean b k and a variance and 1∕s k is denoted by We study a statistical model and a prior given in the following equations, where w = (a, b, s) is a parameter, a = {a k } , b = {b k } , and s = {s k }, Here K is the number of components in the statistical model. Note that the mixture ratio a = {a k } satisfies a k ≥ 0 and ∑ k a k = 1 , and ( , r, , ) is a hyperparameter. A normal mixture is a typical singular model, because if a model is redundant, then W 0 is a set with singularities as pointed out by Yamazaki and Watanabe (2003). Let us show simple experimental results for regular and singular cases as is shown in Watanabe (2021). An experiment was set as N = 3 , n = 100 , r = N∕2 , = 1 , and = = 1 . A distribution q(x) was set as p(x|w 0 ), . For a regular case, K 0 = 4 was used, whereas, for a singular case, K 0 = 2 was used. These two distributions were estimated by a normal mixture of K = 4 components. The posterior distributions were approximated by using the Gibbs sampler of the simultaneous distribution of the parameter and latent variable. In Table 1, means and standard deviations of G n − S , AIC − S n , DIC − S n , WAIC − S n and ISCV − S n are shown which were calculated by using 100 independent trials, where S = L(w 0 ) and S n = L n (w 0 ) are entropy and empirical entropy of p(x|w 0 ) . In the regular case, all information criteria and ISCV could approximate the generalization loss, whereas, in the singular case, averages of AIC and DIC were larger and smaller than the generalization loss, respectively, whereas averages of WAIC and ISCV were almost equal to the generalization loss. The sample size n = 100 is not sufficiently large according to the asymptotic theory, however, information criteria and cross validation loss were almost equal to each other.

Influential observation
It is often assumed that a sample is independently and identically distributed (i.i.d.) in theoretical studies of information criteria. Such an i.i.d. condition is necessary for the use of the cross validation loss, however, information criteria can be employed in several not i.i.d. cases, for example, conditional independent cases. In this subsection, we study the statistical inference of a conditional probability density for dependent inputs. Assume that x n = {x i ;i = 1, 2, … , n} is not a set of random variables but a constant sequence and that {Y i ;i = 1, 2, … , n} is a set of an Table 1 Means and standard deviations of information criteria in regular and singular cases A normal mixture of four components was used as a statistical model, and two different distributions were estimated. In a singular case, AIC and DIC were larger and smaller than the generalization loss, respectively. WAIC and ISCV could estimate the generalization loss even if the posterior distribution was singular independent random variables whose probability distribution is ∏ n i=1 q(y i �x i ) . In this situation, the posterior average w [ ] and the variance w [ ] are defined by the posterior distribution The generalization and training losses are defined by The information criteria and leave-one-out cross-validation loss are defined by the same way as foregoing equations, Note that, if x n is independent and identically distributed (i.i.d.), then [ISCV] = [G n−1 ] by definition, whereas, if otherwise, [ISCV] ≠ [G n−1 ] . For example, if x n is fixed or dependent, ISCV is not an unbiased estimator of the generalization loss in general. On the other hand, even for a fixed x n , information criteria are asymptotically unbiased estimators of the generalization loss if the sample size is sufficiently large. If an unknown distribution is realizable by and regular for a statistical model, then AIC, DIC, and WAIC can be employed, which is shown by .
Example 2 (Influential observation in linear regression) We study a simple regression model of x, y ∈ ℝ and a prior (w), (12) p(y|x, a, s) =(s∕(2 )) 1∕2 exp(−(s∕2)(y − ax) 2 ), 1 3 where a ∈ ℝ, s > 0 are parameters and = 0.01 is a hyperparamater. Let the sample size n = 10 , and we set the x n as x i = i∕10 (i = 1, 2, … , 9) and We compare information criteria and importance sampling leave-one-out cross-validation loss for two cases x 10 = 4 , which is a leverage sample point, and x 10 = 1 , which is not. A regression problem in which a leverage sample point is contained is often called an influential observation problem. The output Y n are independently taken from ∏ i p(y i �x i , a 0 , s 0 ) , where a 0 = 0.1 , 1∕s 0 = 0.01 . Then S and S n are conditional entropy and empirical one, respectively. The average and standard deviations of the generalization loss, information criteria, and importance sampling leave-oneout cross-validation loss are shown in Table.2. In an influential observation case, the posterior average and variance of w [1∕p(Y i |X i , w)] are larger or infinite by Peruggia (1997); Epifani et al. (2008), and both the average and variance of the [ISCV] were larger than that of the generalization loss, where the averages of information criteria were asymptotically equal to that of generalization loss. In general, whether a sample point is leverage one or not depends on not only a sample point itself but also a statistical model. The author would like to recommend that both information criteria and cross validation had better be calculated and compared. If they are different, a leverage sample point is contained in the data, which can be found because it makes the functional variance larger than the other sample points, which is shown in .

Table 2 Means and standard deviations of information criteria for fixed inputs
The case x 10 = 4 is a leverage sample point, which is the influential observation, but x 10 = 1 is not. Cross-validation loss can not be used if an input is a leverage sample point In this example, we study the difference of two generalization losses, The first one G (1) n is the generalization loss for the case {(X i , Y i )} is the set of independent random variables subject to q(x)q(y|x), whereas the second one G (2) n is that for the case {(Y i |x i )} is the set of conditionally independent random variables subject to q(y i |x i ) for given n . Two entropies are defined by Information criteria AIC, DIC, WAIC, and ISCV are defined by the same way as Eqs. (8)

-(11). LOO is defined by
In the experiment, {x i } was independently taken from and fixed. {Y i } was independently taken from p(y|x i , a 0 ) , where a 0 = (1, 1, ..., 1) ∈ ℝ M . The sample size n = 50 was fixed, and the dimension M was set as M = 10, 20, 50, 100, 200 . In high dimensional cases, we can not assume that n is sufficiently large and there are many leverage sample points. Two different generalization losses were calculated. In Table 3, GE1 and GE2 denote G (1) n − S (1) and G (2) n − S (2) , respectively. LOO, ISCV, WAIC, DIC, and AIC show the values of LOO-S n , ISCV-S n , WAIC-S n , DIC-S n , and AIC-S n , respectively. Note that the empirical entropy does not depned on the assumption of inputs. In the experiment, LOO and WAIC estimated averages of GE1 and GE2, respectively. The standard deviations of LOO were larger than other criteria. In a high dimensional case, ISCV was not equal to LOO, because the posterior average in ISCV by MCMC was not calculated precisely. a)]dxdy,

Prior optimization problem
In this subsection, we study a prior optimization problem in regular cases. That is to say, for a fixed statistical model, the generalization loss, information criteria, and cross-validation loss are understood as functionals of priors, and their minimization problem is studied. In this section, a candidate prior (w) is assumed to be a general nonnegative function of a parameter w ∈ ℝ d which may be improper, that is to say ∫ (w)dw may be infinite. Even in such cases, the posterior and predictive distributions can be defined by the same equations as Eqs. (1), (2). Let G( ) , WAIC( ) , and LOO( ) be the generalization loss, WAIC, and the leave-one-out cross-validation loss which are functionals of a candidate prior (w).

Remark
In this subsection, we assume the regularity condition that an unknown distribution is regular for a statistical model, however, it may be unrealizable by a statistical model. The general theory for singular cases is not yet constructed. If a posterior distribution is singular, then it has the phase transition according to hyperparameter control as is explained in , in other words, the asymptotic support of the posterior distribution drastically changes at some critical value of the hyperparameter.
Let 0 (w) be an arbitrary fixed prior. For example, one can choose 0 (w) = 1 for an arbitrary w. For a given candidate prior (w) , we define (w) = (w)∕ 0 (w). If 0 (w) ≡ 1 is chosen, then (w) = (w) . Assume that the posterior distribution can be approximated by a normal distribution. Then it was proved in  that there exists a function M( , w) of (w) and w which satisfies Also it was proved in  that there exists a function M( , w) of (w) and w which satisfies The concrete forms of M( , w) and M( , w) are defined by using higher order differentials of the log density function log p(x|w) , which are given in . They satisfy It should be emphasized that the generalization loss as a random variable has different behavior from the information criteria and leave-one-out cross-validation loss.
Hence minimization of WAIC or LOO makes the average generalization loss asymptotically, however, it does not minimize the generalization loss as a random variable. Note that minimization of AIC or DIC does not minimize the generalization loss as either the average or the random variable.
Example 4 (Prior Optimization in linear regression) A simple regression model of x, y ∈ ℝ same as Eqs. (12) and (13) in Example.2 is studied. In this example, x n is independently taken from the standard normal distribution, and Y n is also independently taken from p(y i |x i , a 0 , s 0 ) , where n = 10 , a 0 = 0.1 , and 1∕s 0 = 0.01 . The generalization loss, AIC, DIC, WAIC, and importance samping leave-one-out crossvalidation loss of 50 candidate hyperparameters 0 < ≤ 0.2 were compared, and the hyperparameters which made them smallest were obtained for 100 independent samples. They are shown in Table 4. Note that which were different. By using WAIC or ISCV, the former was estimated, whereas it was not estimated by AIC or DIC. Based on the higher order equivalence, WAIC and ISCV can be used for prior optimization for the average generalization loss. In the experiment, the variance by WAIC was smaller than that of ISCV.

Discussion: WBIC as an extension of BIC
In the previous sections, we studied the information criteria and leave-one-out crossvalidation loss as estimators of the generalization loss. In this discussion, we compare the AIC type criteria with BIC ones.
Before discussion, we summarize the statistical property of the marginal likelihood. The marginal likelihood Z n is defined by Eq.(1), which can be understood as an estimated probability density of X n by a statistical model and a prior, because it satisfies ∫ Z n dx n = 1 . The free energy F n , which is equal to the minus log marginal likelihood, is also defined by Then its expectation value is equal to where S is the entropy of q(x), q(x n ) = ∏ n i=1 q(x i ) and p(x n ) = Z n . Hence the expectation value of the free energy is smallest if and only if the Kullback-Leibler divergence between q(x n ) and p(x n ) is smallest. Moreover, it is immediately derived that Its asymptotic value is given by where (w 0 ) is the real log canonical threshold. If the posterior distribution can be approximated by a normal distribution, then (w 0 ) = d∕2 and F n is asymptotically approximated by BIC in Schwarz (1978). Even in singular cases, F n can be estimated by sBIC in Drton and Plummer (2017) and WBIC in Watanabe (2013).
Let us compare the information criteria and the marginal likelihood from two different points of view. For short description, the information criteria AIC, DIC, WAIC, and LOO are referred to as AICs, whereas BIC, sBIC, WBIC, and the free energy are BICs.
First, we compare AICs and BICs from the viewpoint of model selection consistency. Let P = {p k (x|w k ) ; k = 1, 2, … , K} be a set of candidate statistical models. Assume that there exists a statistical model in P by which q(x) is realizable. The set of models P 0 is defined by the set of all statistical models in P by which q(x) is realizable. For a given q(x), the statistical model p ∈ P 0 that has the smallest real log canonical threshold is called the smallest true model. A model selection criterion is called to have the model selection consistency if and only if the probability that the smallest true model is chosen converges to one as the sample size tends to infinity. The asymptotic properties of AICs and BICs are given by If q(x) is unrealizable by p(x|w), then the convergence in probability holds, If q(x) is realizable by p(x|w), then in both regular and singular cases, n(T n − S n ) converges to a random variable in distribution in Watanabe (2009. In AICs, the penalty term has the same order as the standard deviation of T n − S n , whereas in BICs, it is larger than that of T n − S n . It follows that AICs do not have model selection consistency but BICs have, in both regular and singular cases. Second, we study AICs and BICs from the viewpoint of asymptotic efficiency. Assume that q(x) is not realizable by any model in P . Then the optimal model for efficiency is defined by the model that minimizes [G n ] , which depends on the sample size n. Remark that consistency is defined for the case n → ∞ whereas efficiency is for finite n. A model selection criterion is said to have more efficiency if and only if [G n ] is made smaller. If BICs are employed, then the penalty term is too large for the smaller generalization loss, hence AICs have more efficiency than BICs, in both regular and singular cases. In practical applications of data science, appropriate information criteria had better be employed according to the purpose of the data analysis. T n − S n → min w KL(q(x)||p(x|w)) > 0.

Conclusion
Information criteria and leave-one-out cross-validation loss in Bayesian modeling are studied, and three points are clarified. First, if the posterior distribution is far from any normal distribution, then AIC and DIC cannot be used for model evaluation, whereas WAIC and leave-one-out cross-validation loss can be employed. Second, if inputs in a sample are not independent, then leave-one-out cross-validation loss is not an asymptotic unbiased estimator of the generalization loss. And last, in the hyperparameter optimization problem, WAIC and leave-one-out cross-validation loss can be used for minimum average generalization loss.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.