• 제목/요약/키워드: metropolis-hastings algorithm

검색결과 43건 처리시간 0.018초

ASSESSING POPULATION BIOEQUIVALENCE IN A $2{\times}2$ CROSSOVER DESIGN WITH CARRYOVER EFFECT IN A BAYESIAN PERSPECTIVE

  • Oh Hyun-Sook
    • Journal of the Korean Statistical Society
    • /
    • 제35권3호
    • /
    • pp.239-250
    • /
    • 2006
  • A $2{\times}2$ crossover design including carryover effect is considered for assessment of population bioequivalence of two drug formulations in a Bayesian framework. In classical analysis, it is complex to deal with the carryover effect since the estimate of the drug effect is biased in the presence of a carryover effect. The proposed method in this article uses uninformative priors and vague proper priors for objectiveness of priors and the posterior probability distribution of the parameters of interest is derived with given priors. The posterior probabilities of the hypotheses for assessing population bioequivalence are evaluated based on a Markov chain Monte Carlo simulation method. An example with real data set is given for illustration.

BAYESIAN HIERARCHICAL MODEL WITH SKEWED ELLIPTICAL DISTRIBUTION

  • Chung, Youn-Shik;Dipak K. Dey;Yang, Tae-Young;Jang, Jung-Hoon
    • Journal of the Korean Statistical Society
    • /
    • 제32권4호
    • /
    • pp.425-448
    • /
    • 2003
  • Meta-analysis refers to quantitative methods for combining results from independent studies in order to draw overall conclusions. We consider hierarchical models including selection models under a skewed heavy tailed error distribution proposed originally by Chen et al. (1999) and Branco and Dey (2001). These rich classes of models combine the information of independent studies, allowing investigation of variability both between and within studies, and incorporate weight function. Here, the testing for the skewness parameter is discussed. The score test statistic for such a test can be shown to be expressed as the posterior expectations. Also, we consider the detail computational scheme under skewed normal and skewed Student-t distribution using MCMC method. Finally, we introduce one example from Johnson (1993)'s real data and apply our proposed methodology. We investigate sensitivity of our results under different skewed errors and under different prior distributions.

Bayesian analysis for the bivariate Poisson regression model: Applications to road safety countermeasures

  • Choe, Hyeong-Gu;Lim, Joon-Beom;Won, Yong-Ho;Lee, Soo-Beom;Kim, Seong-W.
    • Journal of the Korean Data and Information Science Society
    • /
    • 제23권4호
    • /
    • pp.851-858
    • /
    • 2012
  • We consider a bivariate Poisson regression model to analyze discrete count data when two dependent variables are present. We estimate the regression coefficients as sociated with several safety countermeasures. We use Markov chain and Monte Carlo techniques to execute some computations. A simulation and real data analysis are performed to demonstrate model fitting performances of the proposed model.

변곡 S-형 소프트웨어 신뢰도성장모형의 베이지안 모수추정 (Bayesian Estimation for Inflection S-shaped Software Reliability Growth Model)

  • 김희수;이종형;박동호
    • 품질경영학회지
    • /
    • 제37권4호
    • /
    • pp.16-22
    • /
    • 2009
  • The inflection S-shaped software reliability growth model (SRGM) proposed by Ohba(1984) is one of the most commonly used models and has been discussed by many authors. The main purpose of this paper is to estimate the parameters of Ohba's SRGM within the Bayesian framework by applying the Markov chain Monte Carlo techniques. While the maximum likelihood estimates for these parameters are well known, the Bayesian method for the inflection S-shaped SRGM have not been discussed in the literature. The proposed methods can be quite flexible depending on the choice of prior distributions for the parameters of interests. We also compare the Bayesian methods with the maximum likelihood method numerically based on the real data.

평균회귀확률과정을 이용한 2요인 사망률 모형 (A Two Factor Model with Mean Reverting Process for Stochastic Mortality)

  • 이강수;조재훈
    • 응용통계연구
    • /
    • 제28권3호
    • /
    • pp.393-406
    • /
    • 2015
  • 본 논문은 2요인(two-factor) 사망률 모형에 평균회귀모형(mean reverting process)을 적용하여 2요인의 확률적 변동을 모형화하여 사망률리스크(mortality risk)와 장수리스크(longevity risk)를 분석하였다. 최근 고령사회로 진입한 국가들에서 사망률 개선의 둔화가 관측되고 있는 시점에서 기존의 선형증가 또는 감소의 사망률 개선 모형을 보완함에 그 목적을 두었다. 영국의 1991~2015년 사망률 자료를 이용하여 제시한 모형의 모수를 메트로폴리스 알고리듬을 이용해 추정하였고 추정된 모수 값을 이용하여 다수 시뮬레이션을 통하여 장기간의 미래 사망률 예측값을 계산하였다. 평균회귀 모형의 특성으로 인해 약 60년의 시간이 지난 뒤부터는 사망률 개선이 거의 사라져 사망률이 일정한 값에 근접하였다. 사망률 개선이 둔화되는 현상이 관측되는 특정 집단(국가, 사회)의 경우 2요인 평균회귀 모형은 장기간 사망률 예측방법의 대안으로 간주될 것으로 기대되며, 모형의 응용으로서 평균회귀율의 추정결과로부터 사망률 개선의 속도를 계량화하는 기준을 제시하였다. 끝으로, 2014년~2040 기간의 사망률 예측값을 이용하여 25년 만기 장수채권의 발행가격을 산출하였다.

Grid Method 기법을 이용한 베이지안 비정상성 확률강수량 산정 (Bayesian Nonstationary Probability Rainfall Estimation using the Grid Method)

  • 곽도현;김광섭
    • 한국수자원학회논문집
    • /
    • 제48권1호
    • /
    • pp.37-44
    • /
    • 2015
  • 본 연구에서는 Grid method를 사용하여 베이지안 비정상성 확률강우량 산정 모형을 확립하였다. 강우 극치자료의 분포로 Gumbel 분포를 채택하였으며, 분포형의 매개변수에 사전분포를 적용하고, 사전분포에 포함된 매개변수에는 초사전 분포를 적용하여 계층적 베이지안 모형을 구성하였다. Grid method는 매개변수의 발생가능 전 구간에 대하여 확률적으로 더 높은 뒷받침이 있는 하위 구간에서 난수를 직접 생성하여 집합을 구성함으로써 잘못된 결과를 도출할 수 가능성이 높은 상황에서도 보다 정확한 매개변수의 추정을 가능케 하므로 매개변수의 추정과정에서 비표준분포로 나타나는 조건부 확률밀도함수를 통한 난수의 추출은 기존에 사용해 온 Metropolis Hastings 알고리즘이 아닌 Grid method를 사용하였다. 개발된 모형은 서울의 1973년부터 2012년까지의 시강우자료를 이용하여 미래에 대한 재현기간에 따른 확률강수량을 산정하였으며, 그 결과로 기존 정상성 가정에 비해 목표연도에 따라 5%에서 8%정도의 증가율을 나타냈다.

불확실성이 고려된 비용-편익분석 기법을 도입한 최적설계홍수량 산정 (Calculation of optimal design flood using cost-benefit analysis with uncertainty)

  • 김상욱;최광배
    • 한국수자원학회논문집
    • /
    • 제55권6호
    • /
    • pp.405-419
    • /
    • 2022
  • 홍수피해를 최소화하기 위한 수공구조물의 적정 규모 결정을 위해 사용되는 홍수빈도분석에는 통계적 분석절차에 따른 불확실성이 포함된다. 따라서 불확실성이 포함된 범주 내에서 가장 적절한 설계홍수량(design flood)을 결정하는 과정은 수공구조물의 최종단계에서 중요하게 다루어져야 하는 부분이나 이를 제시한 연구는 많지 않다. 비용-편익 분석기법을 홍수빈도분석 절차에 도입하여 구성되는 총 기대비용함수(total expected cost function)는 설계홍수량 중 최적설계홍수량(optimal design flood)을 결정하기 위한 새로운 접근방식이다. 이 절차는 UNCODE (UNcertainty COmpliant DEsign)로 명명되어 사용된 바 있으나, 국내에서는 아직 적용 결과가 소개되지 않고 있다. 따라서 본 연구에서는 UNCODE의 수학적 구성 절차를 소개함과 함께 북한강수계에 위치한 수력발전댐(화천댐, 춘천댐, 의암댐, 청평댐)의 연최대유입량을 사용하여 최적설계홍수량을 산정하고 이 결과를 기존 홍수빈도분석 결과와 비교하였다. 불확실성이 고려된 총 기대비용함수로부터 확률분포함수들(Gumbel 및 GEV)의 모수를 추출하는 과정에서 Metropolis-Hastings 알고리즘을 사용하여 불확실성의 범위를 추정하였으며, 비용-편익 분석기법에 사용되는 비용 및 피해함수는 수학적 구성의 편의성을 위하여 1차 선형함수로 가정되었다. 4개의 발전용댐, 2개의 확률분포 및 2개의 재현기간에 대하여 최적설계홍수량의 중앙값이 기존 홍수빈도분석 절차에 의해 산정된 설계홍수량보다 일정 정도 큰 값으로 산정됨을 알 수 있었다. 향후에는 본 연구에서 적용된 절차를 간단한 수식형태로 함수화하여 발전용댐 운영의 실무업무나 하천기본계획의 수립 등에 있어 비용-편익분석 기법의 적용성을 높이기 위한 연구가 필요할 것으로 판단된다.

How Does Internal Control Affect Bank Credit Risk in Vietnam? A Bayesian Analysis

  • PHAM, Hai Nam
    • The Journal of Asian Finance, Economics and Business
    • /
    • 제8권1호
    • /
    • pp.873-880
    • /
    • 2021
  • The purpose of this study is to investigate the impact of internal control on credit risk of joint stock commercial banks in Vietnam from 2007 to 2018. Furthermore, we specify bank-specific characteristics and macroeconomic conditions, and analyze how these factors affect credit risk of banks: the number of board members, the number of board members with banking or finance background as ratio of total board members, loans to total assets ratio, loans to deposit ratio, the number of days between the year-end and the publication of the financial statements, and the use of top four auditing firms proxy for five elements of internal control. By using the dataset of 30 Vietnamese joint stock commercial banks and Bayesian linear regression via Random-walk Metropolis Hastings algorithm, the results of this study show that five elements of internal control have a impact on bank credit risk, namely, control environment, risk assessment, control activities, information and communication, and monitoring activities. For factors of banks' characteristics, bank size and financial leverage have a negative impact on banks' credit risk, and bank age has a positive effect. For macroeconomic factors, inflation has a positive impact and economic growth has a negative impact on banks' credit risk.

그리드 단체 위의 디리슐레 분포에서 마르코프 연쇄 몬테 칼로 표집 (MCMC Algorithm for Dirichlet Distribution over Gridded Simplex)

  • 신봉기
    • 정보과학회 컴퓨팅의 실제 논문지
    • /
    • 제21권1호
    • /
    • pp.94-99
    • /
    • 2015
  • 비모수 베이스 통계학, 확률적 표집에 기반한 추론 등이 기계학습의 주요 패러다임으로 등장하면서 디리슐레(Dirichlet) 분포는 최근 다양한 그래프 모형 곳곳에 등장하고 있다. 디리슐레 분포는 일변수 감마 분포를 벡터 분포로 확장한 형태의 하나이다. 본 논문에서는 감마 분포를 갖는 임의의 자연수 X를 K개의 자연수의 합으로 임의 분할 할 때 각 부분의 크기 비율을 디리슐레 분포에서 표집하는 방법을 제안한다. 일반적으로 디리슐레 분포는 연속적인 (K-1)-단체(simplex) 위에 정의 되지만 자연수로 분할하는 표본은 자연수라는 조건 때문에 단체 내부의 이산 그리드 점에만 정의된다. 본 논문에서는 단체 위의 그리드 상의 이웃 점들의 확률 분포로부터 마르코프연쇄 몬테 칼로(MCMC) 제안 분포를 정의하고 일련의 표본들의 마르코프 연쇄를 구현하는 알고리듬을 제안한다. 본 방법은 마르코프 모델, HMM 및 준-HMM 등에서 각 상태별 시간 지속 분포를 표현하는데 활용 가능하다. 나아가 최근 제안된 전역-지역(global-local) 상태지속 분포를 동시에 모형화하는 감마-디리슐레 HMM에도 응용가능하다.

Structural modal identification and MCMC-based model updating by a Bayesian approach

  • Zhang, F.L.;Yang, Y.P.;Ye, X.W.;Yang, J.H.;Han, B.K.
    • Smart Structures and Systems
    • /
    • 제24권5호
    • /
    • pp.631-639
    • /
    • 2019
  • Finite element analysis is one of the important methods to study the structural performance. Due to the simplification, discretization and error of structural parameters, numerical model errors always exist. Besides, structural characteristics may also change because of material aging, structural damage, etc., making the initial finite element model cannot simulate the operational response of the structure accurately. Based on Bayesian methods, the initial model can be updated to obtain a more accurate numerical model. This paper presents the work on the field test, modal identification and model updating of a Chinese reinforced concrete pagoda. Based on the ambient vibration test, the acceleration response of the structure under operational environment was collected. The first six translational modes of the structure were identified by the enhanced frequency domain decomposition method. The initial finite element model of the pagoda was established, and the elastic modulus of columns, beams and slabs were selected as model parameters to be updated. Assuming the error between the measured mode and the calculated one follows a Gaussian distribution, the posterior probability density function (PDF) of the parameter to be updated is obtained and the uncertainty is quantitatively evaluated based on the Bayesian statistical theory and the Metropolis-Hastings algorithm, and then the optimal values of model parameters can be obtained. The results show that the difference between the calculated frequency of the finite element model and the measured one is reduced, and the modal correlation of the mode shape is improved. The updated numerical model can be used to evaluate the safety of the structure as a benchmark model for structural health monitoring (SHM).