論文
Common anti-cancer therapies induce somatic mutations in stem cells of healthy tissue
https://www.nature.com/articles/s41467-022-33663-5#Sec18
今天的推文我們重復(fù)一下論文中的Figure1a，看起來(lái)非常像一個(gè)表格万牺，很有意思
image.png
我們使用ggplot2來(lái)實(shí)現(xiàn)缓熟，主要分為三個(gè)部分鱼响，文本，形狀憨栽，和背景的灰色條
首先是文本的部分?jǐn)?shù)據(jù)截圖
image.png
形狀的三個(gè)數(shù)據(jù)截圖
image.png
背景灰色條數(shù)據(jù)截圖
image.png
總共5份數(shù)據(jù)式镐，放在同一個(gè)excel文件里，5個(gè)不同的子表格
讀取數(shù)據(jù)并作圖
library(readxl)
library(ggplot2)
library(stringr)
library(tidyverse)
dat01<-read_excel("data/20221009/example_data.xlsx",
sheet = "Sheet1")
dat01 %>%
mutate(label01=case_when(
length(label) == 7 ~ str_pad(label,8,"right"),
TRUE ~ label
)) -> new.dat01
new.dat01
dat02<-read_excel("data/20221009/example_data.xlsx",
sheet = "Sheet2")
dat02
dat03<-read_excel("data/20221009/example_data.xlsx",
sheet = "Sheet3")
dat03
dat04<-read_excel("data/20221009/example_data.xlsx",
sheet = "Sheet4")
dat04
dat05<-read_excel("data/20221009/example_data.xlsx",
sheet = "Sheet5")
dat05
x_labels<-c("","Age","Normal\nTissue","Tissue\nSubtype",
"Treatment\ntype","Treatment\ncycli",
"Sampling after\ntreatment\n(months)",
"Seqenced\nsamples")
ggplot()+
scale_x_continuous(limits = c(0.5,8),
breaks = 1:8,
labels = x_labels,
position = "top")+
geom_text(data=new.dat01,aes(x=x,y=y,label=label01))+
geom_point(data=dat02,aes(x=x,y=y,color=`Normal tissue type`),
size=5)+
scale_color_manual(values = c("Colon"="#d38e91",
"Liver"="#1f639a"))+
ggnewscale::new_scale_color()+
geom_point(data=dat03,aes(x=x,y=y,color=`Tissue subtype`),
size=5)+
scale_color_manual(values = c("Descending colon"="#f6d65b",
"Sigmoid"="#3eada2",
"Rectum"="#eb553a",
"Liver"="#206599"))+
ggnewscale::new_scale_color()+
geom_point(data=dat04,aes(x=x,y=y,color=Treatment),
size=5,shape=15)+
scale_color_manual(values = c("5-FU + platinum"="#f49b5e",
"5-FU + radiation"="#c44657",
"5-FU + platinum + radiation"=
"#89520e"))+
theme_bw()+
theme(axis.title = element_blank(),
panel.grid = element_blank(),
panel.border = element_blank(),
axis.text.y = element_blank(),
axis.ticks = element_blank(),
axis.text.x = element_text(hjust=0.5,vjust=0.5))+
geom_rect(data=dat05,
aes(xmin=-Inf,xmax=Inf,ymin=ymin,ymax=ymax),
fill="gray",alpha=0.2)
image.png
示例數(shù)據(jù)和代碼可以給推文點(diǎn)贊 點(diǎn)擊在看 最后留言獲取
歡迎大家關(guān)注我的公眾號(hào)
小明的數(shù)據(jù)分析筆記本
小明的數(shù)據(jù)分析筆記本 公眾號(hào) 主要分享：1奋蔚、R語(yǔ)言和python做數(shù)據(jù)分析和數(shù)據(jù)可視化的簡(jiǎn)單小例子；2烈钞、園藝植物相關(guān)轉(zhuǎn)錄組學(xué)泊碑、基因組學(xué)、群體遺傳學(xué)文獻(xiàn)閱讀筆記毯欣；3蛾狗、生物信息學(xué)入門學(xué)習(xí)資料及自己的學(xué)習(xí)筆記！
