Effluent concentration estimation equations for wetland system were developed throught statistical analysis of treatment wetland experimental data. Existin g empirical equations were reviewed for thier accuracy with experimental data, and compared with the estimatin equations. About 70 experimental data sets were used for multiple regression, and variables include influent concentration, hydraulic loading rate, average daily air temperature , and plant coverage. The estimatin equations developed for BOD5 , SS ,T-P, and T-N predicted effluent concentrations moderately well, and coefficient fo determination ($R^2$) for them was 0.74 , 0.60, 0.59 and 0.58 respectively. The equations obtained from same data but excluding plant coverage showed relatively lower $R^2$ than the former case, and it was 0.66, 0.52, 0.41 and 0.57 respectively. The EPA, WPCF , and Kadlec and Knight equations worked poorly and $R^2$ for them was significantly lower than the estimation equation developed in the study. The reason might be that the existing equations were oversimplified that they did ot include important parameters such as air temperature and plant coverage. Therefore, developing reasonable estimation equations from experiment under realistic condition is highly recommended rather than using exiting estimation equations.