• 제목/요약/키워드: Metropolis-Hastings 알고리즘

검색결과 18건 처리시간 0.024초

총 기대비용함수를 이용한 최적설계홍수량 결정 (Determination of optimal flood using total expected cost function)

  • 김상욱;최광배;서동일;전영일
    • 한국수자원학회:학술대회논문집
    • /
    • 한국수자원학회 2022년도 학술발표회
    • /
    • pp.287-287
    • /
    • 2022
  • 홍수빈도분석의 실용적 측면의 목적은 특정 재현기간에 대하여 발생 가능한 홍수량을 설계홍수량(design flood)으로 설정함으로써 댐, 제방, 배수시설, 하수관거 등의 치수기능을 가진 치수시설물이 설계홍수량 내에서 홍수로 인한 피해를 발생시키지 않도록 그 규모와 기능을 설계함에 있다. 특히 우리나라의 경우 유량자료의 부족으로 강우빈도분석을 수행하여 재현기간별 확률강우량을 먼저 산정하고 이를 강우-유출모형을 통해 확률홍수량으로 전환한 뒤 하천등급에 따른 재현기간 기준에 따라 설계홍수량을 산정하고 있다. 그러나 이와 같이 결정된 설계홍수량이 특정유역에서 발생될 수 있는 피해규모에 대해 얼마나 적정한 지의 여부를 과학적으로 판단하기 위한 연구는 국내·외에서 찾아보기 어려우며, 이러한 문제를 개선하기 위한 기초 이론을 제공하는 것이 본 연구의 가장 중요한 목표이다. 홍수빈도분석을 통해 산정된 설계홍수량의 적정성 여부를 과학적으로 판단하기 위해 최근에 진행된 해외의 몇몇 연구에서는 총 기대비용함수(total expected cost function)의 개발에 근거한 최적설계홍수량을 활용할 수 있음을 제안한 바 있다. 이 개념은 계획된 설계홍수량 이상에서 발생될 수 있는 피해함수(damage function) 및 기대피해함수(expected damage function)와 비용함수(cost function)가 결정되면, 이로부터 총 비용을 나타내는 총 기대비용함수(total expected cost function)을 도출하고 총 기대비용함수가 최소가 되는 최적설계홍수량(optimal design flood)을 산정하여 이를 계획된 설계홍수량(tentative design flood) 비교함으로써 계획된 설계홍수량의 적정성을 판단하는 과정을 기초이론으로 활용한다. 본 연구에서는 불확실성으로 발생되는 범위를 고려한 최적설계홍수량을 산정하기 위하여 Metropolis-Hastings 알고리즘을 사용하였으며, 자료의 종류에 따른 홍수량의 변화를 분석하기 위하여 년최대계열 및 부분시계열 자료를 각각 적용하였다. 한강유역에서 가평대성, 여주 및 한강대교 수위표 지점에서 측정된 자동관측유량장치에 의한 홍수량 자료를 활용하였으며, 최적설계홍수량이 기존 설계홍수량에 비해 크게 산정됨을 알 수 있었다.

  • PDF

보조 혼합 샘플링을 이용한 베이지안 로지스틱 회귀모형 : 당뇨병 자료에 적용 및 분류에서의 성능 비교 (Bayesian logit models with auxiliary mixture sampling for analyzing diabetes diagnosis data)

  • 이은희;황범석
    • 응용통계연구
    • /
    • 제35권1호
    • /
    • pp.131-146
    • /
    • 2022
  • 로지스틱 회귀 모형은 다양한 분야에서 범주형 종속 변수를 예측하거나 분류하기 위한 모형으로 많이 사용되고 있다. 로지스틱 회귀 모형에 대한 전통적인 베이지안 추론 기법으로 메트로폴리스-헤이스팅스 알고리즘이 많이 사용되었지만, 수렴의 속도가 느리고 제안 분포에 대한 적절성을 보장하기 어렵다. 따라서, 본 논문에서는 모형에 대한 베이지안 추론 방법으로 Frühwirth-Schnatter와 Frühwirth (2007)에서 제안된 보조 혼합 샘플링(auxiliary mixture sampling) 기법을 사용하였다. 이 방법은 모형의 선형성과 정규성을 만족시키기 위해 두 단계에 거쳐 잠재변수를 도입하며, 결과적으로 깁스 샘플링을 통한 추론을 가능하게 한다. 제안한 모형의 효과를 검증하기 위해 2020년 지역사회 건강조사 당뇨병 자료에 적용하여 메트로폴리스-헤이스팅스를 사용한 모형과 추론 결과를 비교 분석하였다. 또한, 다양한 분류 모형들과 본 논문에서 제안한 모형의 분류 성능을 비교한 결과 제안된 모형이 분류 분석에서도 좋은 성능을 보이는 것을 확인할 수 있었다.

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개의 재현기간에 대하여 최적설계홍수량의 중앙값이 기존 홍수빈도분석 절차에 의해 산정된 설계홍수량보다 일정 정도 큰 값으로 산정됨을 알 수 있었다. 향후에는 본 연구에서 적용된 절차를 간단한 수식형태로 함수화하여 발전용댐 운영의 실무업무나 하천기본계획의 수립 등에 있어 비용-편익분석 기법의 적용성을 높이기 위한 연구가 필요할 것으로 판단된다.

베이지안 포아송 모형을 적용한 자기-대조 환자군 연구에서의 약물상호작용 위험도 분석 (A Bayesian Poisson model for analyzing adverse drug reaction in self-controlled case series studies)

  • 이은채;황범석
    • 응용통계연구
    • /
    • 제33권2호
    • /
    • pp.203-213
    • /
    • 2020
  • 자기-대조 환자군(self-controlled case series; SCCS) 연구는 별도의 대조군 없이 환자의 비노출기간을 대조기간으로 설정하여 노출기간에 대한 상대적인 발생 위험도를 측정하는 역학 연구의 한 방법이다. 이 방법은 대조군을 선정할 때 발생하는 편의를 최소화할 수 있는 장점이 있어서 약물 복용 후 이상반응 발생 위험도를 측정하기 위한 방법으로 전통적으로 많이 사용되어왔다. 본 연구는 SCCS 연구를 바탕으로 두 개 이상의 약물을 동시에 사용했을 때 그 부작용의 위험이 어떻게 증가하는지 살펴보고자 한다. 마약성 진통제 유사체인 tramadol과 다빈도 병용 약물 간 약물상호작용에 대해 조건부 포아송 모형을 가정하고 분석하였다. 이때 베이지안 추론법을 사용하여 최대가능도추정량이 지니고 있는 과대적합 문제를 해결하며, 사전분포의 민감도를 측정하기 위해 정규 사전분포와 라플라스 사전분포를 가정하여 모형화하였다.

북한강 수계 수력발전댐의 최적설계홍수량 추정 (Estimating optimal flood of the hydroelectric dams on the Bukhangang River)

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

  • PDF

오차 계산 방식에 따른 사료용 벼 품종의 품종모수 추정치 불확도 비교 (Comparison between Uncertainties of Cultivar Parameter Estimates Obtained Using Error Calculation Methods for Forage Rice Cultivars)

  • 조영상;현신우;김광수
    • 한국농림기상학회지
    • /
    • 제25권3호
    • /
    • pp.129-141
    • /
    • 2023
  • 작물 모형은 작물의 유전적 특성을 나타내는 품종모수를 요구하며, 품종모수는 작물의 개별 품종별로 추정되어야 한다. 품종모수의 추정에는 고품질의 많은 생육 자료가 요구되지만, 자료의 생산에 상당한 비용이 필요하다. 비교적 낮은 품질의 가용성이 높은 자료를 활용하는 대신, 대량의 랜덤 모수를 생성하고 이를 평가하여 품종모수를 추정할 수 있다. 본 연구에서는 SIMPLE 작물 모델의 불확도를 최소화하기 위해 품종모수 추정 방식을 비교하고, 두 앙상블 방식과 대한 비교를 하였다. 모수 추정을 위한 Metropolis-Hastings (MH) 알고리즘에 대한 목적함수로 로그 가능도(log-likelihood: LL)와 generic composite similarity measure (GCSM)를 사용하였다. 또한 품종모수의 평균값을 사용한 예측(Epm)과 개별 모수들로부터 얻어진 추정값의 평균값(Eem)의 일치도를 분석하여 앙상블 방식에 따른 불확도 변화를 파악하였다. 국내에서 재배되는 사료용 벼 품종인 조우 벼와 영우 벼를 대상으로 품종모수를 추정하였다. 2013년, 2014년, 2016년에 대한 수원, 전주, 나주, 익산에 위치한 실험포장에서 얻은 수량 관측 자료를 사용하였다. 또한 2016년부터 2018년까지 수원에서 보고된 별도의 수량 관측 자료를 사용하였다. 목적함수에 따라 추정된 품종모수의 분포에 차이가 있었다. LL을 통해 얻은 품종모수는 GCSM으로 얻은 품종모수보다 좁은 범위에 분포하였다. 두 가지 앙상블 접근법은 통계적으로 유의한 차이가 나타나지 않음을 확인하였다. GCSM의 상대적으로 높은 불확도는 수용확률을 조정하여 낮출 수 있다고 사료되고, Epm의 결과는 기존과 다른 앙상블 방식을 통해 적은 연산을 통해 불확도를 낮출 수 있음을 보인다.

t-링크를 갖는 마코프 이항 회귀 모형을 이용한 인도네시아 어린이 종단 자료에 대한 베이지안 분석 (Bayesian inference of longitudinal Markov binary regression models with t-link function)

  • 심보현;정윤식
    • 응용통계연구
    • /
    • 제33권1호
    • /
    • pp.47-59
    • /
    • 2020
  • 본 논문에서는 마코프 이항 회귀 모형의 시차가 알려져 있거나 그렇지 않은 경우일 때, t-링크 함수를 갖는 종단적 마코프 이항 회귀 모형을 제시한다. 일반적으로, 이항 회귀 모형에서는 로직 모형이나 프로빗 모형이 주로 사용된다. t-링크 함수는 t 분포가 자유도가 커질수록 정규분포로 근사하기 때문에 프로빗 모형을 대신 더 많은 유연성을 위해 사용될 수 있다. 게다가 마코프 회귀모형은 종단 자료에 대해 사용될 수 있다. 우리는 마코프 회귀 모형의 시차를 결정하기 위해 베이지안 방법을 제시하고자 한다. 특히, 각 모델의 차수에 대해 알고 있는 경우에는 DIC를 기준으로 모델 비교를 실시하였다. 모델의 차수에 대해 모르는 경우에는 가능한 모델들의 사후 확률을 이용하였다. 복잡한 베이지안 계산을 해결하기 위하여 Albert와 Chib (1993), Kuo와 Mallick (1998)과 Erkanli 등 (2001)의 방법을 이용하여 모델을 재설정하였다. 제안하는 방법은 시뮬레이션 데이터와 Somer 등 (1984)에 의해 조사된 인도네시아 어린이 종단 데이터에 적용했다. 마코프 이항 회귀모형의 순서에 대해서 아는 경우와 모르는 경우를 각각 가정하여 최적의 모델을 알아보기 위해 MCMC 방법을 사용하였다. 또한, 매트로폴리스 해스팅 알고리즘의 수렴성을 점검하기 위해 Gelman과 Rubin의 진단을 이용했다.