• Journal of the Korean Data and Information Science Society
• /
• v.7 no.2
• /
• pp.247-256
• /
• 1996

In this paper, we applied bayesian hierarchical model to analyze the field mice example introduced by Demster et al.(1981). For this example, we use Gibbs sampler method to provide the posterior mean and compared it with LSE(Least Square Estimator) and MLR(Maximum Likelihood estimator with Random effect) via the EM algorithm.

• Journal of the Korean Data and Information Science Society
• /
• v.20 no.5
• /
• pp.787-796
• /
• 2009

A random coefficient model is applied when times of the repeated measurements are not fixed in experiments with respect to the subjects. The procedures of the inference of a random coefficient model are same as those of a mixed model. Diet experimental data was used for applying the random coefficient model. Various random coefficient models are investigated for the experimental data, and are compared each other. Finally, optimal random coefficient model would be selected. It resulted from the analysis that for the fixed effect factor, the baseline, treatment, height, and time effect were very significant. The treatment effect of the diet foods and exercises were more effective in losing weight than the effect of the diet foods only. The fixed cubic time effect was very significant. The variance components corresponding to the subject effect, linear time effect, quadratic time effect, and cubic time effect of the random coefficients are all positive. When quartic time effect was added as random coefficients the model did not converge. Thus random coefficients up to the cubic terms was considered as the optimal model.

• Journal of Navigation and Port Research
• /
• v.37 no.5
• /
• pp.453-461
• /
• 2013

This paper describes evaluation procedures and experimental results for the estimation of Cumulative Distribution Functions (CDF) giving best-fit to the sample data in the Probability based risk Evaluation Techniques (PET) which is to assess the risks of a small-sized sea floater. The CDF in the PET is to provide the reference values of risk acceptance criteria which are to evaluate the risk level of the floater and, it can be estimated from sample data sets of motion response functions such as Roll, Pitch and Heave in the floater model. Using Maximum Likelihood Estimates and with the eight kinds of regulated distribution functions, the evaluation tests for the CDF having maximum likelihood to the sample data are carried out in this work. Throughout goodness-of-fit tests to the distribution functions, it is shown that the Beta distribution is best-fit to the Roll and Pitch sample data with smallest averaged probability errors $\bar{\delta}(0{\leq}\bar{\delta}{\leq}1.0)$ of 0.024 and 0.022, respectively and, Gamma distribution is best-fit to the Heave sample data with smallest $\bar{\delta}$ of 0.027. The proposed method in this paper can be expected to adopt in various application areas estimating best-fit distributions to the sample data.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2015.05a
• /
• pp.186-186
• /
• 2015

본 연구에서는 해안도시(부산광역시 수영구) 지역의 토지이용도와 미기상인자를 고려하여 증발산량을 산정하였으며, 증발산량 변동에 대한 미기상인자의 영향성을 구명하였다. 수영구 지역의 토지이용도와 미기상인자는 2001년 12월부터 2011년 11월에 관측된 일별 자료를 사용하였다. 토지이용도는 불투수(건물, 도로 등) 및 산림(임야), 초지(논밭, 공원 등), 수계(하천, 호수 등) 지역으로 분류하였으며, 4개 지역 특성을 고려한 최적의 추정식을 적용하여 증발산량을 산정하였다. 수영구 지역의 전체 증발산량은 4개 지역에서 산정된 증발산량에 토지이용 비율을 곱하여 구하였다. 연간 증발산량 변동은 1월부터 7월까지 증가하다가 8월부터 12월까지 감소하는 형태를 보였다. 수영구 지역에서 증발산량은 강수량의 약 13.3% 정도이었으며, 이는 연구지역의 72%에 해당하는 불투수 지역에서 배수로를 통한 물의 유출이 강우 발생 후 짧은 시간 동안 다량 발생하였기에 지속적인 증발산이 가능한 잠재수량의 저유량이 적었기 때문이다. 증발산량과 미기상인자 간의 상관분석을 수행하였으며, 증발산량과 이슬점 온도의 상관계수가 0.63으로 가장 높았다. 증발산량에 대한 기온 및 강수량, 순복사 인자의 상관계수는 0.5 이상으로 양의 상관성을, 기압 및 일조시간은 0.5 이상의 음의 상관성을 보였다. 증발산량에 대한 상관계수가 0.5 이상인 미기상인자(이슬점온도와 기온, 순복사, 기압, 강수량)에 대한 회귀 분석을 수행하였다. 이슬점온도와 기온, 순복사, 기압에 대한 증발산량 회귀함수 그래프는 강수의 유무에 따라 2가지 경향을 보였다. 이슬점온도에 따른 증발산량 회귀함수는 강수 발생 시에는 $ET=0.004x+0.7$, 무강수 시에는 $ET=0.25{\times}e^{0.04x}$로 추정되었으며, 결정계수는 각각 0.48과 0.96 정도로서 무강수 시에 높게 나타났다. 기온에 따른 증발산량 회귀함수는 강수 발생 시에는 $ET=0.004x+0.53$, 무강수 시에는 $ET=0.13{\times}e^{0.06x}$로 추정되었으며, 결정계수는 각각 0.39와 0.89 정도로서 무강수 시에 높게 나타났다. 순복사에 따른 증발산량 회귀함수는 강수 발생 시에는 $ET=0.79x+0.49$, 무강수 시에는 $ET=0.22x+0.03$로 추정되었으며, 결정계수는 각각 0.34와 0.89 정도로서 무강수 시에 높게 나타났다. 기압에 따른 증발산량 회귀함수는 강수 발생 시에는 $ET=-0.04x+37.91$, 무강수 시에는 $ET=5.18{\times}10^{22}{\times}e^{-0.05x}$로 추정되었으며, 결정계수는 각각 0.25와 0.45 정도로 나타났다. 강수량에 따른 증발산량 회귀함수는 $ET=0.23lnx+0.90$으로 추정되었으며, 결정계수 0.61정도 나타났다.

• Environmental and Resource Economics Review
• /
• v.13 no.1
• /
• pp.33-47
• /
• 2004

Even though the price of extracted but unprocessed coal has been available in Korea, the use of it as scarcity index would be inappropriate because of price subsidy. Following Halvorsen and Smith(1984), Kim and Lee(2002) derived estimates of the shadow price of unextracted coal by estimating the restricted cost function and differentiating with respect to the quantity of coal extracted. In Korea, however, due to the limited data the capital prices have been computed inconsistently case by case without relying on the robust formula like the Christensen-Jorgenson methodology used in US, which could result in biased estimators of the restricted cost function. In the paper the shadow prices of the resources in situ are obtained by measuring an input distance function defined by Shephard (1970), which requires only the data on the quantities of inputs and output. Empirical results for the Korean coal mining industry show that these shadow prices as a coal scarcity have increased fast by approximately three times in comparisons with those obtained by Kim and Lee.

• The Korean Journal of Applied Statistics
• /
• v.36 no.4
• /
• pp.309-322
• /
• 2023

The Hawkes process is a point process with self-exciting characteristics. It has been mainly used to describe seismic phenomena in which aftershocks occur due to the main earthquake. Recently, it has been used to explain various phenomena with self-exciting properties, such as the spread of infectious diseases and the spread of news on SNS. The Hawkes process can be flexibly modified according to the characteristics of events by using various types of excitation functions. Since it is difficult to implement a maximum likelihood estimator numerically, estimation methods have been improved until recently. In this paper, the conditional intensity function and excitation function are explained to describe the Hawkes process. Then, existing examples of Hawkes processes used in seismic, epidemiological, criminal, and financial fields are described and estimation methods are introduced. I analyze earthquakes that occurred in gyeongsang-do, Korea from November 2017 to December 2022, using R package ETAS.

• KDI Journal of Economic Policy
• /
• v.17 no.3
• /
• pp.3-54
• /
• 1995

외환(外換) 및 자본자유화(資本自由化)의 진전에 따라 환율(換率)의 변동폭이 커지면서 환율(換率)이 국제수지(國際收支)에 미치는 영향에 대한 관심이 높아지고 있다. 환율(換率)이 국제수지(國際收支)에 미치는 효과는 일차적으로 우리나라 수출입함수(輸出入函數)의 가격탄력성(價格彈力性)에 의해 결정되므로 정확한 탄성치(彈性値)를 추정하는 일은 이론적으로나 정책적으로 매우 중요한 의미를 갖는다. 본고에서는 우리나라 수입함수(輸入函數)의 소득(所得) 및 가격탄성치(價格彈性値) 추정에 있어서 사용 가능한 두가지의 가격지수(價格指數), 즉 수입단가(輸入單價)와 수입물가(輸入物價) 중 어느 것을 사용하는 것이 바람직한 것인지를 살펴보고 이중 이론적으로 우월한 지수를 이용하여 수입함수(輸入函數)의 소득(所得) 및 가격탄성치(價格彈性値)를 추정하고자 하였다. 분석으로부터 얻는 중요한 결론(結論)은, 우리나라의 경우 수입함수(輸入函數)의 분석에 있어서 수입단가(輸入單價)보다는 수입물가(輸入物價)를 이용하는 것이 더 적합하리라는 점이다. 이는 또한 우리나라 실질(實質)GNP의 추계에 있어서 현재의 수입단가(輸入單價) 이용방식에서 수입물가(輸入物價) 이용방식으로 편제방법을 바꾸는 것이 바람직하다는 것을 의미한다. 이같은 결론은 이론적, 실증적인 변에서 아래와 같은 네가지 이유에 근거한다. 먼저 이론적인 면에서는 첫째, 수입단가(輸入單價)가 물가지수(物價指數)로서 갖추어야 할 통계적 일관성을 갖추지 못하고 있다는 점을 지적할 수 있다. 즉 수입단가(輸入單價)는 어느 일정 규격조건(規格條件)의 상품가격(商品價格)을 일관성 있게 추적하는 것이 아니므로 잡다한 품질규격의 평균가격(平均價格)으로 나타날 가능성이 높으며 이러한 변에서 가격변화(價格變化)를 정확히 반영하지 못한다. 둘째, 수입단가지수(輸入單價指數)는 가격변화(價格變化)가 전혀 없을 때에도 총수입(總輸入)의 품목별 구성(構成)이 바뀌면 변동하는 속성을 지니고 있어 가격변동(價格變動) 상황(狀況)을 오도할 가능성을 지니고 있다. 다음으로 실증적인 면에서는 첫째, 수입단가(輸入單價)를 이용한 수입함수(輸入函數)에서는 가격변수(價格變數)의 계수(係數)가 이론상 기대되는 것과는 반대되는 부호를 가짐으로써 수입함수(輸入函數))로서의 기본적인 요건을 만족시키지 못하는 것으로 나타나고 있다. 둘째, 단가지수(單價指數)를 이용한 경우 소득(所得)과 상대가격(相對價格)을 회귀변수로 하는 수입함수(輸入函數)에서 장기안정적인 균충관계(均衝關係)가 존재하지 않은 것으로 판명되었다. 반면 수입물가(輸入物價)를 이용한 수입함수(輸入函數)에서는 각 변수의 계수(係數)가 이론적 기대치에 부합하였으며 장기안정적인 균충관계(均衝關係)가 성립하는 것으로 나타나고 있다. 한편 미국에서는 수출입단가지수(輸出入單價指數)의 문제점에 대한 오랜 연구 결과, 1989년부터 단가지수(單價指數)의 편제를 아예 중단하는 조치를 취한 바 있으며, GNP추계에 있어서도 종래의 수출입단가(輸出入單價) 대신 수출입물가(輸出入物價)를 이용하는 방식으로 그 편제방법을 변경한 바 있다.

• Proceedings of the Earthquake Engineering Society of Korea Conference
• /
• 2001.04a
• /
• pp.128-135
• /
• 2001

본 논문에서는 지반-구조물 상호작용계의 입력하중과 강성에 관련된 물성값들을 구조물에서 계측된 지진응답만을 사용하여 효과적으로 추정할 수 있는 새로운 방법을 제안하였는데, 제안된 방법은 비지계수 추정을 위한 목적함수가 입력하중에 독립적이기 때문에 구조물에서 계측된 지진응답만으로 미지계수의 추정이 가능하도록 되어있다. 본 연구에서 제안된 방법의 검증은 국제공동 연구의 일환으로 최근 대만의 화련에 건설된 대형지진시험 구조물에서 계측된 지진 응답을 사용하여 수행하였다. 추정된 입력하중과 지반과 구조물의 강성에 관련된 물성값들을 사용하여 계산된 지진응답의 계측치와 매우 잘 일치하여, 추정결과의 타당성을 검증할 수 있었다.

• Journal of Korean Society of Transportation
• /
• v.18 no.6
• /
• pp.47-61
• /
• 2000

차량추종모형은 차선변경모형과 함께 미시적 시뮬레이션 모형의 핵심 모형으로서 정확한 차량추종모형의 적용은 시뮬레이션 모형 분석결과의 신뢰성에 영향을 미친다. 기존의 차량추종모형은 대부분 가속상황과 감속상황에 동일한 모형을 적용하거나 자극부문과 민감도 부분은 거의 동일한 형태를 취하고 있으며, 가상자료를 토대로 모형이 개발되었다. 본 연구의 목적은 첫째, 모형을 추정하기 위한 연구체계를 구성하고, 둘째, 이러한 연구체계를 이용하여 신호교차로 미시적 시뮬레이션에 직접 적용이 가능한 차량추종모형을 정립하는 데 있다. 본연구에서는 이를 위해서 신호교차로에서 미시적 실측자료를 수집하였으며, 수집된 자료에서 통계적 기법을 통해서 모형 추정에 적합한 자료를 구축하였다. 현장에서 실측된 이산적인 원시자료는 차량추종모형 추정에 직접 이용할 수 없다. 따라서, 이산적 원시자료를 모형추정에 적합한 자료로 구성하기 위하여 비선형 곡선적합 알고리즘을 이용하여 연속적인 궤적함수를 구성하였고, 이를 통해서 모형추정에 필요한 선두차량과 후행차량의 현재시점과 반응시간 이후 시점에서의 위치, 순간 가속도와 순간속도 등을 도출하였다. 본 연구에서는 차량추종상황을 가속상황, 감속상황, 출발상황, 그리고 정지상황으로 구분하여 차량추종모형을 구축하였다. 가속상황과 감속상황에 대해서는 실측자료를 이용하여 모형을 추정하였으며, 출발상황과 정지상황에 대해서는 추정된 가/감속상황의 모형을 바탕으로 출발상황과 정지상황에서의 차량행태를 설명할 수 있는 모형을 구축하였다 또한, 민감도와 자극부분을 새롭게 정의하여 각 상황별로 미시적 실측자료를 이용하여 모형을 추정하였다. 이렇게 추정된 모형들 중에서 통계적 기법과 실측치와의 비교를 통해서 가장 적합한 모형을 선택하였다. 선택된 모형은 통계적 검정, 가상자료 그리고 실측치와의 비교를 통해서 분석하였으며, 분석결과 구축모형의 적용성은 우수한 것으로 분석되었다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2017.05a
• /
• pp.171-171
• /
• 2017

강우의 시간분포는 다양한 시간규모에 걸쳐 프랙털 또는 멀티프랙털 특성을 가지고 있음이 알려졌다. 기존의 연구는 주로 시간단위 이상의 프랙털 특성에 관한 것이었다. 실제로 극한 홍수를 가져오는 집중호우는 짧은 시간 규모에서 발생함에도, 이것에 대해서는 관측 자료가 제한되어 극소수의 실험적 연구만 가능했다. 본 연구에서는 기상청에서 제공한 고해상도(1분 단위) ASOS(Automated Synoptic Observation System) 자료를 이용하여, 강우 사상 안에서의 프랙털 특성을 분석해보았다. 대부분의 사상에서 단일 멱함수보다는 2개의 멱함수로 나누어지는 것이 밝혀졌으며, 나뉘는 시간 규모(T*)는 $3{\times}10$ 분으로 파악되었다. 이 시간 규모는 한 단위의 집중호우를 가져올 수 있는 구름크기의 물리적 상한과 관련이 있는 것으로 보인다. T*보다 작은 시간 규모에서의 멱함수 지수는 그 이후의 값보다 대체로 작은 것으로 나타났다. 이는 호우가 집중되는 기간의 변동성이, 강수가 물리적 한계에 도달한 이후보다 훨씬 작기 때문으로 보인다. 구체적인 멱함수의 지수는 강수의 발생과정과도 관련이 있을 것으로 추정된다.