R語言怎么分析群體遺傳進化進化樹

本篇內容主要講解“R語言怎么分析群體遺傳進化進化樹”，感興趣的朋友不妨來看看。本文介紹的方法操作簡單快捷，實用性強。下面就讓小編來帶大家學習“R語言怎么分析群體遺傳進化進化樹”吧!

群體遺傳進化進化樹分析方法：

cd $workdir  #回到工作目錄
mkdir 01.phylo_tree
cd 01.phylo_tree
#文件格式轉換
run_pipeline.pl  -Xmx5G -importGuess  $workdir/00.filter/clean.vcf.gz  \
    -ExportPlugin -saveAs supergene.phy -format Phylip_Inter
#最大似然法構建進化樹
#方法1：fasttree 構建進化樹
fasttree -nt -gtr  supergene.phy   >  fasttree.nwk
#方法2：iqtree 構建進化樹  設置boots值，最大似然法
iqtree2 -s supergene.phy -st DNA -T 2  -mem 8G \
    -m  GTR  -redo \
    -B 1000 -bnni \
    --prefix iqtree 
    
#方法3：raxml 構建進化樹 最大似然法
#raxml-ng  -msa supergene.phy --model GTR  --prefix raxml_tree \
    #    --threads 2 --seed 123   1>raxml.log 2>raxml.err
## with bootstrap
#raxml-ng -all  -msa supergene.phy --model GTR --bs-trees 1000 \
    #    --prefix raxml_tree_bootstrap --threads 2 --seed 123   1>raxml_bs.log 2>raxml_bs.err
#方法4：phylip NJ法構建進化樹  Phylip 格式對樣品ID字符要求不超過10個，不然會截斷
cp supergene.phy infile
echo -e "Y\n" |dnadist
cp infile.dist infile
echo -e "Y\n"  |neighbor
mv outtree outtree.nwk

到此，相信大家對“R語言怎么分析群體遺傳進化進化樹”有了更深的了解，不妨來實際操作一番吧！這里是億速云網站，更多相關內容可以進入相關頻道進行查詢，關注我們，繼續學習！