Combining data assimilation and machine learning to infer unresolved scale parametrisation