The concept of slice regular function over the real algebra H of quaternions is a generalization of the notion of holomorphic function of a complex variable. Let Ω be an open subset of H, which intersects R and is invariant under rotations of H around R. A function f : Ω → H is slice regular if it is of class C 1 and, for all complex planes C I spanned by 1 and a quaternionic imaginary unit I, the restriction f I of f to Ω I = Ω ∩ C I satisfies the Cauchy-Riemann equations associated to I, i.e., ∂ I f I = 0 on Ω I , whereWe define global slice polyanalytic functions of order n as the functions f : Ω → H, which admit a decomposition of the form f (x) = n−1 h=0 x h f h (x) for some slice regular functions f 0 , . . . , f n−1 . Global slice polyanalytic functions of any order n are slice polyanalytic of the same order n. The converse is not true: for each n ≥ 2, we give examples of slice polyanalytic functions of order n, which are not global.The aim of this paper is to study the continuity and the differential regularity of slice regular and global slice polyanalytic functions viewed as solutions of the slice-by-slice differential equations ∂ n I f I = 0 on Ω I and as solutions of their global version ϑ n f = 0 on Ω \ R.Our quaternionic results extend to the monogenic case.