[论文解读] Phylogenetic ideals and varieties for the general Markov model

Elizabeth S. Allman, John A. Rhodes|ArXiv.org|Oct 28, 2004

Advanced Combinatorial Mathematics参考文献 18被引用 12

一句话总结

本文證明了Pachter-Sturmfels的猜想：對於任意二叉系統發育樹上的兩狀態一般馬可夫模型，系統發育不變量的完整素理想由沿每條內部邊的聯合機率張量的二維展平（flattening）之3×3子式生成。關鍵方法在於利用三物種多樣體填滿其周圍空間（ℙ⁷）的事實，從而能使用矩陣作用與展平運算進行歸納論證，最終證明僅靠邊基不變量即可定義整個多樣體。

ABSTRACT

The general Markov model of the evolution of biological sequences along a tree leads to a parameterization of an algebraic variety. Understanding this variety and the polynomials, called phylogenetic invariants, which vanish on it, is a problem within the broader area of Algebraic Statistics. For an arbitrary trivalent tree, we determine the full ideal of invariants for the 2-state model, establishing a conjecture of Pachter-Sturmfels. For the $κ$-state model, we reduce the problem of determining a defining set of polynomials to that of determining a defining set for a 3-leaved tree. Along the way, we prove several new cases of a conjecture of Garcia-Stillman-Sturmfels on certain statistical models on star trees, and reduce their conjecture to a family of subcases.

研究动机与目标

• 確定任意二叉樹上一般馬可夫模型的完整素理想系統發育不變量。
• 解決Pachter與Sturmfels關於兩狀態模型的第13號猜想，該猜想主張沿邊展平的3×3子式生成完整的不變理想。
• 將κ狀態模型尋找不變量的問題簡化至三葉樹的情形，從而簡化一般分類問題。
• 證明Garcia-Stillman-Sturmfels關於星形樹統計模型的猜想之新情形，並將其簡化為有限多組子情形。

提出的方法

• 利用κ=2時，三物種多樣體V(T₃)填滿其周圍空間ℙ⁷的事實，此一特殊代數結構在κ>2時不成立。
• 對物種數量使用歸納法，將樹Tₙ分解為Tₙ₋₁*T₃，其中aₙ₋₁與aₙ形成一對櫻桃（cherry）。
• 對聯合機率張量P應用矩陣展平運算，將其視為以內部邊兩側的物種分組後的矩陣。
• 利用GL(4,ℂ)在櫻桃結構上的作用，使Lemma 12得以應用，進而連結不同大小樹的不變理想。
• 將理想成員資格問題簡化為證明形如f(Flatₑ(P)*B)的多項式屬於由3×3子式生成的理想，並利用展平運算的秩條件。
• 利用多樣體CV(Tₙ)在矩陣乘法下的封閉性，證明理想由邊基不變量生成。

实验结果

研究问题

• RQ1 沿每條內部邊的二維展平中，3×3子式是否生成任意二叉n物種樹上兩狀態一般馬可夫模型的完整素理想系統發育不變量？
• RQ2 能否將κ狀態模型不變理想之確定問題簡化至三葉樹的情形？
• RQ3 由三物種多樣體V(T₃)的二維展平所產生的不變量是否如Pachter與Sturmfels所猜想的那樣，生成兩狀態模型的完整理想？
• RQ4 當κ>2時，一般馬可夫模型的系統發育多樣體具有何種結構？與κ=2的情形有何差異？
• RQ5 Garcia-Stillman-Sturmfels關於星形樹模型的猜想能否在特定子情形下證明？若可，是哪些子情形？

主要发现

• 任意二叉n物種樹上兩狀態一般馬可夫模型的完整素理想系統發育不變量，由沿每條內部邊的聯合機率張量二維展平之3×3子式生成。
• 當κ=2時，三物種樹V(T₃)填滿其周圍空間ℙ⁷，此為關鍵結構性質，使歸納論證得以成立。
• 對於任意二叉樹Tₙ，多樣體CV(Tₙ)是CV(Tₙ₋₁) × M(4,ℂ)在矩陣乘法下像的閉包，此性質使證明中的歸納步驟得以實現。
• κ狀態模型的不變理想可簡化為對三葉樹不變量的研究，大幅簡化分類問題。
• 本文證明了Garcia-Stillman-Sturmfels猜想在星形樹模型上的新情形，並將完整猜想簡化為有限多組子情形。
• 該方法確立所有理想中的不變量均可表示為邊基3×3子式的組合，從而確認Pachter與Sturmfels對兩狀態模型的猜想。