The objective of this study was to analyze the flood stage considering the uncertainty caused by the river roughness coefficients and discharge. The methodology of this study involved the GLUE (Generalized Likelihood Uncertainty Estimation) to quantify the uncertainty bounds applying three different storm events. The uncertainty range of the roughness was 0.025~0.040. In case of discharge, the uncertainty stemmed from parameters in stage-discharge rating curve, if represents stage for discharge , which can be written as . Parameters in rating curve (, and ) were estimated by non-linear regression model and assumed by distribution. The range of parameters in rating curve was 5.138~18.442 for , -0.524~0.104 for and 2.427~2.924 for . By sampling 10,000 parameter sets, Monte Carlo simulations were performed. The simulated stage value was represented by 95% confidence interval. In storm event 1~3, the average bound was 0.39 m, 0.83 m and 0.96 m, respectively. The peak bound was 0.52 m, 1.36 m and 1.75 m, respectively. The recurrence year of each storm event applying the frequency analysis was 1-year, 10-year and 25-year, respectively.