距離法構(gòu)建進化樹演示文稿

距離法構(gòu)建進化樹演示文稿目前一頁\總數(shù)三十頁\編于十九點變形距離法如果假定速率恒定在各譜系間不成立則UPGMA法可能會給出錯誤的拓撲圖目前二頁\總數(shù)三十頁\編于十九點目前三頁\總數(shù)三十頁\編于十九點不過該拓撲圖錯誤也許能用被稱為變形距離法(transformeddistancemethod)的校正方法來訂正。簡單地說此法是用一個組外單位outgroup作參照物來對被研究譜系間進化速率不等的情況作一些校正，然后對新得到的距離矩陣應用UPGMA法從而推測出該樹的拓撲圖。目前四頁\總數(shù)三十頁\編于十九點組外單位是一個我們對其有外在認識的(OTU)。這些外在認識，如分類學的或古生物學的知識能清楚地表明該(OTU)已先于所有其他被研究的(OTU)而從其共同祖先分岐出來。在現(xiàn)在的例子中，我們假定分類單位對所有別的分類單位來說是一個組外單位，那么即可當作參照物，并用下式來變換距離。目前五頁\總數(shù)三十頁\編于十九點這里dij是變形距離i=A,B或C，d=(dAD+dBD+dCD)/3d項的引入是為了保證所有dij的值都是正的在我們的例子里d=37/3而新距離矩陣中關(guān)于分類單位A、B和C的值為所以A和B最先被聚類在一起,然后,C再加入到該樹中。目前六頁\總數(shù)三十頁\編于十九點Least-squaresmethod最小二乘法（LS）

包括一般LS法和加權(quán)LS法容許各個分支核苷酸替代速率有所不同算法：實踐中，可能的拓撲結(jié)構(gòu)非常多，逐個計算不可能。一般采用計算一個拓撲結(jié)構(gòu)，然后利用分支交換法搜索。目前七頁\總數(shù)三十頁\編于十九點在標準的LS法中，所有可能的拓樸結(jié)構(gòu)的相比后，最小的拓樸結(jié)構(gòu)即為最終的樹目前八頁\總數(shù)三十頁\編于十九點目前九頁\總數(shù)三十頁\編于十九點分支長度的估計目前十頁\總數(shù)三十頁\編于十九點目前十一頁\總數(shù)三十頁\編于十九點有時某一被估出的枝長可能會是負的由于實際枝長絕不可能是負的所以最好用來代替這樣的估值目前十二頁\總數(shù)三十頁\編于十九點目前十三頁\總數(shù)三十頁\編于十九點目前十四頁\總數(shù)三十頁\編于十九點目前十五頁\總數(shù)三十頁\編于十九點目前十六頁\總數(shù)三十頁\編于十九點最小進化（EM）法原理在本方法中，所有分支長度估計的和為：所有可能的拓撲結(jié)構(gòu)都要計算S值，具有最小S值的拓撲結(jié)構(gòu)被挑選作為最優(yōu)樹我們已提及分支長度的LS估計是由距離估計dij的函數(shù)給出的，即目前十七頁\總數(shù)三十頁\編于十九點目前十八頁\總數(shù)三十頁\編于十九點目前十九頁\總數(shù)三十頁\編于十九點Neighborjoiningmethod在一個無根兩分叉樹中,如果兩個OTU通過一個內(nèi)部節(jié)點聯(lián)結(jié),則它們就被說成是近鄰(neighbors)比如4個OTU，i、j、m、n。、假定第一個和的值最小，那么我們把i和j對以及m和n對都記1分,其余對記0分。在對所有可能的四單位組都評過分以后,得總分最高的對子即被選為第一個近鄰對目前二十頁\總數(shù)三十頁\編于十九點算法用NJ方法構(gòu)建一棵星狀樹，這就是假設不存在聚合群如果估計一棵星狀樹的分支長度并計算出所有分支的總和（S0），那么這個數(shù)目應該比真實樹或最終的樹的Sf要大目前二十一頁\總數(shù)三十頁\編于十九點目前二十二頁\總數(shù)三十頁\編于十九點目前二十三頁\總數(shù)三十頁\編于十九點目前二十四頁\總數(shù)三十頁\編于十九點目前二十五頁\總數(shù)三十頁\編于十九點目前二十六頁\總數(shù)三十頁\編于十九點目前二十七頁\總數(shù)三十頁\編于十九點Neighborjoiningmethod鄰接法的一般步驟如下1．采用了Jukes-Cautor單參數(shù)模型，計算出第i條和第j條序列距離dij計算公式：其中q為兩個序列中相應位置上相同堿基的概率目前二十八頁\總數(shù)三十頁\編于十九點Neighborjoiningmethod2．計算第i個葉子結(jié)點（即第i個序列）的凈分歧度ri其中N是葉子結(jié)點的個數(shù)，dik為葉子結(jié)點i和葉子結(jié)點k之間的距離目前二十九頁\總數(shù)三十頁\編于十九點Neighborjoiningmethod3．計算任意兩兩結(jié)點i和j之間的速率校正距離（rate-correcteddistance）Mi