JZXBK的个人博客分享 http://blog.sciencenet.cn/u/JZXBK

博文

R语言 堆叠柱状图 物种丰度直方图

已有 17878 次阅读 2019-12-20 21:04 |系统分类:科研笔记

#使用ggplot2包绘制物种丰度直方图

#首先安装ggplot2包,并载入

install.packages(ggplot2)

library(ggplot2)


#OTUs相对丰度矩阵表(以97%的相似度为聚类阈值的分析结果)

image.png

#安自己的实验设计整理,第一列是行名(row.names=1)(#门水平top10,others是将剩余低丰度门的进行#了求和)

image.png

#读取数据

sss<-read.csv("sptop10.csv",row.names=1)

#使用factor()函数,并指定levels选项进行排序(这个levels和表格里面的level没有半毛钱关系,它是factor()函数里的一个参数)

sss$level<-factor(sss$level,levels=c(“Robinia”,“Clover”,“Control”,“In-situ”)) #设置顺序

sss$sample<-factor(sss$sample,levels=c(“1”,“2”,“3”)) #设置顺序

#设置了一个颜色向量colors

colors<-c("#984EA3","cyan3","brown3","burlywood2","darkorange","chartreuse4","darkseagreen",

"dodgerblue4","gold","darkolivegreen3","grey50")  #颜色


#绘图

#aes()函数的功能是指定每个变量扮演的角色,这里变量sample的值映射到沿x轴的距离上,变量value的值映#射到沿y轴的距离上。

#ggplot()函数指定要绘制的数据源和变量

#几何函数geom_bar()指定这些变量以柱状图的形式展示,参数position="fill"垂直的堆叠条形图并规范其高

#度相等

ggplot(sss,aes(x=sample,y=value,fill=species))+

geom_bar(position = "fill",stat="identity")+

theme_bw()+

facet_grid(.~level)+ #设置分组;每个level水平的独立图,配置成一个单行。

scale_fill_manual(values=colors)+

labs(x="Depth(cm)",y="Relative Abundance",fill="Phylum")#添加标题

结果如下图


image.png




https://blog.sciencenet.cn/blog-3419243-1210867.html

上一篇:R 语言 常用的小操作 小技巧
下一篇:α多样性指数变化趋势
收藏 IP: 124.23.134.*| 热度|

0

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

数据加载中...
扫一扫,分享此博文

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

GMT+8, 2024-12-26 20:02

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部