有限個の状態の離散的世代経過を状態間遷移関数でシミュレートすれば、遺伝的ドリフトを試せる。 #染色体N本中p本が着目アレル #それがk倍になった後に、N本を抜き出し #そこにi本の着目アレルがある確率 probDrift<-function(N,p,k,i){ kN<-k*N kp<-k*p kNp…