読者です 読者をやめる 読者になる 読者になる

INFINITY -数学とかプログラミングとか-

統計とプログラムを使って役に立たせたい

TeX用コマンド入力を支援するための辞書をご利用ください。
sanctuary's blogは,適当なことが書いてあります。

平面上の一様分布の直交変換

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)