• The Korean Journal of Applied Statistics
• /
• v.26 no.5
• /
• pp.697-712
• /
• 2013

When analyzing repeated binary data, the generalized estimating equations(GEE) approach produces consistent estimates for regression parameters even if an incorrect working correlation matrix is used. However, time-varying covariates experience larger changes in coefficients than time-invariant covariates across various working correlation structures for finite samples. In addition, the GEE approach may give biased estimates under missing at random(MAR). Weighted estimating equations and multiple imputation methods have been proposed to reduce biases in parameter estimates under MAR. This article studies if the two methods produce robust estimates across various working correlation structures for longitudinal binary data with time-varying covariates under different missing mechanisms. Through simulation, we observe that time-varying covariates have greater differences in parameter estimates across different working correlation structures than time-invariant covariates. The multiple imputation method produces more robust estimates under any working correlation structure and smaller biases compared to the other two methods.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2010.05a
• /
• pp.357-361
• /
• 2010

미래의 연별 최대 강수량 예측의 정확성을 향상시키는데 역사적 자료가 도움이 된다는 많은 연구 결과가 있었다. 관측의 오차와 자료의 손실로 역사자료를 이용한 강수 예측 방법은 절단자료의 분석을 중심으로 연구되었다. 대표적인 역사자료의 이용방법으로 조건부 적률을 이용한 B17B [Interagency Committee in Water Data, 1982], 조건부적률과적률 관계식을 이용한 Expected Moment Algorithm(EMA) [Cohn et al.;1997], 조건부 확률가중적률을 이용한 Partial Probability Weighted Moment (PPWM)[Wang ; 1991] 방법이 있다. 본 연구에서는 역사적 자료를 반영하는 방법에 있어 B17B와 EMA의 관계를 밝히고 그러한 관계가 PPWM에 동일하게 적용할 수 있음을 보였다. 우리는 B17B와 EMA의 관계를 적률방정식으로 표현하였고 PPWM에서 확률가중 적률 방정식을 정의함으로써 PPWM을 확장하였다. 본 연구에서 제안한 새로운 역사 자료를 이용한 강수예측 방법론을 Expected Probability Weighted Momemt (EPWM) 방법이라고 부르고 그 예측 방법의 성능을 다른 예측방법과 시뮬레이션 결과를 통해 비교하였다. 역사 자료 방법론의 비교는 Generalized Extreme Value (GEV) 분포를 이용하여 이루어졌으며, 각 방법론은 GEV분포의 형태모수(shape parameter)따라 다른 특성을 나타난다는 것을 보였다. 뿐만 아니라 여기서 제안한 EPWM 방법은 대부분의 경우에 좋은 추정량을 준다는 것을 보였다.

• The Korean Journal of Applied Statistics
• /
• v.27 no.2
• /
• pp.331-343
• /
• 2014

We propose two estimating procedures to analyze clustered interval-censored data with an informative cluster size based on a marginal model and investigate their asymptotic properties. One is an extension of Cong et al. (2007) to interval-censored data and the other uses the within-cluster resampling method proposed by Hoffman et al. (2001). Simulation results imply that the proposed estimators have a better performance in terms of bias and coverage rate of true value than an estimator with no adjustment of informative cluster size when the cluster size is related with survival time. Finally, they are applied to lymphatic filariasis data adopted from Williamson et al. (2008).

• Proceedings of the Korea Water Resources Association Conference
• /
• 2007.05a
• /
• pp.1209-1213
• /
• 2007

하천 수치모델링을 통한 흐름, 오염물질 거동, 지형변화 해석 등은 효율적인 하천 수질 관리를 위해서 상당히 중요한 부분은 차지한다. 수질이나 지형변화를 보다 정확하게 예측하기 위해서는 하천 흐름 예측의 정확도 향상이 중요한 역할을 하게 된다. 본 연구는 평면 이차원 하상변동 및 수질예측 수치모형인 KU-RLMS 모형을 이용하여 낙동강 상류의 반변천 합류부의 흐름 특성을 규명하고, 수질 모형을 수행하기 위한 흐름 계산 결과를 제공하기 위해 수행하였다. KU-RLMS 모형은 하천 및 저수지의 국부적인 수리, 수질, 유사이동 해석을 위해 개발된 평면 이차원 비정상 수치모형이다. 직사각형 격자를 사용하는 유한차분법의 단점을 보완하기 위해, 흐름 계산을 위한 지배방정식은 3차원 Reynolds 방정식으로부터 수심적분된 2차원 연속방정식과 운동량방정식을 불규칙한 경계를 현실적으로 모사할 수 있는 직교곡선 좌표계로 변환한 방정식을 사용한다. 수치모형 적용을 위한 현황분석으로 안동 및 임하 조정지댐의 방류량, 안동 수위관측소의 자료를 분석하였다. 흐름 모형을 보정하기 위해 안동대교 지점에서 횡유속 분포를 측정하였으며, 이 결과를 사용하여 흐름 모형의 매개변수인 Manning 계수와 공간가중계수를 추정 및 검증하였다. 안동다목적댐과 임하다목적댐의 방류량을 고려하여 수치모의조건을 결정하였으며, 각 조건에 대한 흐름 변화 특성을 분석하였다.

• The Korean Journal of Applied Statistics
• /
• v.29 no.7
• /
• pp.1155-1163
• /
• 2016

This paper deals with an estimation procedure of variance components in a mixed effects model by projections. Projections are used to obtain sums of squares instead of using reductions in sums of squares due to fitting both the assumed model and sub-models in the fitting constants method. A projection matrix can be obtained for the residual model at each step by a stepwise procedure to test the hypotheses. A weighted least squares method is used for the estimation of fixed effects. Satterthwaite's approximation is done for the confidence intervals for variance components.

• Journal of the Korean Data and Information Science Society
• /
• v.8 no.2
• /
• pp.153-161
• /
• 1997

A Korean exporience life table(male) is constructed by using a mixture of weighted moving average(WMA) model and Gompertz' parametric survival model based on 25,000,000 insured of major 6 life insurance companies from 1988 to 1992. The graduated values are taken as those which minimize the composite measure of fittness and smoothness. Moreover, we propose closed form estimators for three parameters of Gompertz' model.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2012.05a
• /
• pp.464-464
• /
• 2012

수문분야에 있어서 빈도해석의 목적은 특정 재현기간에 대한 발생 가능한 수문량의 규모를 파악하는데 있으며, 빈도해석의 정확도는 적합한 확률분포모형의 선택과 매개변수 추정방법에 의존하게 된다. 일반적으로 각 확률분포모형의 특성을 대표하는 매개변수를 추정하기 위해서는 모멘트 방법, 확률가중 모멘트 방법, 최대우도법 등을 이용하게 된다. 모멘트 방법에 의한 매개변수 추정은 해를 구하기 위한 과정이 단순한 반면, 비대칭형의 왜곡된 분포를 갖는 자료들에 대해서는 부정확한 결과를 나타내게 된다. 확률가중 모멘트 방법은 표본의 크기가 작거나 왜곡된 자료일 경우에도 비교적 안정적인 결과를 제공하는 반면, 확률 가중치가 정수로만 제한되는 단점을 갖고 있다. 그리고 대수 우도함수를 이용하여 매개변수를 추정하게 되는 최우도법은 가장 효율적인 매개변수 추정치를 얻을 수 있는 것으로 알려져 있으나, 비선형 연립방정식으로 표현되는 해를 구하기 위해서는 Newton-Raphson 방법을 사용하는 등 절차가 복잡하며, 때로는 수렴이 되지 않아 해룰 구하지 못하는 경우가 발생되게 된다. 이에 반해, 최근의 Genetic Algorithm, Ant Colony Optimization 및 Simulated Annealing과 같은 Meta-Heuristic Algorithm들은 복잡합 공학적 최적화 문제 있어서 효율적인 대안으로 주목받고 있으며, Hassanzadeh et al.(2011)에 의해 수문학적 빈도해석을 위한 매개변수 추정에 있어서도 그 적용성이 검증된바 있다. 본 연구의 목적은 연 최대강수 자료의 빈도해석에 적용되는 확률분포모형들의 매개변수 추정을 위해 Meta-Heuristic Algorithm을 적용하고자 함에 있다. 따라서 본 연구에서는 매개변수 추정을 위한 방법으로 Genetic Algorithm 및 Harmony Search를 적용하였고, 그 결과를 최우도법에 의한 결과와 비교하였다. GEV 분포를 이용하여 Simulation Test를 수행한 결과 Genetic Algorithm을 이용하여 추정된 매개변수들은 최우도법에 의한 결과들과 비교적 유사한 분포를 나타내었으나 과도한 계산시간이 요구되는 것으로 나타났다. 하지만 Harmony Search를 이용하여 추정된 매개변수들은 최우도법에 의한 결과들과 유사한 분포를 나타내었을 뿐만 아니라 계산시간 또한 매우 짧은 것으로 나타났다. 또한 국내 74개소의 강우관측소 자료와 Gamma, Log-normal, GEV 및 Gumbel 분포를 이용한 실증연구에 있어서도 Harmony Search를 이용한 매개변수 추정은 효율적인 매개 변수 추정치를 제공하는 것으로 나타났다.

• Water for future
• /
• v.24 no.1
• /
• pp.109-117
• /
• 1991

The governing equation of flow in porous media is developed on the bases of the continuity equation of fluid for transient flow through a saturated-unsaturated zone, and substitution of Darcy's law. The numerical solutions are obtained by finite element method based on the Galerkin principles weighted residuals. The analysis are carried out by using the unsteady storm data observed and rainfall intensities which are obtained by using the rainfall excess model in considering of the initial losses. The functional relationships between the hydraulic conductivity, capillary pressure head and volumetric water content are applied to the flow of water through unsaturated soil varied with changes of water content.

• The Journal of the Acoustical Society of Korea
• /
• v.21 no.4
• /
• pp.387-393
• /
• 2002

The enhancement technique of noise signal using mixture HFM (Midden Filter Model) are proposed. Given the parameters of the clean signal and noise, noisy signal is modeled by a linear state-space model with Markov switching parameters. Estimation of state vector is required for estimating original signal. The estimation procedure is based on mixture interacting multiple model (MIMM) and the estimator of speech is given by the weighted sum of parallel Kalman filters operating interactively. Simulation results showed that the proposed method offers performance gains relative to the previous results with slightly increased complexity.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2007.05a
• /
• pp.1915-1919
• /
• 2007

하천관리 및 수공구조물의 설계 등을 위한 빈도별 홍수량 산정 시 이용되는 방법에는 Clark 단위도법, Snyder 단위도법, KICT 단위도법 등이 있다. 이 중 대표적인 합성단위도법인 Clark 단위도법은 실무에서 가장 많이 사용되는 방법으로 꼽을 수 있으나, Clark 단위도법 매개변수인 집중시간$(T_c)$과 저류상수(K)를 추정하는 것은 단순하지 않다. Clark 단위도법의 매개변수 추정은 계측 유역인 경우와 미계측 유역인 경우로 나눌 수 있는데 대부분의 경우 미계측 유역의 도달시간 및 저류상수를 추정하는 방법으로 홍수량을 산정하고 있는 실정이다. 본 연구에서는 감천유역을 대상유역으로 선정하였으며, 김천수위표 지점을 유출지점으로 선정하여 홍수량을 산정하였다. 호우사상은 김천수위국의 수위-유량 곡선식으로 유출량의 검정이 가능한 년도를 고려하여 1998${\sim}$2005년 사이의 호우사상 중 분석이 가능한 것으로 판단되는 부항2와 김천 강우관측소의 시강우 자료를 선정하였다. 선정된 호우사상은 1999년, 2002년, 2004년에 대한 총 6개의 시우량 자료로 WAMIS(국가수자원종합정보시스템)와 한강홍수통제소(www.hrfco.go.kr)의 자료실에서 제공하는 한국수문조사연보 자료를 이용하였다. 선정된 호우사상과 WMS모형을 이용하여 산정된 유역면적 및 강우관측소 별 Thiessen 가중계수를 이용하여 HEC-1모형의 최적화기법으로 Clark 단위도법의 집중시간$(T_c)$과 저류상수(K)를 산정하였으며, 계측유역의 매개변수를 산정하는 방법 중 저류방정식을 이용하는 방법을 이용하여 저류상수를 산정하여 비교하였다. 집중시간의 경우 1999년과 2002년에 비해 2004년의 결과값이 작게 산정되었으며, 저류상수의 경우 2004년의 결과값이 크게 산정되었다. 그러나, 저류방정식을 이용하여 저류상수를 결정하는 방법으로 산정한 저류상수는 비교적 일관성 있게 산정되었다.