• Communications for Statistical Applications and Methods
• /
• v.19 no.4
• /
• pp.597-605
• /
• 2012

Evidence framework enables us to determine the tuning constant in a penalized likelihood formula. We apply the framework to the estimating parameters of normal mixtures. Evidence, which is a solely data-dependent measure, can be evaluated by Laplace approximation. According to a synthetic data simulation, we found that the proper values of the tuning constant can be systematically obtained.

• Communications for Statistical Applications and Methods
• /
• v.16 no.6
• /
• pp.971-978
• /
• 2009

Poisson generalized linear mixed models(GLMMs) have been widely used for the analysis of clustered or correlated count data. For the inference marginal likelihood, which is obtained by integrating out random effects is often used. It gives maximum likelihood(ML) estimator, but the integration is usually intractable. In this paper, we propose how to obtain the ML estimator via Laplace approximation based on hierarchical-likelihood (h-likelihood) approach under the Poisson GLMMs. In particular, the h-likelihood avoids the integration itself and gives a statistically efficient procedure for various random-effect models including GLMMs. The proposed method is illustrated using two practical examples and simulation studies.

• Journal of the Korean Data and Information Science Society
• /
• v.26 no.2
• /
• pp.495-504
• /
• 2015

The Gaussian geostatistical model has been widely used for modeling spatial data. However, this model suffers from a severe difficulty in computation because inference requires to invert a large covariance matrix in evaluating log-likelihood. In addressing this computational challenge, three strategies have been employed: likelihood approximation, lower dimensional space approximation, and Markov random field approximation. In this paper, we reviewed statistical approaches attacking the computational challenge. As an illustration, we also applied integrated nested Laplace approximation (INLA) technology, one of Markov approximation approach, to real data to provide an example of its use in practice dealing with large spatial data.

• Communications for Statistical Applications and Methods
• /
• v.19 no.3
• /
• pp.451-457
• /
• 2012

Baumgartner et al. (1998) proposed a novel statistical test for the null hypothesis that two independently drawn samples of data originate from the same population, and Murakami (2006) generalized the test statistic for more than two samples. Whereas the expressions of the exact density and distribution functions of the generalized Baumgartner statistic are not yet found, the characteristic function of its limiting distribution has been obtained. Due to the development of computational power, the Fourier series approximation can be readily utilized to accurately and efficiently approximate its density function based on its Laplace transform. Numerical examples show that the Fourier series method provides an accurate approximation for statistical quantities of the generalized Baumgartner statistic.

• IEIE Transactions on Smart Processing and Computing
• /
• v.4 no.4
• /
• pp.202-208
• /
• 2015

In this study, we propose a new inference algorithm for a multiclass Gaussian process classification model using a variational EM framework and the Laplace approximation (LA) technique. This is performed in two steps, called expectation and maximization. First, in the expectation step (E-step), using Bayes' theorem and the LA technique, we derive the approximate posterior distribution of the latent function, indicating the possibility that each observation belongs to a certain class in the Gaussian process classification model. In the maximization step, we compute the maximum likelihood estimators for hyper-parameters of a covariance matrix necessary to define the prior distribution of the latent function by using the posterior distribution derived in the E-step. These steps iteratively repeat until a convergence condition is satisfied. Moreover, we conducted the experiments by using synthetic data and Iris data in order to verify the performance of the proposed algorithm. Experimental results reveal that the proposed algorithm shows good performance on these datasets.

• Communications for Statistical Applications and Methods
• /
• v.2 no.1
• /
• pp.137-144
• /
• 1995

Hierarchical Bayes estimations of exchangeable mean vector of a multivariate Poisson distribution are obtained. Since sophiscated analytic integration procedures are needed, the Laplace method is employed in order tocompute these estimations approximately. An example is presented.

• Communications for Statistical Applications and Methods
• /
• v.14 no.2
• /
• pp.389-399
• /
• 2007

In this paper, Bayes estimates for the parameters k, c and reliability function of the Burr type XII model based on a type II censored samples under asymmetric loss functions viz., LINEX and SQUAREX loss functions are obtained. An approximation based on the Laplace approximation method (Tierney and Kadane, 1986) is used for obtaining the Bayes estimators of the parameters and reliability function. In order to compare the Bayes estimators under squared error loss, LINEX and SQUAREX loss functions respectively and the maximum likelihood estimator of the parameters and reliability function, Monte Carlo simulations are used.

• Communications for Statistical Applications and Methods
• /
• v.12 no.1
• /
• pp.19-27
• /
• 2005

The present paper obtains Bayes estimators for the probability p(Y

• The Korean Journal of Applied Statistics
• /
• v.15 no.1
• /
• pp.139-151
• /
• 2002

This study is concerned with model selection and diagnostics for nonlinear regression model through Bayes factor. In this paper, we use informative prior and simulate observations from the posterior distribution via Markov chain Monte Carlo. We propose the Laplace approximation method and apply the Laplace-Metropolis estimator to solve the computational difficulty of Bayes factor.

• Bulletin of the Korean Mathematical Society
• /
• v.58 no.6
• /
• pp.1419-1443
• /
• 2021

The asymmetric Laplace distribution arises as a limiting distribution of geometric summations of independent and identically distributed random variables with finite second moments. The main purpose of this paper is to study the weak limit theorems for geometric summations of independent (not necessarily identically distributed) random variables together with convergence rates to asymmetric Laplace distributions. Using Trotter-operator method, the orders of approximations of the distributions of geometric summations by the asymmetric Laplace distributions are established in term of the "large-𝒪" and "small-o" approximation estimates. The obtained results are extensions of some known ones.