わざとボールを投げているかどうかにどれくらい早く気づけるか Npt<-500 mx<-0.2 my<-0 sdx<-1/5 sdy<-1/5 r<-1 # 2 は正規分布 # 正規分布が「自然」と感じるか？ # もっと距離での減衰の緩い分布(距離の１乗に比例して確率が減じるとか：指数分布)の方が自…

円形の的 xlim<-c(-1,1) ylim<-xlim ori<-c(0,0) plot(ori[1],ori[2],cex=60,pch=19,col=1,xlim=xlim,ylim=ylim) par(new=TRUE) plot(ori[1],ori[2],cex=50,pch=19,col=2,xlim=xlim,ylim=ylim) par(new=TRUE) plot(ori[1],ori[2],cex=40,pch=19,col=3,xlim=x…

確率事象をビジュアルに説明する企画(こちら)のためのメモ k<-24 col<-rainbow(k) dur<-10 # 継続時間 rateParam<-0.01 # 1メモリ分の時間 currentTime<-Sys.time() while(Sys.time()-currentTime < dur){ col<-c(col[2:k],col[1]) pie(rep(1, k), col = col…

ポアッソン過程、そこからの逸脱を 「映画」的に時間を使ってみせたり 空間に配置してみせたりしたい # ぽつりぽつりと起きる現象 m<-0.5 Ne<-20 events<-rexp(Ne*(1-proportionRunif-proportionRnorm),m) evPows<-c(0.01,0.1,0.5,1,2,3,5) for(ep in evPows…