連鎖平衡

連鎖平衡を代数統計する

2アレル多型が2つある。そのアレル頻度をそれぞれp1,p2,q1,q2とする 連鎖平衡の下での4ハプロタイプの頻度をh1,h2,h3,h4とする Singularを使ってやってみる ring p=0,(p1,p2,q1,q2),lp; ring h=0,(h1,h2,h3,h4),lp; setring p; map f = h,p1*q1,p1*q2,p2*q…