## Not run:
# VB5 = dnc(Y,M=5,p=3,"VB")
# Gibbs5 = dnc(Y,M=5,p=3,"Gibbs")
# print(VB5)
# print(Gibbs5)
# BIC(VB5)
# BIC(Gibbs5)
# plot(VB5)
# plot(Gibbs5,aggregated=FALSE,plotRGL=FALSE)
# ## End(Not run)
Run the code above in your browser using DataLab