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