We calculate the shear viscosity over entropy density ratio η/s in Yang-Mills theory from the Kubo formula using an exact diagrammatic representation in terms of full propagators and vertices using gluon spectral functions as external input. We provide an analytic fit formula for the temperature dependence of η/s over the whole temperature range from a glueball resonance gas at low temperatures, to a high-temperature regime consistent with perturbative results. Subsequently we provide a first estimate for η/s in QCD.PACS numbers: 12.38. Aw, 11.10.Wx , 11.15.Tk Introduction -The experimental heavy-ion programs at RHIC [1,2] and at the LHC [3] explore the physics of the quark-gluon plasma (QGP). It turns out that the dynamics of the hot plasma created in heavy-ion collisions is well-described by hydrodynamics. Therefore, the determination of transport coefficients in the QGP is of great interest. One aspect is that the inference of the initial state physics requires a precise description of the hydrodynamical evolution, which in turn depends on transport coefficients as microscopic input, [4]. In particular, the viscosity over entropy ratio η/s governs the efficiency of the conversion of the initial spatial anisotropy into a momentum anisotropy of the final state.For the determination of η/s and its temperature dependence in the quark-gluon plasma, theoretical approaches face several challenges.The temperature regimes below and above the critical temperature T c are characterised by different degrees of freedom, and for temperatures T 2T c non-perturbative effects become important. Of particular interest is the vicinity of T c , where the minimum for η/s is expected [5,6]. A universal lower bound for η/s of 1/4π was conjectured in [7] using the AdS/CFT correspondence. Indeed, measurements of the elliptic flow v 2 indicate a value for η/s which is of the order of this lower bound [8]. The bound has been tested theoretically with several methods for the QGP [9-15], but also for other potentially perfect liquids, such as ultracold atoms [16][17][18].The Kubo formulae relate η to the energy-momentum tensor (EMT) [19]. Spectral functions are real-time quantities and cannot be obtained directly from Euclidean correlation functions. However, the direct calculation of real-time correlation functions represents a notoriously difficult problem in non-perturbative approaches to quantum field theory. Even though first computations in this direction have been performed e.g. in [20,21], we shall utilise Euclidean correlation functions within a numerical analytic continuation.In this work we study the shear viscosity over entropy ratio η/s in pure SU (3) Landau gauge Yang-Mills (YM) theory within the approach set-up in [9]. In the present