特に、何か面白いことをしたわけではないですが、 ベータ関数を描いたので、コード載せとく。 logf = function(x){ rr = 0 for(i in 1:x){ rr = rr + log(i) } return(rr) } #beta distibutionの点の値を求める関数 d_beta = function(theta,n=129,y=118){ r…