We introduce the concepts of complex polyanalytic weighted Bergman spaces and of quaternionic polyanalytic weighted Bergman spaces of first and second kind. In these spaces we then prove qualitative and quantitative results in approximation by polyanalytic polynomials. The quantitative approximation results are given in terms of higher order $L^{p}$-moduli of smoothness and in terms of the best approximation quantity.