沉闷科学的掘墓人分享 http://blog.sciencenet.cn/u/Bearjazz

博文

在R语言中使用遗传距离矩阵构建简单的系统发育树

已有 19941 次阅读 2013-9-11 20:56 |个人分类:我的研究|系统分类:科研笔记| R语言, 聚类, 系统发育树, 遗传距离矩阵, NJ树

熊荣川

xiong rongchuan

六盘水师范学院生物信息学实验室

xiongrongchuan@126.com

http://blog.sciencenet.cn/u/Bearjazz

一般的遗传距离矩阵如下图

对这样的矩阵稍作修改,删除标题行列,保存为“Fst.csv”

library(ape)  #需要用到ape程序包
M2 = read.csv("Fst.csv",head = F)  #输入表格
dimnames(M2) <- list(1:8, 1:8)  #为表格设置行列名称

M2<-as.matrix(M2)   #将表格转化成矩阵格式
tr2 <- bionj(M2)  #构建NJ树
pl树ot(tr2, "u")   #输出系统发育

 



https://blog.sciencenet.cn/blog-508298-724117.html

上一篇:解决photoshop保存不了jpg格式问题
下一篇:解决线粒体37个基因下载问题
收藏 IP: 119.78.81.*| 热度|

0

该博文允许注册用户评论 请点击登录 评论 (4 个评论)

数据加载中...

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2024-11-8 00:31

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部