diff --git a/inst/stan/glm_multi_beta_binomial.stan b/inst/stan/glm_multi_beta_binomial.stan
index badc1d8..b882559 100755
--- a/inst/stan/glm_multi_beta_binomial.stan
+++ b/inst/stan/glm_multi_beta_binomial.stan
@@ -757,11 +757,11 @@ model{
   
   // Random intercept
   if(is_random_effect>0){
-    for(m in 1:(M-1)) random_effect_raw[,m] ~ std_normal();
+    for(m in 1:(M-1)) random_effect_raw[,m] ~ normal ( prior_mean_coefficients[1], prior_mean_coefficients[2]); 
     for(m in 1:(M-1)) random_effect_sigma_raw[m] ~ std_normal();
     for(m in 1:(M-1)) sigma_correlation_factor[m] ~ lkj_corr_cholesky(2);   // LKJ prior for the correlation matrix
     
-    for(m in 1:(M-1)) random_effect_raw_2[,m] ~ std_normal();
+    for(m in 1:(M-1)) random_effect_raw_2[,m] ~ normal ( prior_mean_coefficients[1], prior_mean_coefficients[2]);
     for(m in 1:(M-1)) random_effect_sigma_raw_2[m] ~ std_normal();
     if(is_random_effect>1) for(m in 1:(M-1)) sigma_correlation_factor_2[m] ~ lkj_corr_cholesky(2);   // LKJ prior for the correlation matrix