library(deldir) set.seed(115) x1 <- matrix(rep(array(1:20),10)) x2 <- matrix(rep(array(1:20)+0.5,10)) x <- rbind(x1, x2) y1 <- matrix(rep(array(1:10),each=20)) y2 <- matrix(rep(array(1:10)+0.5,each=20)) y <- rbind(y1*sqrt(3), y2*sqrt(3)) for (i in 1:400){ a=runif(1,0,2*pi) b=runif(1,-1,1) x[i]=x[i]+0.1*cos(a)*b y[i]=y[i]+0.1*sin(a)*b } xy <- data.frame(x,y) dxy <- deldir(xy) w <- tile.list(dxy) plot(w,number=TRUE, col.num="red", cex=0.3) law <- lawSummary(dxy) print(law)