連鎖平衡を代数統計する

• 2アレル多型が２つある。そのアレル頻度をそれぞれp1,p2,q1,q2とする
• 連鎖平衡の下での４ハプロタイプの頻度を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*q1,p2*q2;
ideal i0=0;
setring h;
preimage(p,f,i0);

• • が出ます