03 ggplot基本图形
2022-12-09 本文已影响0人
Jachin111
散点图
library(ggplot2)
df1 <- read.csv("scatter_plot_example.csv",header=T)
head(df1)
![](https://img.haomeiwen.com/i9248907/8eebaf823d749eb0.png)
ggplot(data=df1,aes(x=var1,y=var2)) + geom_point()
![](https://img.haomeiwen.com/i9248907/2e4b43d7bf741ca8.png)
ggplot(data=df1,aes(x=var1,y=var2)) +
geom_point(color="red",size=10,shape=18,alpha=0.5)
![](https://img.haomeiwen.com/i9248907/584c604cc4434320.png)
ggplot(data=df1,aes(x=var1,y=var2,color=var3,shape=var3)) +
geom_point(size=5)
![](https://img.haomeiwen.com/i9248907/ab3f3bf69669d1cd.png)
ggplot(data=df1,aes(x=var1,y=var2,color=var3,shape=var3)) +
geom_point(size=5) +
scale_color_manual(values=c("steelblue","yellowgreen","violetred1"))
![](https://img.haomeiwen.com/i9248907/84aea11e91694fb0.png)
ggplot(data=df1,aes(x=var1,y=var2,color=var3,shape=var3)) +
geom_point(size=5) +
scale_color_manual(values=c("steelblue","yellowgreen","violetred1")) +
scale_shape_manual(values=c(9,10,11))
![](https://img.haomeiwen.com/i9248907/72594f9aa64b0b7b.png)
气泡图
df2 <- read.csv("bubble_plot_example.csv",header=T)
head(df2)
![](https://img.haomeiwen.com/i9248907/ecfad38c77d650dd.png)
ggplot(data=df2,aes(x=var1,y=var2,color=var3,size=var4)) +
geom_point(alpha=0.8) +
scale_size_continuous(range=c(1,50)) +
theme_bw() +
theme(legend.position="none") +
ylim(0,0.4) +
scale_color_manual(values=c("steelblue","yellowgreen","violetred1"))
![](https://img.haomeiwen.com/i9248907/0856703f5f7ab279.png)
折线图
df3 <- read.csv("line_example.csv",header=T)
head(df3)
![](https://img.haomeiwen.com/i9248907/0bb022a0985e2a47.png)
ggplot(data=df3,aes(x=time_point,y=value)) +
geom_line() +
geom_point() +
geom_errorbar(aes(ymin=value-sd,ymax=value+sd),width=0.15) +
theme_bw()
![](https://img.haomeiwen.com/i9248907/b7308d1c2bb02d78.png)
柱形图
df4 <- read.csv("bar_plot_example.csv",header=T)
head(df4)
![](https://img.haomeiwen.com/i9248907/39774a4d93baa85e.png)
ggplot(data=df4,aes(x=var1,y=var2)) + geom_col()
![](https://img.haomeiwen.com/i9248907/646021cabc367a2a.png)
ggplot(data=df4,aes(x=var1,y=var2)) +
geom_col(color="red",fill="lightblue")
![](https://img.haomeiwen.com/i9248907/df78a29184b843d0.png)
ggplot(data=df4,aes(x=var1,y=var2)) +
geom_col(aes(fill=var1)) +
theme_bw() +
scale_fill_manual(values=c("#008fd5","#ff2700","#77ab43","#ffd700")) +
scale_y_continuous(expand=c(0,0),limits=c(0,8))
![](https://img.haomeiwen.com/i9248907/8c6f2bbf97de33f5.png)
df5 <- read.csv("bar_plot_example_1.csv",header=T)
head(df5)
![](https://img.haomeiwen.com/i9248907/a538bd3359594993.png)
ggplot(data=df5,aes(x=var3,y=var2,fill=var1)) +
geom_bar(stat="identity") +
theme_bw() +
scale_fill_manual(values=c("steelblue","yellowgreen","violetred1")) +
scale_y_continuous(expand=c(0,0),limits=c(0,20))
![](https://img.haomeiwen.com/i9248907/70e12df7222e648c.png)
ggplot(data=df5,aes(x=var3,y=var2,fill=var1)) +
geom_bar(stat="identity",position="dodge") +
scale_fill_manual(values=c("steelblue","yellowgreen","violetred1")) +
scale_y_continuous(expand=c(0,0),limits=c(0,10)) +
theme_bw()
![](https://img.haomeiwen.com/i9248907/ca5bcceaf509fd32.png)
箱线图
df6 <- read.csv("ToothGrowth.csv",header=T)
head(df6)
table(df6$dose)
![](https://img.haomeiwen.com/i9248907/b154b80dd7ac00c1.png)
df6$dose <- factor(df6$dose,levels=c("0.5","1","2"))
ggplot(data=df6,aes(x=dose,y=len,fill=dose)) +
geom_boxplot() +
scale_fill_manual(values=c("steelblue","yellowgreen","violetred1")) +
theme_bw() +
stat_boxplot(geom="errorbar",width=0.2)
![](https://img.haomeiwen.com/i9248907/4e3d2190883efba5.png)
小提琴图
ggplot(data=df6,aes(x=dose,y=len,fill=dose)) +
geom_violin() +
theme_bw() +
scale_fill_manual(values=c("steelblue","yellowgreen","violetred1"))
![](https://img.haomeiwen.com/i9248907/7dbc7274ef00b307.png)
ggplot(data=df6,aes(x=dose,y=len,fill=dose)) +
geom_violin() +
geom_boxplot(width=0.3) +
scale_fill_manual(values=c("steelblue","yellowgreen","violetred1")) +
theme_bw() +
stat_boxplot(geom="errorbar",width=0.2)
![](https://img.haomeiwen.com/i9248907/3824004c998a82cb.png)
热图
df7 <- read.csv("pheatmap_example_data.csv",header=T)
head(df7)
![](https://img.haomeiwen.com/i9248907/752e4c5355d00405.png)
library(reshape2)
df7.1 <- melt(df7,id.vars='gene_name')
head(df7.1)
ggplot(data=df7.1,aes(x=variable,y=gene_name,fill=value)) +
geom_tile(color="black") +
theme(axis.text.x=element_text(angle=90,hjust=0.5,vjust=0.5),axis.ticks=element_blank()) +
scale_fill_gradient(low="steelblue",high="violetred1") +
labs(x=NULL,y=NULL) +
geom_text(aes(label=round(value,1)))
![](https://img.haomeiwen.com/i9248907/0fe818c5955dab1a.png)