昨天公眾號后臺有人留言作圖募闲，示例圖如下
image.png
我選擇使用R語言的ggplot2來實現(xiàn)步脓，這個是箱線圖和熱圖的拼接，右側(cè)的熱圖可以借助
geom_point()函數(shù)實現(xiàn)，將點的形狀改為正方塊沪编，數(shù)值按照正負(fù)來映射顏色呼盆，按照一定的數(shù)值來映射大小∫侠基本的思路有了访圃。今天的推文記錄一下如何實現(xiàn)這個熱圖
首先構(gòu)造一份數(shù)據(jù)集
set.seed("20200407")
x<-seq(-1,1,by=0.001)
df<-data.frame(matrix(sample(x,90),nrow = 6))
head(df)
colnames(df)<-paste0("Y",1:15)
rownames(df)<-paste0("X",1:6)
write.csv(df,file = "20210407_example.csv",quote = F)
image.png
這樣是寬格式數(shù)據(jù)，ggplot2作圖需要整理成長格式
df<-read.csv("20210407_example.csv",header=T)
head(df)
df1<-reshape2::melt(df,id.vars="Sample")
head(df1)
image.png
構(gòu)造數(shù)據(jù)
library(dplyr)
df1%>%
mutate(group_1 = case_when(
value <= 0 ~ "A",
TRUE ~ "B"
))%>%
mutate(group_2=case_when(
value >= -1 & value < -0.7 ~ "[-1,-0.7)",
value >= -0.7 & value < -0.5 ~ "[-0.7,-0.5)",
value >= -0.5 & value < -0.3 ~ "[-0.5,-0.3)",
value >= -0.3 & value <= 0 ~ "[-0.3,0]",
value > 0 & value <= 0.3 ~ "(0,0.3)",
value > 0.3 & value <= 0.5 ~ "(0.3,0.5]",
value > 0.5 & value <= 0.7 ~ "(0.5,0.7]",
value > 0.7 & value <= 1 ~ "(0.7,1]",
))%>%
mutate(value_1=case_when(
value >= -1 & value < -0.7 ~ -0.8,
value >= -0.7 & value < -0.5 ~ -0.6,
value >= -0.5 & value < -0.3 ~ -0.4,
value >= -0.3 & value <= 0 ~ -0.2,
value > 0 & value <= 0.3 ~ 0.2,
value > 0.3 & value <= 0.5 ~ 0.4,
value > 0.5 & value <= 0.7 ~ 0.6,
value > 0.7 & value <= 1 ~ 0.8,
)) -> df2
構(gòu)造數(shù)據(jù)用于添加線段
df3<-data.frame(
x = seq(1.5,5.5,1),
xend = seq(1.5,5.5,1),
y = -Inf,
yend = Inf
)
df3
df4<-data.frame(
x = -Inf,
xend = Inf,
y = seq(1.5,14.5,1),
yend = seq(1.5,14.5,1)
)
畫圖
library(ggplot2)
ggplot(data=df2,aes(x=Sample,y=variable))+
geom_point(aes(size=abs(value_1),
color=factor(value_1)),
shape=15)+
scale_color_manual(values = c(rep("#fe0000",4),rep("#009ccc",4)))+
theme_bw()+
theme(panel.grid = element_blank(),
panel.border = element_rect(color="grey"),
axis.ticks = element_blank())+
geom_segment(data=df3,aes(x=x,xend=xend,y=y,yend=yend),
color="grey")+
geom_segment(data=df4,aes(x=x,xend=xend,y=y,yend=yend),
color="grey")+
scale_size_continuous(range = c(2,10))+
scale_y_discrete(position = "right")+
labs(x=NULL,y=NULL)
image.png
這里遇到的問題是如何讓圖例想推文開頭的圖片一樣：一個圖例就能夠既反映大小又有顏色呢相嵌？而且圖例的排布是從打到小然后又從小到大腿时。這個如何用代碼實現(xiàn)呢？暫時不知道了饭宾，手動編輯吧批糟，可能得費一些時間
image.png
這樣手動調(diào)整還挺費時間的，今天就先到這里了看铆，好好想想有沒有用代碼實現(xiàn)的辦法徽鼎，歡迎大家留言討論呀!
歡迎大家關(guān)注我的公眾號
小明的數(shù)據(jù)分析筆記本
小明的數(shù)據(jù)分析筆記本 公眾號 主要分享：1、R語言和python做數(shù)據(jù)分析和數(shù)據(jù)可視化的簡單小例子弹惦；2否淤、園藝植物相關(guān)轉(zhuǎn)錄組學(xué)、基因組學(xué)棠隐、群體遺傳學(xué)文獻(xiàn)閱讀筆記石抡；3、生物信息學(xué)入門學(xué)習(xí)資料及自己的學(xué)習(xí)筆記啰扛！
