• Title/Summary/Keyword: Gibbs sampling

Search Result 168, Processing Time 0.028 seconds

Accelerating Scanline Block Gibbs Sampling Method using GPU (GPU 를 활용한 스캔라인 블록 Gibbs 샘플링 기법의 가속)

  • Zeng, Dongmeng;Kim, Wonsik;Yang, Yong;Park, In Kyu
    • Proceedings of the Korean Society of Broadcast Engineers Conference
    • /
    • 2014.06a
    • /
    • pp.77-78
    • /
    • 2014
  • A new MCMC method for optimization is presented in this paper, which is called the scanline block Gibbs sampler. Due to its slow convergence speed, traditional Markov chain Monte Carlo (MCMC) is not widely used. In contrast to the conventional MCMC method, it is more convenient to parallelize the scanline block Gibbs sampler. Since The main part of the scanline block Gibbs sampler is to calculate message between each edge, in order to accelerate the calculation of messages passing in scanline sampler, it is parallelized in GPU. It is proved that the implementation on GPU is faster than on CPU based on the experiments on the OpenGM2 benchmark.

  • PDF

Sampling Based Approach to Bayesian Analysis of Binary Regression Model with Incomplete Data

  • Chung, Young-Shik
    • Journal of the Korean Statistical Society
    • /
    • v.26 no.4
    • /
    • pp.493-505
    • /
    • 1997
  • The analysis of binary data appears to many areas such as statistics, biometrics and econometrics. In many cases, data are often collected in which some observations are incomplete. Assume that the missing covariates are missing at random and the responses are completely observed. A method to Bayesian analysis of the binary regression model with incomplete data is presented. In particular, the desired marginal posterior moments of regression parameter are obtained using Meterpolis algorithm (Metropolis et al. 1953) within Gibbs sampler (Gelfand and Smith, 1990). Also, we compare logit model with probit model using Bayes factor which is approximated by importance sampling method. One example is presented.

  • PDF

Bayesian Analysis for Categorical Data with Missing Traits Under a Multivariate Threshold Animal Model (다형질 Threshold 개체모형에서 Missing 기록을 포함한 이산형 자료에 대한 Bayesian 분석)

  • Lee, Deuk-Hwan
    • Journal of Animal Science and Technology
    • /
    • v.44 no.2
    • /
    • pp.151-164
    • /
    • 2002
  • Genetic variance and covariance components of the linear traits and the ordered categorical traits, that are usually observed as dichotomous or polychotomous outcomes, were simultaneously estimated in a multivariate threshold animal model with concepts of arbitrary underlying liability scales with Bayesian inference via Gibbs sampling algorithms. A multivariate threshold animal model in this study can be allowed in any combination of missing traits with assuming correlation among the traits considered. Gibbs sampling algorithms as a hierarchical Bayesian inference were used to get reliable point estimates to which marginal posterior means of parameters were assumed. Main point of this study is that the underlying values for the observations on the categorical traits sampled at previous round of iteration and the observations on the continuous traits can be considered to sample the underlying values for categorical data and continuous data with missing at current cycle (see appendix). This study also showed that the underlying variables for missing categorical data should be generated with taking into account for the correlated traits to satisfy the fully conditional posterior distributions of parameters although some of papers (Wang et al., 1997; VanTassell et al., 1998) presented that only the residual effects of missing traits were generated in same situation. In present study, Gibbs samplers for making the fully Bayesian inferences for unknown parameters of interests are played rolls with methodologies to enable the any combinations of the linear and categorical traits with missing observations. Moreover, two kinds of constraints to guarantee identifiability for the arbitrary underlying variables are shown with keeping the fully conditional posterior distributions of those parameters. Numerical example for a threshold animal model included the maternal and permanent environmental effects on a multiple ordered categorical trait as calving ease, a binary trait as non-return rate, and the other normally distributed trait, birth weight, is provided with simulation study.

Estimation of Genetic Parameters for Carcass Traits in Hanwoo Steer (거세한우의 도체형질에 대한 유전모수 추정)

  • Yoon, H.B.;Kim, S.D.;Na, S.H.;Chang, U.M.;Lee, H.K.;Jeon, G.J.;Lee, D.H.
    • Journal of Animal Science and Technology
    • /
    • v.44 no.4
    • /
    • pp.383-390
    • /
    • 2002
  • The data were consisted of 1,262 records for carcass traits observed at Hanwoo steers from 1998 to 2001 at Namwon and Deakwanryung branch of National Livestock Research Institute, Rural Development Administration. Pedigrees of young bulls were traced back to search for magnifying inbreeding. Genetic parameters for carcass traits with Gibbs sampling in a threshold animal model were compared to estimates with REML algorithm in linear model. As the results, most of bulls were not inbred and sire pedigree group was non-inbred population. However, most of the bulls fell in some relationship with each other. Heritability estimates as fully posterior means by Gibbs samplers in threshold model were higher than those by REML in linear model. Furthermore, these estimates in threshold model using GS showed higher estimates than estimates using tested young bulls in previous study and same model. Heritability estimate by GS for marbling score was 0.74 and genetic correlation estimate between marbling score and body weight at slaughter was –0.44. Further study for correlation of breeding values between REML algorithm in linear model and Gibbs sampling algorithm in threshold model was needed.

Genetic Correlation of Reproductive Trait in Pigs by Parity (돼지 번식형질의 산차간 유전상관)

  • Cho, Kyu-Ho;Kim, Myung-Jick;Lee, Il-Joo;Kim, In-Cheul;Jeon, Gwang-Joo
    • Journal of Animal Science and Technology
    • /
    • v.50 no.4
    • /
    • pp.457-464
    • /
    • 2008
  • Heritability and genetic correlation for reproductive traits in Yorkshire pig breed were estimated using Bayesian method via Gibbs sampling. The data set consisted of 9,609 reproductive records at pig breeding farm in Korea. For estimating those parameters using Gibbs sampling, 5,000 cycles of ‘burn-in’ period were discarded among a total of 55,000 samples. Out of the remaining 50,000 samples, 5,000 estimates by each parameter were retained and used for analyses to avoid any correlation among adjacent samples. The reproductive trait considered in this study were total number of born piglets(TNB) and estimated by two different models. The estimated heritability and permanent environmental effect using Gibbs sampler were 0.12±0.020 and 10.9±1.63, respectively. Estimated genetic correlations considered parities as different traits were distributed from 0.99 to -0.13. Such results indicated that reproductive traits for sows should be considered as different traits.

A novel Metropolis-within-Gibbs sampler for Bayesian model updating using modal data based on dynamic reduction

  • Ayan Das;Raj Purohit Kiran;Sahil Bansal
    • Structural Engineering and Mechanics
    • /
    • v.87 no.1
    • /
    • pp.1-18
    • /
    • 2023
  • The paper presents a Bayesian Finite element (FE) model updating methodology by utilizing modal data. The dynamic condensation technique is adopted in this work to reduce the full system model to a smaller model version such that the degrees of freedom (DOFs) in the reduced model correspond to the observed DOFs, which facilitates the model updating procedure without any mode-matching. The present work considers both the MPV and the covariance matrix of the modal parameters as the modal data. Besides, the modal data identified from multiple setups is considered for the model updating procedure, keeping in view of the realistic scenario of inability of limited number of sensors to measure the response of all the interested DOFs of a large structure. A relationship is established between the modal data and structural parameters based on the eigensystem equation through the introduction of additional uncertain parameters in the form of modal frequencies and partial mode shapes. A novel sampling strategy known as the Metropolis-within-Gibbs (MWG) sampler is proposed to sample from the posterior Probability Density Function (PDF). The effectiveness of the proposed approach is demonstrated by considering both simulated and experimental examples.

Bayes Estimation for the Rayleigh Failure Model

  • Ko, Jeong-Hwan;Kang, Sang-Gil;Shin, Jae-Kyoung
    • Journal of the Korean Data and Information Science Society
    • /
    • v.9 no.2
    • /
    • pp.227-235
    • /
    • 1998
  • In this paper, we consider a hierarchical Bayes estimation of the parameter, the reliability and hazard rate function based on type-II censored samples from a Rayleigh failure model. Bayes calculations can be implemented easily by means of the Gibbs sampler. A numerical study is provided.

  • PDF

Monte Carlo Estimation of Multivariate Normal Probabilities

  • Oh, Man-Suk;Kim, Seung-Whan
    • Journal of the Korean Statistical Society
    • /
    • v.28 no.4
    • /
    • pp.443-455
    • /
    • 1999
  • A simulation-based approach to estimating the probability of an arbitrary region under a multivariate normal distribution is developed. In specific, the probability is expressed as the ratio of the unrestricted and the restricted multivariate normal density functions, where the restriction is given by the region whose probability is of interest. The density function of the restricted distribution is then estimated by using a sample generated from the Gibbs sampling algorithm.

  • PDF

Semiparametric Bayesian multiple comparisons for Poisson Populations

  • Cho, Jang Sik;Kim, Dal Ho;Kang, Sang Gil
    • Communications for Statistical Applications and Methods
    • /
    • v.8 no.2
    • /
    • pp.427-434
    • /
    • 2001
  • In this paper, we consider the nonparametric Bayesian approach to the multiple comparisons problem for I Poisson populations using Dirichlet process priors. We describe Gibbs sampling algorithm for calculating posterior probabilities for the hypotheses and calculate posterior probabilities for the hypotheses using Markov chain Monte Carlo. Also we provide a numerical example to illustrate the developed numerical technique.

  • PDF

Bayesian Analysis for Multiple Capture-Recapture Models using Reference Priors

  • Younshik;Pongsu
    • Communications for Statistical Applications and Methods
    • /
    • v.7 no.1
    • /
    • pp.165-178
    • /
    • 2000
  • Bayesian methods are considered for the multiple caputure-recapture data. Reference priors are developed for such model and sampling-based approach through Gibbs sampler is used for inference from posterior distributions. Furthermore approximate Bayes factors are obtained for model selection between trap and nontrap response models. Finally one methodology is implemented for a capture-recapture model in generated data and real data.

  • PDF