We present preliminary results of a search for previously unidentified stellar associations in Gaia DR2. We have modified a Python code originally intended to identify wide stellar binaries through a Bayesian formulation, which includes correlated uncertainties in the proper motions and parallaxes, to expand it to larger stellar groups of 10 or more stars. The search volume includes the whole sky within 500 pc, where stellar parallaxes are more discernible. The code has been optimized by testing over several well known stellar associations, including the Pleiades open cluster and the Scorpius-Centaurus Association. Here we discuss the results of our search for new stellar associations within Gaia DR2.