library(fAsianOptions) cgamma(1i) my.cbeta <- function(a,b){ ret <- cgamma(a)*cgamma(b)/cgamma(a+b) return(ret) } my.dcbeta <- function(x,a,b){ x^(a-1) * (1-x)^(b-1) / my.cbeta(a,b) } x <- seq(from=-1,to=2,length=100) a <- -0.5 b <- -0.5 y …

今、N回の観測をして、○がn1,×がn2回だったとする この後、m回繰り返したら、m回のうち、何回成功するんだろう？ まず観測データから、○確率pの確率密度関数をベータ関数で表して、ただし、はベータ関数 今、○確率がpのときにm回のうちk回成功するのはなので…