library(circular) LightOn=c() LightOff=c() plot.circular(x, stack=T,sep=0.08,bins=100) rayleigh.test(x) mean.circular(x) rho.circular(x) watson.two.test(x,y)