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