导航：首页 > 互联网科技 >

R语言ggplot2两个离散变量画热图分别是怎样的

发表于：2025-01-31 作者：千家信息网编辑

千家信息网最后更新 2025年01月31日，本篇文章为大家展示了R语言ggplot2两个离散变量画热图分别是怎样的，内容简明扼要并且容易理解，绝对能使你眼前一亮，通过这篇文章的详细介绍希望你能有所收获。数据集的部分截图，总共是81行，52列，行

千家信息网最后更新 2025年01月31日R语言ggplot2两个离散变量画热图分别是怎样的

本篇文章为大家展示了R语言ggplot2两个离散变量画热图分别是怎样的，内容简明扼要并且容易理解，绝对能使你眼前一亮，通过这篇文章的详细介绍希望你能有所收获。

数据集的部分截图，总共是81行，52列，行是物种，列是地点，数值代表物种丰度

数据集下载链接 http://userweb.eng.gla.ac.uk/umer.ijaz/bioinformatics/ecological/SPE_pitlatrine.csv

读入数据

abund_table<-read.csv("pitlatrines/SPE_pitlatrine.csv",row.names=1,check.names=FALSE)

这里的参数

row.names=1指定第一列作为数据集的行名
check.names 参数我平时很少用，，查了一下帮助文档，作用是检查每列的名字是否符合规范

比如如下数据集分别设置 check.names参数为T和F大家可以看下效果

数据转置并计算相对频率

abund_table<-t(abund_table)
# Convert to relative frequencies
abund_table <- abund_table/rowSums(abund_table)

宽格式数据转换长格式

library(reshape2)
df<-melt(abund_table)
head(df)
colnames(df)<-c("Samples","Species","Value")

分组求根

library(plyr)
library(scales)

# We are going to apply transformation to our data to make it
# easier on eyes 

#df<-ddply(df,.(Samples),transform,rescale=scale(Value))
df<-ddply(df,.(Samples),transform,rescale=sqrt(Value))

ggplot2画图

library(ggplot2)
p <- ggplot(df, aes(Species, Samples)) + 
  geom_tile(aes(fill = rescale),colour = "white") + 
  scale_fill_gradient(low = "white",high = "steelblue")+
  scale_x_discrete(expand = c(0, 0)) +
  scale_y_discrete(expand = c(0, 0)) + theme(legend.position = "none",axis.ticks = element_blank(),axis.text.x = element_text(angle = 90, hjust = 1,size=5),axis.text.y = element_text(size=5))
p