The investigation of functional neuronal synchronization has recently become a growing field of research. With high temporal resolution, electroencephalography and magnetoencephalography are well-suited measurement techniques to identify networks of interacting sources underlying the recorded data. The analysis of the data in terms of effective connectivity, nevertheless, contains intrinsic issues such as the problem of volume conduction and the non-uniqueness of the inverse solution. Here, we briefly introduce a series of existing methods assessing these problems. To determine the locations of interacting brain sources robust to volume conduction, all computations are solely based on the imaginary part of the cross-spectrum as a trustworthy source of information. Furthermore, we demonstrate the feasibility of estimating causal relationships of systems of neuronal sources with the phase slope index in realistically simulated data. Finally, advantages and drawbacks of the applied methodology are highlighted and discussed.