収束の様子をきれいに描く

  • ポアソン分布を気にしている(こちら)
  • こちらポアソン事象の観測データを積み重ねていき、そのポアソン分布のパラメタが収束していく様子をガンマ関数で視覚的に見せている
  • RのRColorBrewerパッケージというのを用いて色を作っている
library(RColorBrewer)
## create a sequential palette for usage and show colors
mypalette<-brewer.pal(7,"Greens")
image(1:7,1,as.matrix(1:7),col=mypalette,xlab="Greens (sequential)",
       ylab="",xaxt="n",yaxt="n",bty="n")
# k は事象の生起回数の総和、Nは観測の機会数
dgamma(lambda.val, k+1,N+1)
  • 共役事前分布のダイアグラム(こちら)