This paper describes a method for obtaining a model of a single or a set of nonlinear loads (NLL) connected to a certain point of an electrical network. The basic assumption is that the network supplying the NLL has significant series impedances and is disturbed by other parallel, random, and unknown neighbor loads, sharing part of the supply system with the NLL. The main interest for obtaining the model is its further use to predict the amount and flow of harmonic currents generated by the NLL, in the case of adding a filter to reduce the harmonics distortion. The modeling technique used in the paper is based on multivariate multiple outputs regression and leads to a set of equations giving the NLL behavior (one for each of the harmonic currents). The model is obtained from data taken at measuring point and is only valid to predict the NLL behavior when new loads are connected at this point. The modeling method was first tested with V, I data coming from simulations using a MATLAB-Simulink SimPowerSystems toolbox. Finally, the method has been validated using V, I data taken in a real installation with different neighbor loads and under different load conditions.