The multilevel models are used in applied sciences including social sciences, sociology, medicine, economic for analysing correlated data. There are various approaches to estimate the model parameters when the responses are normally distributed. To implement the Bayesian approach, a generalized version of the Markov Chain Monte Carlo algorithm, which has a simple structure and removes the correlations among the simulated samples for the fixed parameters and the errors in higher levels, is used in this article. Because the dimension of the covariance matrix for the new error vector is increased, based upon the Cholesky decomposition of the covariance matrix, two methods are proposed to speed the convergence of this approach. Then, the performances of these methods are evaluated in a simulation study and real life data.