個々のＳＮＰのアレル頻度は、0-1に均一にした上で、ハプロタイプ頻度をそれなりに出したい N<-1000 k<-10 c<-1 xxxx<-rdirichlet(N,rep(1,2^k)*1/(2^(k-c))) yyyy<-xxxx[,1:2^(k-c)] zzzz<-apply(yyyy,1,sum) hist(zzzz)

ディリクレ分布については、Wikipediaの記事がこちら。 このWikipedia記事の図は、３つの独立事象について、 をの場合について、の範囲での範囲について計算し、その結果を、x1,x2の軸についてプロットしたものである。 これをエクセルで計算するとこんな感…

bayesmパッケージをダウンロード library(bayesm) help(rdirichlet) set.seed(66) rdirichlet(c(rep(3,5))) sum(rdirichlet(c(rep(3,5)))) [1] 1

作業は統計環境『Ｒ』を用いるものとする。ガンマ分布や、Dirichlet分布からの乱数発生関数 rdirichlet (MCMCpack)を用いる(その使い方についてはこちら) サイコロ 今、サイコロの「できのよさ」を調べながら、そのサイコロで博打をうつことを想定する サイ…