Laplace近似积分

拉普拉斯方法又称为拉普拉斯近似（Laplace Approximation）。它可以用来计算一元或多元积分。
这里的思想类似于上篇博客中所讲的使用Laplace近似分布一样，这里把$\int_a^b f(x)dx$中的$f(x)$看成是一个分布(需要正规化)，然后使用一个正态分布近似这个分布，这样就可以得到原来积分的近似了。
拉普拉斯方法是有局限性的：被积分的函数有一个峰值，并且和正态曲线长的像。这种情况下的近似才能比较精确。
对于概率密度函数，（无论一维还是多维），大部分都是一个峰值。或者因为中心极限定理，统计量均值的分布和正态函数很像，也是一个峰值。因此拉普拉斯方法的用处还是很多的。此外，拉普拉斯方法计算很快。比如在线性混合效果模型中，拉普拉斯方法是用的最广的方法，也可能是唯一能在实际中使用的方法［2-4］。

［1］维基百科 http://en.wikipedia.org/wiki/Laplaces_method
［2］Approximate Inference in Generalized Linear Mixed Models. N. E. Breslow and D. G. Clayton. Journal of the American Statistical Association Vol. 88, No. 421 (Mar., 1993) , pp. 9-25
［3］Variance component testing in generalised linear models with random effects. XIHONG LIN. Biometrika (1997) 84 (2): 309-326.
［4］lme4 package Douglas Bates et al. http://cran.r-project.org/web/packages/lme4/lme4.pdf (see nAGQ parameter)