## Not run:
# varobj <- szbsvar(Y, p, z = NULL, lambda0, lambda1, lambda3, lambda4,
# lambda5, mu5, mu6, ident, qm = 4)
# A0.posterior <- gibbs.A0(varobj, N1, N2)
# fit <- posterior.fit(varobj, A0.posterior)
# print(fit)
# ## End(Not run)
Run the code above in your browser using DataLab