Complex models typically involve intractable likelihood functions which, from a Bayesian perspective, lead to intractable posterior distributions. In this context, Approximate Bayesian computation (ABC) methods can be used in order to obtain a valid posterior approximation. However, when simulation from the model is computationally demanding, then the ABC approach may be cumbersome. We discuss an alternative method, where the intractable likelihood is approximated by a quasi-likelihood calculated through an algorithm that is reminiscent of the ABC. The proposed approximation method requires less computational effort than ABC. An extension to multiparameter models is also considered and the method is illustrated by several examples