PROC MCMC DATA=ndata NBI=5000 NMC=40000 THIN=2 SEED=17 INIT=RANDOM DIAG=(ESS GEWEKE) STATISTICS(ALPHA=0.05)=(SUMMARY INTERVAL) OUTPOST=histnorm; ODS SELECT PARAMETERS POSTSUMMARIES POSTINTERVALS GEWEKE ESS TADPANEL; ARRAY b[2] L S; ARRAY beta[2]; ARRAY Sigma_b[2,2]; ARRAY beta0[2] (0 0); ARRAY sigma0[2,2] (1000 0 0 1000); ARRAY V[2,2] (1 0 0 1); PARMS beta {5 2}; PARMS sigma_b {1 0 0 1}; PARMS var_e 1; PRIOR beta ~ MVN(beta0, sigma0); PRIOR Sigma_b ~ IWISH(2, V); PRIOR var_e ~ IGAMMA(0.001, SCALE=0.001); RANDOM b ~ MVN(beta, Sigma_b) SUBJECT=id; mu = L + S * time; MODEL y ~ NORMAL(mu, var=var_e); RUN;