1. Abstract Signals of brain electric neuronal activity, either invasively measured or non-invasively estimated, are commonly used for connectivity inference. One popular methodology assumes that the neural dynamics follow a multivariate autoregression, where the autoregressive coefficients represent the couplings among regions. If observation noise is present and ignored, as is common in practice, the estimated couplings are biased, affecting all forms of Granger-causality inference, both in time and in frequency domains. Significant nonsense coupling, i.e., nonsense connectivity, can appear when in reality there is none, since there is always observation noise in two possible forms: measurement noise, and activity from other brain regions due to volume conduction and low spatial resolution. This problem is critical, and is currently not being addressed, calling into question the validity of many Granger-causality reports in the literature. An estimation method that accounts for noise is based on an overdetermined system of high-order multivariate Yule-Walker equations, which give reduced variance estimators for the coupling coefficients of the unobserved signals. Simulation-based comparisons to other published methods are presented, demonstrating its adequate performance. In addition, simulation results are presented for a zero connectivity case with noisy observations, where the new method correctly reports no connectivity while classical analyses (as found in most software packages) report nonsense connectivity. For the sake of reproducible research, the supplementary material includes, in human readable format, all the time series data used here.