Continuum modeling of active nematics via data-driven equation discovery