MalacoWiki

wiki of malacology

RUA | 使用R进行聚类分析

RUA(use R for Unlimitedly Analysis)系列,系统发育树的本质可以理解为聚类分析,此处,我们使用一些连续性数值进行聚类分析,可以方便直观地看出个体之间在所选取数据上的相互关系。

计算数值

对已知数据的计算

setwd("C:/Users/月季花中的星辰/Desktop/R Phylo") #设置路径
Data = read.table("data.txt") #读取文件,定义
dists = dist(Data,method = "euclidean") #设置欧氏距离,来进行分析
hc = hclust(dists, method = "ave") #方法是average

树的展现

第一种展现树的方式

dend = as.dendrogram(hc)
plot(dend, type = "rectangle", ylab="Height", main="Cluster Dendrogram")

第二种展现树的方式

install.packages('ape') #第一次使用请安装
library(ape) #载入ape包
class(hc) #nc是你使用的hclust的定义
tree = phylo(hc)
my_tree =as.phylo(hc)
write.tree(phy=my_tree, file="tree.newick") #设置的路径内找到

原理

施工中

引用

Zhang G. 2020. Clustering by R. malacology.net [Access Time]

Last updated on 6 Dec 2020
Published on 6 Dec 2020

Powered by Arch. Copyright by Guoyi.