91超碰碰碰碰久久久久久综合_超碰av人澡人澡人澡人澡人掠_国产黄大片在线观看画质优化_txt小说免费全本

溫馨提示×

溫馨提示×

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

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

如何使用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

苍溪县| 合水县| 古蔺县| 南昌市| 汕头市| 永德县| 巴中市| 余江县| 永丰县| 阿巴嘎旗| 澳门| 阳山县| 容城县| 什邡市| 无锡市| 亚东县| 霍邱县| 兴仁县| 娄烦县| 宣城市| 惠东县| 贵州省| 新和县| 乐清市| 溆浦县| 安远县| 荥阳市| 舞钢市| 如皋市| 图们市| 菏泽市| 交口县| 甘德县| 张家界市| 永丰县| 肇源县| 宝鸡市| 兴仁县| 长岭县| 乌拉特后旗| 黄大仙区|