中文字幕av专区_日韩电影在线播放_精品国产精品久久一区免费式_av在线免费观看网站

溫馨提示×

溫馨提示×

您好,登錄后才能下訂單哦!

密碼登錄×
登錄注冊×
其他方式登錄
點擊 登錄注冊 即表示同意《億速云用戶服務條款》

如何使用R語言利用SSR數據構建進化樹

發布時間:2021-11-22 14:54:44 來源:億速云 閱讀:455 作者:柒染 欄目:大數據

這期內容當中小編將會給大家帶來有關如何使用R語言利用SSR數據構建進化樹,文章內容豐富且以專業的角度為大家分析和敘述,閱讀完這篇文章希望大家可以有所收獲。

今天有一位讀者留言問有了SSR的數據,如何畫樹狀圖?

SSR的數據我也是第一次接觸,我看了一下他發給我的數據,是csv格式,每行是一個樣本,每列是一個位點。位點的取值是0,和1。這種格式好像叫做二進制的格式。就是下面這種

如何使用R語言利用SSR數據構建進化樹  
image.png

第一列是樣本名稱,后面每一列是一個位點。

我搜索了一下相關教程,找到了一個。需要借助ape這個包,如果是第一次使用的話需要先安裝

install.packages("ape")
 

后面進化樹可視化還需要用到ggtree,如果是第一次使用也要安裝

BiocManager::install("ggtree")
 

BiocManager第一次使用也得先安裝

install.packages("BiocManager")
 

找到的教程的鏈接是 https://www.biostars.org/p/100432/

 首先是讀入數據
df<-read.csv("SSR_example.csv",header=T,row.names = 1)
   將數據框轉化成矩陣
ssr<-as.matrix(df)
   NJ法構建進化樹
library(ape)
tree<-nj(dist.gene(ssr))
   使用ggtree展示進化樹
library(ggtree)
ggtree(tree)+
 geom_tiplab()+
 xlim(0,6.5)
 
如何使用R語言利用SSR數據構建進化樹  
image.png
 樣本太多可以選擇環形樹展示
ggtree(tree,layout="circular")+
 geom_tiplab2(size=3)
 
如何使用R語言利用SSR數據構建進化樹  
image.png
 去掉枝長看起來可能會好看一點
ggtree(tree,layout="circular",branch.length = "none")+
 geom_tiplab2(size=3)
 
如何使用R語言利用SSR數據構建進化樹  
image.png
 

遇到的問題

這位讀者的數據過多,最后畫樹狀圖如果帶上樣本名字的效果如下

如何使用R語言利用SSR數據構建進化樹  
image.png

這個圖應該如何美化我還真沒有思路。大家如果做過這種圖歡迎留言討論 如何美化會好看一點

我記得好像MEGA也可以利用這種二進制數據構建進化樹了,但是一時想不起來如何做了。

上述就是小編為大家分享的如何使用R語言利用SSR數據構建進化樹了,如果剛好有類似的疑惑,不妨參照上述分析進行理解。如果想知道更多相關知識,歡迎關注億速云行業資訊頻道。

向AI問一下細節

免責聲明:本站發布的內容(圖片、視頻和文字)以原創、轉載和分享為主,文章觀點不代表本網站立場,如果涉及侵權請聯系站長郵箱:is@yisu.com進行舉報,并提供相關證據,一經查實,將立刻刪除涉嫌侵權內容。

AI

确山县| 新津县| 邢台县| 礼泉县| 治县。| 揭阳市| 高州市| 榆中县| 莲花县| 天全县| 若尔盖县| 喀什市| 融水| 盱眙县| 杨浦区| 徐州市| 杭州市| 介休市| 都江堰市| 平远县| 蛟河市| 同江市| 房产| 松潘县| 利辛县| 塘沽区| 无棣县| 连南| 平阴县| 遂平县| 泸州市| 哈尔滨市| 安多县| 台南市| 怀宁县| 彰化县| 祁阳县| 宁津县| 霍林郭勒市| 尉氏县| 乳源|