的外れ

  • わざとボールを投げているかどうかにどれくらい早く気づけるか
Npt<-500

mx<-0.2
my<-0
sdx<-1/5
sdy<-1/5
r<-1 # 2 は正規分布
# 正規分布が「自然」と感じるか?
# もっと距離での減衰の緩い分布(距離の1乗に比例して確率が減じるとか:指数分布)の方が自然?→ノーコン・ピッチャー??
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){
		
	}
}