NJ樹構(gòu)建(splitstree)

cd /mnt/d/peng/climate/paper-climate/climate-1/NJtree
plink --file 72breeds --indep-pairwise 100 25 0.05 --chr-set 29 --out goatld0.2
plink --file 72breeds --indep-pairwise 100 50 0.05 --chr-set 29 --out goatld50
plink --file 72breeds --extract goatld50.prune.in --chr-set 29 --recode --out goatforpca

用spider 軟件map,ped轉(zhuǎn)換成Arlequin格式

計(jì)算距離,轉(zhuǎn)換成矩陣,在splitstree中構(gòu)建進(jìn)化樹
格式為:
3
ABR 0 0.08517 0.0922
ANK 0.08517 0 0.07226
ARG 0.0922 0.07226 0

?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
【社區(qū)內(nèi)容提示】社區(qū)部分內(nèi)容疑似由AI輔助生成,瀏覽時(shí)請結(jié)合常識與多方信息審慎甄別。
平臺聲明:文章內(nèi)容(如有圖片或視頻亦包括在內(nèi))由作者上傳并發(fā)布,文章內(nèi)容僅代表作者本人觀點(diǎn),簡書系信息發(fā)布平臺,僅提供信息存儲(chǔ)服務(wù)。

相關(guān)閱讀更多精彩內(nèi)容

友情鏈接更多精彩內(nèi)容