平面上の一様分布の直交変換
theta=3.14/6 mm <- matrix(c(cos(theta),-sin(theta), sin(theta),cos(theta)),2,2) mm=t(mm) n=200 x=runif(n) y=runif(n) xt=matrix(0,ncol=n,nrow=2) for(i in 1:n){ for(p in 1:2){ for(q in 1:2){ xt[,i]=mm%*%c(x[i],y[i]) } } } plot(x,y,xlim=c(-1.5,1.5),ylim=c(-1.5,1.5)) par(new=T) plot(xt[1,],xt[2,],xlim=c(-1.5,1.5),ylim=c(-1.5,1.5),col=2)