メモ

個々のＳＮＰのアレル頻度は、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)