The cosmic distance duality relation (CDDR), D L (1+z) −2 /D A = η = 1, with D L and D A , being the luminosity and angular diameter distances, respectively, is a crucial premise in cosmological scenarios. Many investigations try to test CDDR through observational approaches, even some of these also consider a deformed CDDR, i.e., η = η(z). In this paper, we use type Ia supernovae luminosity distances and galaxy cluster measurements (their angular diameter distances and gas mass fractions) in order to perform a Bayesian model comparison between η(z) functions. We show that the data here used is unable to pinpoint, with a high degree of Bayesian evidence, which η(z) function best captures the evolution of CDDR.