これも書き換える予定 いろいろ試したんだけど、スプライン補間がうまくできなくて、値が微妙に異なる感じになってしまった... ちょっと後回しなう。 > beta.select function (quantile1, quantile2) { betaprior1 = function(K, x, p) { m.lo = 0 m.hi = 1 …

Rのpdisc()関数をpythonで書いた。 参考"yuji_nkのメモ" "ハリ・セルダンになりたくて" from numpy import * def pdisc(p,prior, data): p1=[] s=data[0] f=data[1] p1 = p + 0.5*(p==0) - 0.5*(p==1) p1=array(p1) like= s * log(p1) + f * log(1.0-p1) for…