The inverse problem method is tested for a class of monomer-dimer statistical mechanics models that contain also an attractive potential and display a mean-field critical point at a boundary of a coexistence line. The inversion is obtained by analytically identifying the parameters in terms of the correlation functions and via the maximum-likelihood method. The precision is tested in the whole phase space and, when close to the coexistence line, the algorithm is used together with a clustering method to take care of the underlying possible ambiguity of the inversion. 2 Definition of the model Let G = (V, E) be a finite simple graph with vertex set V and edge set E = {uv ≡ {u, v}|u = v ∈ V }. Definition 2.1. A dimer configuration D on the graph G is a set of dimers 7