In previous works, an approach to the study of cyclic functions in reproducing kernel Hilbert spaces has been presented, based on the study of so called optimal polynomial approximants. In the present article, we extend such approach to the (non-Hilbert) case of spaces of analytic functions whose Taylor coefficients are in β p (Ο), for some weight Ο. When Ο = {(k + 1) Ξ± } kβN , for a fixed Ξ± β R, we derive a characterization of the cyclicity of polynomial functions and, when 1 < p < β, we obtain sharp rates of convergence of the optimal norms.