In this study, we proposed a residential load modeling method based on neuro-fuzzy inference system by considering of various harmonics. The developed method was implemented by using harmonic information, fundamental frequency and voltage which are essential input factors in conventional method. Thus, the proposed method makes it possible to effectively estimate load characteristics in power lines with harmonics. To show the effectiveness, the proposed method has been intensively tested with various dataset acquired under the different frequency and voltage and compared it with a conventional method based on neural networks.