A stable and robust method to identify modules of functionally coherent genes