data(friendship)
set.seed(123)
dncObj <- dnc(friendship,M=4,p=3,method="Gibbs",
controls=list(nDraws=250,burnin=50,
MaxItStg2=25,epsilonStg2=1e-15))
par(mar=rep(0,4)+0.05)
plot(dncObj,plotRGL=FALSE,pch=16,phi=60,lwd=2,cex=1.5)
Run the code above in your browser using DataLab