(9) ２次元の的当て(２次元正規分布)

• 対立仮説が１点

# ２次元の的当て
Npt<-50
xlim<-c(-1,1)
ylim<-xlim
ori<-c(0,0)

mx<-0.2
my<-0
sdx<-1/5
sdy<-1/5
r<-1 # 2 は正規分布
rateParam<-0.1
XYs<-matrix(10,Npt,2)
for(i in 1:Npt){
	plot(ori[1],ori[2],cex=60,pch=19,col=1,xlim=xlim,ylim=ylim,main="",xlab="",ylab="")
par(new=TRUE)
plot(ori[1],ori[2],cex=50,pch=19,col=2,xlim=xlim,ylim=ylim,main="",xlab="",ylab="")
par(new=TRUE)
plot(ori[1],ori[2],cex=40,pch=19,col=3,xlim=xlim,ylim=ylim,main="",xlab="",ylab="")
par(new=TRUE)
plot(ori[1],ori[2],cex=30,pch=19,col=4,xlim=xlim,ylim=ylim,main="",xlab="",ylab="")
par(new=TRUE)
plot(ori[1],ori[2],cex=20,pch=19,col=5,xlim=xlim,ylim=ylim,main="",xlab="",ylab="")
par(new=TRUE)
plot(ori[1],ori[2],cex=10,pch=19,col=6,xlim=xlim,ylim=ylim,main="",xlab="",ylab="")
par(new=TRUE)

	par(new=TRUE)
	tmpx<-rnorm(1,mean=mx,sd=sdx)
	tmpy<-rnorm(1,mean=my,sd=sdy)
	tmpR<-sqrt(tmpx^2+tmpy^2)
	XYs[i,1]<-tmpx*tmpR^(r/2)
	XYs[i,2]<-tmpy*tmpR^(r/2)
	plot(XYs[,1],XYs[,2],cex=1,pch=19,col=gray(100/100),xlim=xlim,ylim=ylim,main="",mex=30,xlab="",ylab="")
	n0<-i%%10
	n1<-((i-n0)/10)%%10
	n2<-((i-n0-n1*10)/100)%%10
	par(new=TRUE)
	plot(0.9,0.9,pch=n0+48,cex=4.5,xlim=xlim,ylim=ylim,main="")
	par(new=TRUE)
	plot(0.77,0.9,pch=n1+48,cex=4.5,xlim=xlim,ylim=ylim,main="")
	par(new=TRUE)
	plot(0.64,0.9,pch=n2+48,cex=4.5,xlim=xlim,ylim=ylim,main="")
	tmpCurrent<-Sys.time()
	while(Sys.time()-tmpCurrent<rateParam){
		
	}
}

• 中心からある距離の点(ならどこでもよい)を狙っている

# 半径が一定Rtの点を的にする
Npt<-500
xlim<-c(-1,1)
ylim<-xlim
ori<-c(0,0)

mx<-0.2
my<-0
sdx<-1/16
sdy<-1/16

Rts<-seq(from=0,to=0.6,by=0.01)

for(j in 1:length(Rts)){
	Rt<-Rts[j]
r<-1 # 2 は正規分布
rateParam<-0.1
XYs<-matrix(10,Npt,2)
#for(i in 1:Npt){
i<-Npt
	plot(ori[1],ori[2],cex=60,pch=19,col=1,xlim=xlim,ylim=ylim,main="",xlab="",ylab="")
par(new=TRUE)
plot(ori[1],ori[2],cex=50,pch=19,col=2,xlim=xlim,ylim=ylim,main="",xlab="",ylab="")
par(new=TRUE)
plot(ori[1],ori[2],cex=40,pch=19,col=3,xlim=xlim,ylim=ylim,main="",xlab="",ylab="")
par(new=TRUE)
plot(ori[1],ori[2],cex=30,pch=19,col=4,xlim=xlim,ylim=ylim,main="",xlab="",ylab="")
par(new=TRUE)
plot(ori[1],ori[2],cex=20,pch=19,col=5,xlim=xlim,ylim=ylim,main="",xlab="",ylab="")
par(new=TRUE)
plot(ori[1],ori[2],cex=10,pch=19,col=6,xlim=xlim,ylim=ylim,main="",xlab="",ylab="")
par(new=TRUE)

	par(new=TRUE)
	for(i in 1:Npt){
	tmpr<-runif(1)*2*pi
	mtx2<-Rt*cos(tmpr)
	mty2<-Rt*sin(tmpr)

	tmpx<-rnorm(1,mean=mtx2,sd=sdx)
	tmpy<-rnorm(1,mean=mty2,sd=sdy)
	tmpR<-sqrt(tmpx^2+tmpy^2)
	XYs[i,1]<-tmpx*tmpR^(r/2)
	XYs[i,2]<-tmpy*tmpR^(r/2)
	}
	plot(XYs[,1],XYs[,2],cex=1,pch=19,col=gray(100/100),xlim=xlim,ylim=ylim,main="",mex=30,xlab="",ylab="")
	#n0<-i%%10
	#n1<-((i-n0)/10)%%10
	#n2<-((i-n0-n1*10)/100)%%10
	#par(new=TRUE)
	#plot(0.9,0.9,pch=n0+48,cex=4.5,xlim=xlim,ylim=ylim,main="")
	#par(new=TRUE)
	#plot(0.77,0.9,pch=n1+48,cex=4.5,xlim=xlim,ylim=ylim,main="")
	#par(new=TRUE)
	#plot(0.64,0.9,pch=n2+48,cex=4.5,xlim=xlim,ylim=ylim,main="")
	#tmpCurrent<-Sys.time()
	#while(Sys.time()-tmpCurrent<rateParam){
		
	#}
#}
}