PRML 2章 ベータ分布のグラフを描く

Rに慣れるため本にあるグラフをひたすら描いてみる。

png("beta01.png",height=600,width=800)

fbeta <- function(x) { dbeta(x,5,5) } 
curve(fbeta,0,1,col="1")
fbeta <- function(x) { dbeta(x,2,2) } 
curve(fbeta,0,1,col="2",add=T)
fbeta <- function(x) { dbeta(x,1,1) } 
curve(fbeta,0,1,col="3",add=T)
fbeta <- function(x) { dbeta(x,0.5,0.5) } 
curve(fbeta,0,1,col="4",add=T)
fbeta <- function(x) { dbeta(x,0.1,0.1) } 
curve(fbeta,0,1,col="5",add=T)

legend(0,2.5,c("a=5,b=5","a=2,b=2","a=1,b=1","a=0.5,b=0.5","a=0.1,b=0.1"), cex=0.6,col=1:5,lty=rep(1,5))
dev.off()

某所からの単なるパクリだが気にしてはいけない。

hatenablog markdown モードで書いてるけど数式の入れ方が分からない。Beta函数の数式を入れたかった。