亚洲激情专区-91九色丨porny丨老师-久久久久久久女国产乱让韩-国产精品午夜小视频观看

溫馨提示×

溫馨提示×

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

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

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

苏尼特左旗| 柏乡县| 泾阳县| 罗源县| 顺昌县| 鸡东县| 永济市| 云南省| 邹城市| 花垣县| 桃园县| 天水市| 习水县| 乾安县| 孝昌县| 疏勒县| 德化县| 罗城| 梨树县| 札达县| 隆化县| 潜江市| 宜宾县| 买车| 房山区| 新民市| 衡水市| 化隆| 庆元县| 郑州市| 鄂托克前旗| 瓦房店市| 平度市| 自治县| 石屏县| 凌云县| 溧水县| 浑源县| 保康县| 繁昌县| 江津市|