- ポアッソン過程、そこからの逸脱を
- 「映画」的に時間を使ってみせたり
- 空間に配置してみせたりしたい
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){
evPow<-ep
events2<-events^evPow
cumeve<-cumsum(events2)
dur<-20
cumeve<-cumeve/max(cumeve)*dur
currentTime<-Sys.time()
for(i in 1:length(cumeve)){
while(Sys.time()-currentTime<cumeve[i]){
}
plot(cumeve[1:i],rep(0,i),pch=19,cex=0.5,xlim=c(0,max(cumeve)))
}
}
for(ep in evPows){
evPow<-ep
events2<-events^evPow
cumeve<-cumsum(events2)
dur<-10
cumeve<-cumeve/max(cumeve)*dur
currentTime<-Sys.time()
xlim<-ylim<-c(-1,1)
for(i in 1:length(cumeve)){
plot(c(10),c(10),pch=19,cex=4,xlim=xlim,ylim=ylim)
while(Sys.time()-currentTime<cumeve[i]){
}
plot(c(0),c(0),pch=19,cex=10,xlim=xlim,ylim=ylim)
}
}