ryamada22 https://blog.hatena.ne.jp/ryamada22/ ryamadaの遺伝学・遺伝統計学メモ https://ryamada22.hatenablog.jp/ R 遺伝的浮動 不動 有限個の状態の離散的世代経過を状態間遷移関数でシミュレートすれば、遺伝的ドリフトを試せる。 #染色体N本中p本が着目アレル #それがk倍になった後に、N本を抜き出し #そこにi本の着目アレルがある確率 probDrift<-function(N,p,k,i){ kN<-k*N kp<-k*p kNp<-kN-kp k1N<-(k-1)*N n11<-i n12<-kp-n11 n21<-N-n11 n22<-k1N-n12 return(exp(lgamma(kp+1)+lgamma(kNp+1)+lgamma(k1N+1)+lgamma(N+1)-lgamma(kN+1)-lgamma(n1… 190 <iframe src="https://hatenablog-parts.com/embed?url=https%3A%2F%2Fryamada22.hatenablog.jp%2Fentry%2F20100119%2F1263893351" title="状態遷移 - ryamadaの遺伝学・遺伝統計学メモ" class="embed-card embed-blogcard" scrolling="no" frameborder="0" style="display: block; width: 100%; height: 190px; max-width: 500px; margin: 10px 0px;"></iframe> Hatena Blog https://hatena.blog 2010-01-19 18:29:11 rich https://ryamada22.hatenablog.jp/entry/20100119/1263893351 1.0 100%