set.seed(1234567)
ccoef<-axialnntsrandominitial(M=4)
data<-axialnntssimulation(100,ccoef,M=4)$simulations
hist(data,breaks=seq(0,pi,pi/7),xlab="Orientations (radians)",freq=FALSE,
ylab="",main="",ylim=c(0,1))
axialnntsplot(ccoef,M=4,add=TRUE)
Run the code above in your browser using DataLab