Modeling return to education in heterogeneous populations. An application to Italy