ggplot correlation with values shown¶

[1]:

df = read.csv("Cas9_affinity_function_HS2_7gRNAs_reformat.txt",sep="\t")

[2]:

head(df)

Cas9_affinity	HBE	cat
264.1604	0.4222148	HBE
1122.6194	0.6088005	HBE
591.9303	0.3578722	HBE
986.7332	0.5000359	HBE
489.4101	0.4793091	HBE
2333.3105	0.5794054	HBE

[3]:

library(ggpubr)

Warning message:
"package 'ggpubr' was built under R version 3.5.3"Loading required package: ggplot2
Warning message:
"package 'ggplot2' was built under R version 3.5.3"Loading required package: magrittr

[4]:

library(ggplot2)
ggplot(df,aes(x = Cas9_affinity, y = HBE)) +
geom_smooth(data = subset(df, cat %in% c("HBE")),fill="blue",method='lm')+
geom_point(data = subset(df, cat %in% c("HBE")),color="blue")+
geom_smooth(data = subset(df, cat %in% c("HBG")),fill="green",method='lm')+
geom_point(data = subset(df, cat %in% c("HBG")),color="green")+
geom_smooth(data = subset(df, cat %in% c("HBB")),fill="red",method='lm')+
geom_point(data = subset(df, cat %in% c("HBB")),color="red")
# +
# stat_cor(data = subset(df, cat %in% c("HBE")),method = "pearson", label.x = 3, label.y = 30,color="red")

../_images/jupyter_notebooks_ggplot2.correlation_4_0.png

[5]:

library(ggplot2)
ggplot(df,aes(x = Cas9_affinity, y = HBE)) +
geom_smooth(data = subset(df, cat %in% c("HBE")),fill="blue",color="blue",method='lm')+
geom_point(data = subset(df, cat %in% c("HBE")),color="blue")+
stat_cor(data = subset(df, cat %in% c("HBE")),method = "pearson", label.x = 250, label.y = 0.95,color="blue")+
geom_smooth(data = subset(df, cat %in% c("HBG")),fill="green",color="green",method='lm')+
geom_point(data = subset(df, cat %in% c("HBG")),color="green")+
stat_cor(data = subset(df, cat %in% c("HBG")),method = "pearson", label.x = 250, label.y = 0.85,color="green")+
geom_smooth(data = subset(df, cat %in% c("HBB")),fill="red",color="red",method='lm')+
geom_point(data = subset(df, cat %in% c("HBB")),color="red")+
stat_cor(data = subset(df, cat %in% c("HBB")),method = "pearson", label.x = 250, label.y = 0.75,color="red")

../_images/jupyter_notebooks_ggplot2.correlation_5_0.png

[6]:

library(ggplot2)
library(ggpubr)
ggplot(df,aes(x = Cas9_affinity, y = HBE)) +
geom_smooth(data = subset(df, cat %in% c("HBE")),fill="blue",se=FALSE,color="blue",method='lm')+
geom_point(data = subset(df, cat %in% c("HBE")),color="blue")+
stat_cor(data = subset(df, cat %in% c("HBE")),method = "pearson", label.x = 250, label.y = 0.95,color="blue")+
geom_smooth(data = subset(df, cat %in% c("HBG")),fill="green",se=FALSE,color="green",method='lm')+
geom_point(data = subset(df, cat %in% c("HBG")),color="green")+
stat_cor(data = subset(df, cat %in% c("HBG")),method = "pearson", label.x = 250, label.y = 0.85,color="green")+
geom_smooth(data = subset(df, cat %in% c("HBB")),fill="red",se=FALSE,color="red",method='lm')+
geom_point(data = subset(df, cat %in% c("HBB")),color="red")+
stat_cor(data = subset(df, cat %in% c("HBB")),method = "pearson", label.x = 250, label.y = 0.75,color="red")+
theme_classic()+
ggsave("figure.png",dpi=300)

Saving 6.67 x 6.67 in image

../_images/jupyter_notebooks_ggplot2.correlation_6_1.png

[10]:

library(ggplot2)
library(ggpubr)
ggplot(df,aes(x = Cas9_affinity, y = HBE)) +
geom_smooth(data = subset(df, cat %in% c("HBE")),fill="blue",se=FALSE,color="blue",method='lm')+
geom_point(data = subset(df, cat %in% c("HBE")),color="blue")+
stat_cor(data = subset(df, cat %in% c("HBE")),method = "pearson", label.x = 250, label.y = 0.95,color="blue",aes(label = paste('HBE',..r.label.., ..p.label.., sep = "~`,`~")))+
geom_smooth(data = subset(df, cat %in% c("HBG")),fill="green",se=FALSE,color="green",method='lm')+
geom_point(data = subset(df, cat %in% c("HBG")),color="green")+
stat_cor(data = subset(df, cat %in% c("HBG")),method = "pearson", label.x = 250, label.y = 0.85,color="green",aes(label = paste('HBG',..r.label.., ..p.label.., sep = "~`,`~")))+
geom_smooth(data = subset(df, cat %in% c("HBB")),fill="red",se=FALSE,color="red",method='lm')+
geom_point(data = subset(df, cat %in% c("HBB")),color="red")+
stat_cor(data = subset(df, cat %in% c("HBB")),method = "pearson", label.x = 250, label.y = 0.75,color="red",aes(label = paste('HBB',..r.label.., ..p.label.., sep = "~`,`~")))+
theme_classic()+
ggsave("figure.png",dpi=300)

Saving 6.67 x 6.67 in image

../_images/jupyter_notebooks_ggplot2.correlation_7_1.png