SOCR ≫ | TCIU Website ≫ | TCIU GitHub ≫ |
set.seed(1234)
x <- rvonmises(n=1000, mu=circular(pi/5), kappa=3)
y <- rvonmises(n=1000, mu=circular(-pi/3), kappa=5)
z <- rvonmises(n=1000, mu=circular(0), kappa=10)
resx <- density(x, bw=25)
res <- plot(resx, points.plot=TRUE,
xlim=c(-1.5,1), ylim=c(-1.1, 1.5), offset=1.1, shrink=1.2, lwd=3)
resy <- density(y, bw=25)
lines(resy, points.plot=TRUE, col=2, points.col=2, plot.info=res, offset=1.1, shrink=1.45, lwd=3)
resz <- density(z, bw=25)
lines(resz, points.plot=TRUE, col=3, points.col=3, plot.info=res, offset=1.1, shrink=1.2, lwd=3)
## quartz_off_screen
## 2