メモ

チーム対決ソース(ポアソン分布で近似する)

# チーム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)