• Title/Summary/Keyword: Markov chain Monte Carlo (MCMC)

Search Result 121, Processing Time 0.028 seconds

Bayesian Analysis of Dose-Effect Relationship of Cadmium for Benchmark Dose Evaluation (카드뮴 반응용량 곡선에서의 기준용량 평가를 위한 베이지안 분석연구)

  • Lee, Minjea;Choi, Taeryon;Kim, Jeongseon;Woo, Hae Dong
    • The Korean Journal of Applied Statistics
    • /
    • v.26 no.3
    • /
    • pp.453-470
    • /
    • 2013
  • In this paper, we consider a Bayesian analysis of the dose-effect relationship of cadmium to evaluate a benchmark dose(BMD). For this purpose, two dose-response curves commonly used in the toxicity study are fitted based on Bayesian methods to the data collected from the scientific literature on cadmium toxicity. Specifically, Bayesian meta-analysis and hierarchical modeling build an overall dose-effect relationship that use a piecewise linear model and Hill model, where the inter-study heterogeneity and inter-individual variability of dose and effect such as gender, age and ethnicity are accounted. Estimation of the unknown parameters is made by using a Markov chain Monte Carlo algorithm based user-friendly software WinBUGS. Benchmark dose estimates are evaluated for various cut-offs and compared with different tested subpopulations with with gender, age and ethnicity based on these two Bayesian hierarchical models.

Bayesian reliability estimation of bivariate Marshal-Olkin exponential stress-strength model

  • Chandra, N.;Pandey, M.
    • International Journal of Reliability and Applications
    • /
    • v.13 no.1
    • /
    • pp.37-47
    • /
    • 2012
  • In this article we attempted reliability analysis of a component under the stress-strength pattern with both classical as well as Bayesian techniques. The main focus is made to develop the theory for dealing the reliability problems in various circumstances for bivariate environmental set up in context of Bayesian paradigm. A stress-strength based model describes the life of a component which has strength (Y) and is subjected to stress(X). We develop the Bayes and moment estimators of reliability of a component for each of the three possible conditions, under the assumption that the two stresses (i.e. $X_1$ and $X_2$) on a component are dependent and follow a Bivariate exponential (BVE) of Marshall-Olkin distribution, the strength of a component (Y) following exponential distribution is independent of the stresses. The simulation study is performed with Markov Chain Monte Carlo technique via Gibbs sampler to obtain the estimates of Bayes estimators of reliability, are compared with moment estimators of reliabilities on the basis of absolute biases.

  • PDF

Seismic Reliability Analysis of Offshore Wind Turbine with Twisted Tripod Support using Subset Simulation Method (부분집합 시뮬레이션 방법을 이용한 꼬인 삼각대 지지구조를 갖는 해상풍력발전기의 지진 신뢰성 해석)

  • Park, Kwang-Yeun;Park, Wonsuk
    • Journal of the Computational Structural Engineering Institute of Korea
    • /
    • v.32 no.2
    • /
    • pp.125-132
    • /
    • 2019
  • This paper presents a seismic reliability analysis method for an offshore wind turbine with a twisted tripod support structure under earthquake loading. A three dimensional dynamic finite element model is proposed to consider the nonlinearity of the ground-pile interactions and the geometrical characteristics of the twisted tripod support structure where out-of-plane displacement occurs even under in-plane lateral loadings. For the evaluation of seismic reliability, the failure probability was calculated for the maximum horizontal displacement of the pile head, which is calculated from time history analysis using artificial earthquakes for the design return periods. The application of the subset simulation method using the Markov Chain Monte Carlo(MCMC) sampling is proposed for efficient reliability analysis considering the limit state equation evaluation by the nonlinear time history analysis. The proposed method can be applied to the reliability evaluation and design criteria development of the offshore wind turbine with twisted tripod support structure in which two dimensional models and static analysis can not produce accurate results.

Sensitivity analysis in Bayesian nonignorable selection model for binary responses

  • Choi, Seong Mi;Kim, Dal Ho
    • Journal of the Korean Data and Information Science Society
    • /
    • v.25 no.1
    • /
    • pp.187-194
    • /
    • 2014
  • We consider a Bayesian nonignorable selection model to accommodate the selection bias. Markov chain Monte Carlo methods is known to be very useful to fit the nonignorable selection model. However, sensitivity to prior assumptions on parameters for selection mechanism is a potential problem. To quantify the sensitivity to prior assumption, the deviance information criterion and the conditional predictive ordinate are used to compare the goodness-of-fit under two different prior specifications. It turns out that the 'MLE' prior gives better fit than the 'uniform' prior in viewpoints of goodness-of-fit measures.

Multivariable Bayesian curve-fitting under functional measurement error model

  • Hwang, Jinseub;Kim, Dal Ho
    • Journal of the Korean Data and Information Science Society
    • /
    • v.27 no.6
    • /
    • pp.1645-1651
    • /
    • 2016
  • A lot of data, particularly in the medical field, contain variables that have a measurement error such as blood pressure and body mass index. On the other hand, recently smoothing methods are often used to solve a complex scientific problem. In this paper, we study a Bayesian curve-fitting under functional measurement error model. Especially, we extend our previous model by incorporating covariates free of measurement error. In this paper, we consider penalized splines for non-linear pattern. We employ a hierarchical Bayesian framework based on Markov Chain Monte Carlo methodology for fitting the model and estimating parameters. For application we use the data from the fifth wave (2012) of the Korea National Health and Nutrition Examination Survey data, a national population-based data. To examine the convergence of MCMC sampling, potential scale reduction factors are used and we also confirm a model selection criteria to check the performance.

Bayesian Methods for Wavelet Series in Single-Index Models

  • Park, Chun-Gun;Vannucci, Marina;Hart, Jeffrey D.
    • 한국데이터정보과학회:학술대회논문집
    • /
    • 2005.04a
    • /
    • pp.83-126
    • /
    • 2005
  • Single-index models have found applications in econometrics and biometrics, where multidimensional regression models are often encountered. Here we propose a nonparametric estimation approach that combines wavelet methods for non-equispaced designs with Bayesian models. We consider a wavelet series expansion of the unknown regression function and set prior distributions for the wavelet coefficients and the other model parameters. To ensure model identifiability, the direction parameter is represented via its polar coordinates. We employ ad hoc hierarchical mixture priors that perform shrinkage on wavelet coefficients and use Markov chain Monte Carlo methods for a posteriori inference. We investigate an independence-type Metropolis-Hastings algorithm to produce samples for the direction parameter. Our method leads to simultaneous estimates of the link function and of the index parameters. We present results on both simulated and real data, where we look at comparisons with other methods.

  • PDF

Bayesian Variable Selection in the Proportional Hazard Model with Application to DNA Microarray Data

  • Lee, Kyeon-Eun;Mallick, Bani K.
    • Proceedings of the Korean Society for Bioinformatics Conference
    • /
    • 2005.09a
    • /
    • pp.357-360
    • /
    • 2005
  • In this paper we consider the well-known semiparametric proportional hazards (PH) models for survival analysis. These models are usually used with few covariates and many observations (subjects). But, for a typical setting of gene expression data from DNA microarray, we need to consider the case where the number of covariates p exceeds the number of samples n. For a given vector of response values which are times to event (death or censored times) and p gene expressions (covariates), we address the issue of how to reduce the dimension by selecting the significant genes. This approach enable us to estimate the survival curve when n < < p. In our approach, rather than fixing the number of selected genes, we will assign a prior distribution to this number. The approach creates additional flexibility by allowing the imposition of constraints, such as bounding the dimension via a prior, which in effect works as a penalty. To implement our methodology, we use a Markov Chain Monte Carlo (MCMC) method. We demonstrate the use of the methodology to diffuse large B-cell lymphoma (DLBCL) complementary DNA(cDNA) data.

  • PDF

Bayesian Inferences for Software Reliability Models Based on Beta-Mixture Mean Value Functions

  • Nam, Seung-Min;Kim, Ki-Woong;Cho, Sin-Sup;Yeo, In-Kwon
    • The Korean Journal of Applied Statistics
    • /
    • v.21 no.5
    • /
    • pp.835-843
    • /
    • 2008
  • In this paper, we investigate a Bayesian inference for software reliability models based on mean value functions which take the form of the mixture of beta distribution functions. The posterior simulation via the Markov chain Monte Carlo approach is used to produce estimates of posterior properties. Its applicability is illustrated with two real data sets. We compute the predictive distribution and the marginal likelihood of various models to compare the performance of them. The model comparison results show that the model based on the beta-mixture performs better than other models.

Multiple Comparisons for a Bivariate Exponential Populations Based On Dirichlet Process Priors

  • Cho, Jang-Sik
    • Journal of the Korean Data and Information Science Society
    • /
    • v.18 no.2
    • /
    • pp.553-560
    • /
    • 2007
  • In this paper, we consider two components system which lifetimes have Freund's bivariate exponential model with equal failure rates. We propose Bayesian multiple comparisons procedure for the failure rates of I Freund's bivariate exponential populations based on Dirichlet process priors(DPP). The family of DPP is applied in the form of baseline prior and likelihood combination to provide the comparisons. Computation of the posterior probabilities of all possible hypotheses are carried out through Markov Chain Monte Carlo(MCMC) method, namely, Gibbs sampling, due to the intractability of analytic evaluation. The whole process of multiple comparisons problem for the failure rates of bivariate exponential populations is illustrated through a numerical example.

  • PDF

Simulating phase transition phenomena of the unitary cell model

  • Kim, Dong-Hoh
    • Journal of the Korean Data and Information Science Society
    • /
    • v.20 no.1
    • /
    • pp.225-235
    • /
    • 2009
  • Lattice process models are used to explain phase transitions in statistical mechanics, a branch of physics. The Ising model, a specific form of lattice process model, was proposed by Ising in 1925. Since then, variants of the Ising model such as the Potts model and the unitary cell model have been proposed. Like the Ising model, it is believed that the more general models exhibit phase transitions on the critical surface, which is based on the mathematical equation. In statistical sense, phase transitions can be simulated through Markov Chain Monte Carlo (MCMC). We applied Swendsen-Wang algorithm, a block Gibbs algorithm, to a general lattice process models and we simulate phase transition phenomena of the unitary cell model.

  • PDF