We propose a method to compute provably invariant subsets of the region-of-attraction for the asymptotically stable equilibrium points of uncertain nonlinear dynamical systems. We consider polynomial dynamics with perturbations that either obey local polynomial bounds or are described by uncertain parameters multiplying polynomial terms in the vector field. This uncertainty description is motivated by both incapabilities in modeling, as well as bilinearity and dimension of the sum-of-squares programming problems whose solutions provide invariant subsets of the region-of-attraction. We demonstrate the method on three examples from the literature and a controlled short period aircraft dynamics example.