Adhesively bonded joints are used in several industrial sectors. Cohezive Zone Modes can be used to predict the adhesive mechanical behaviour. This work presents an approach to calibrate Cohesive Zone Models (CZM) by means of Statistical Inverse Analysis. The Bayesian framework for Inverse Problems is used to infer about the CZM model parameters. The solution corresponds to the exploration of the posterior probability density function of the model parameters. The exploration of the posterior density is performed by means of Markov Chain Monte Carlo (MCMC) methods mixing Population-Based MCMC with Adaptive Metropolis (AD) strategies. The assessment of the approach is performed using measured data from a single-lap shear experimental set-up. Measured data from 5 test-specimens is used for calibration and measured data from five other test-specimens is used for model validation. It is proposed a stochastic effective model for the CZM parameters. The predictions of maximum force and maximum displacement that are provided by the effective model are in accordance with measured data that is used for validation.