尤度

二項分布とベータ分布

Rでやれば 確率pで起きる事象をn=a+b回実施して、a回起きる確率は dbinom(a,n,p) #もしくは choose(n,a)*p^a*(1-p)^b 他方、n回試行して、a回成功したときの、成功確率をpとみなしたときの尤度は dbeta(p,a+1,b+1)/(n+1) であって、等しい 尤度なので、の範…

ソースを書いてもらうための説明文

2SNP間の連鎖不平衡の評価値LODを算出するためのプログラムソースを書いてもらうための説明文 以前の記事(別テイスト)はこちら