We consider a class of SO(10) models with flavor symmetries in the Yukawa sector and investigate their viability by performing numerical fits to the fermion masses and mixing parameters. The fitting procedure involves a top-down approach in which we solve the renormalization group equations from the scale of grand unification down to the electroweak scale. This allows the intermediate scale right-handed neutrinos and scalar triplet, involved in the type I and II seesaw mechanisms, to be integrated out at their corresponding mass scales, leading to a correct renormalization group running. The result is that, of the 14 models considered, only two are able to fit the known data well. Both these two models correspond to Z(2) symmetries. In addition to being able to fit the fermion masses and mixing parameters, they provide predictions for the sum of light neutrino masses and the effective neutrinoless double beta decay mass parameter, which are both within current observational bounds.