メモ
チーム対決ソース(ポアソン分布で近似する)
# チーム2の得点を動かして # 勝敗アイコ率をとる a1<-15 N<-1000 a2s<-seq(from=0,to=30,length.out=N) win1<-win2<-tie<-rep(0,N) for(i in 1:N){ a2<-a2s[i] p1<-dpois(0:100,a1) p2<-dpois(0:100,a2) ylim<-c(0,max(p1,p2)) plot(p1,ylim=ylim,type="b") par(new=TRUE) plot(p2,ylim=ylim,col=2,type="b") # (1:n)x(1:n)の同時確率は p12<-(p1)%*%t(p2) win1[i]<-sum(p12[lower.tri(p12)]) win2[i]<-sum(p12[upper.tri(p12)]) tie[i]<-sum(diag(p12)) } ylim<-c(0,max(win1,win2,tie)) plot(a2s,win1,ylim=ylim) par(new=TRUE) plot(a2s,win2,ylim=ylim,col=2) par(new=TRUE) plot(a2s,tie,ylim=ylim,col=3)