R语言做主坐标举例分析

发布时间:2021-11-22 15:54:42 作者:iii
来源:亿速云 阅读:270

本篇内容介绍了“R语言做主坐标举例分析”的有关知识,在实际案例的操作过程中,不少人都会遇到这样的困境,接下来就让小编带领大家学习一下如何处理这些情况吧!希望大家仔细阅读,能够学有所成!

使用R语言vegan包里的varespec数据集

首先加载vegan这个包

library(vegan)
 

然后通过data()函数来获得这个数据

data("varespec")
 

这个数据集是24行44列,每行是一个样地,每列是样地里的植物种类。数值具体代表什么我还没有看明白,现在的理解就是度量这种植物在这个样地丰富度的一个标准吧。

 下面开始实际操作

df<-varespec
rownames(df)<-paste0("site",1:24)
 
bray_dist<-vegdist(df,method = "bray")
 
library(ape)
df.pcoa<-pcoa(bray_dist,correction = "cailliez")
 

df.pcoa$vectors能够获得用于画图的数据

R语言做主坐标举例分析  
image.png

df.pcoa$values可以获得坐标轴上显示的百分比

R语言做主坐标举例分析  
image.png
 最后用ggplot2来画这个图
df.plot<-data.frame(df.pcoa$vectors)
head(df.plot)
library(ggplot2)
x_label<-round(df.pcoa$values$Rel_corr_eig[1]*100,2)
y_label<-round(df.pcoa$values$Rel_corr_eig[2]*100,2)
x_label
y_label
ggplot(data=df.plot,aes(x=Axis.1,y=Axis.2))+
  geom_point()+
  theme_bw()+
  theme(panel.grid = element_blank())+
  geom_vline(xintercept = 0,lty="dashed")+
  geom_hline(yintercept = 0,lty="dashed")+
  labs(x=paste0("PCoA1 ",x_label,"%"),
       y=paste0("PCoA2 ",y_label,"%"))
 
R语言做主坐标举例分析  
image.png

通过上图我们可以看到这些样地大体上可以分为两组,如果自己手头有样地的分组数据就可以看看这个结果是不是和自己的分组数据一致。

下面人为的给他分个组,然后添加一个表示分组的椭圆

df.plot$group<-ifelse(df.plot$Axis.1<0,"AAA","BBB")
ggplot(data=df.plot,aes(x=Axis.1,y=Axis.2,
                        color=group,shape=group))+
  geom_point(size=5)+
  theme_bw()+
  theme(panel.grid = element_blank())+
  geom_vline(xintercept = 0,lty="dashed")+
  geom_hline(yintercept = 0,lty="dashed")+
  labs(x=paste0("PCoA1 ",x_label,"%"),
       y=paste0("PCoA2 ",y_label,"%"))+
  stat_ellipse(data=df.plot,
               geom = "polygon",
               aes(fill=group),
               alpha=0.3)+
  scale_fill_manual(values = c("#e31a1c","#1f78b4"))
 
R语言做主坐标举例分析  
image.png

“R语言做主坐标举例分析”的内容就介绍到这里了,感谢大家的阅读。如果想了解更多行业相关的知识可以关注亿速云网站,小编将为大家输出更多高质量的实用文章!

推荐阅读:
  1. mysql举例分析
  2. mysqldump流程举例分析

免责声明:本站发布的内容(图片、视频和文字)以原创、转载和分享为主,文章观点不代表本网站立场,如果涉及侵权请联系站长邮箱:is@yisu.com进行举报,并提供相关证据,一经查实,将立刻删除涉嫌侵权内容。

r语言

上一篇:怎样用R语言ggplot2+ggtree+aplot画气泡图组合聚类树图

下一篇:c语言怎么实现含递归清场版扫雷游戏

相关阅读

您好,登录后才能下订单哦!

密码登录
登录注册
其他方式登录
点击 登录注册 即表示同意《亿速云用户服务条款》