有的時(shí)候我們用ggplot2畫(huà)了圖冲粤,比如說(shuō)散點(diǎn)圖,那么如何在一堆點(diǎn)里標(biāo)出我們關(guān)注的點(diǎn)呢菜拓?
> dat = read.csv("xxxx.csv",header = TRUE,sep = ",")
> colnames(dat)
[1] "gene_id" "geneModuleMembership" "geneTraitSignificanc
比如:我想把后兩列用散點(diǎn)圖表示出來(lái),那么:
> library(ggplot2)
> ggplot(dat,aes(x=geneModuleMembership,y=geneTraitSignificance)) +
geom_point(shape=1,colour= "turquoise",size = 2) + xlab("x_axis title")+ ylab("y_axis title")+ ggtitle("picture title")
各種參數(shù):
ggplot(數(shù)據(jù)dataframe笛厦,aes(你想畫(huà)的列)) +
geom_point(點(diǎn)的形狀纳鼎,點(diǎn)的顏色,點(diǎn)的大猩淹埂)+ xlab(x軸標(biāo)題) + ylab(y軸標(biāo)題) + ggtitle(整個(gè)圖的標(biāo)題)
如果想改變點(diǎn)的形狀贱鄙,這里是數(shù)字對(duì)應(yīng)的形狀:ggplot2 point shapes
覺(jué)得圖的背景灰色格子很難看嗎?可以去掉:
> ggplot(dat,aes(x=geneModuleMembership,y=geneTraitSignificance)) +
geom_point(shape=1,colour= "turquoise",size = 2) + xlab("x_axis title")+ ylab("y_axis title")+ ggtitle("picture title")+
theme_classic() #這個(gè)參數(shù)是去掉背景灰色格子的
圖畫(huà)的差不多了姨谷,但是我想在圖里標(biāo)注一個(gè)特定的點(diǎn):
ggplot(dat,aes(x=geneModuleMembership,y=geneTraitSignificance)) +
geom_point(shape=1,colour="turquoise",size = 2) + xlab("x_axis title")+ ylab("y_axis title")+ ggtitle("picture title") +
theme_classic()+
geom_point(aes(0.9071973,0.89379209),color="red",size=2) #這一行是對(duì)特定點(diǎn)標(biāo)注的逗宁,你需要指出它對(duì)應(yīng)的x,y軸的坐標(biāo),還可以指定顏色和大小
給特定點(diǎn)添加注釋:
>ggplot(dat,aes(x=geneModuleMembership,y=geneTraitSignificance)) +
geom_point(shape=1,colour="turquoise",size = 2) + xlab("x_axis title")+ ylab("y_axis title")+ ggtitle("picture title") +
theme_classic()+
geom_point(aes(0.9071973,0.89379209),color="red",size=2) +
geom_text(aes(0.9071973,0.89379209,label="Junb"),hjust=-1,vjust=1)#這一行是給特定點(diǎn)加注釋的
NOTE: 如果你添加了很多個(gè)特定點(diǎn)菠秒,而這些點(diǎn)離的很近疙剑,注釋信息重疊了,你可以通過(guò)改變注釋信息的位置來(lái)調(diào)整践叠。調(diào)整注釋信息位置的參數(shù)是hjust和vjust
hjust: 調(diào)整左右位置言缤。 以紅點(diǎn)為中心,該數(shù)值越大禁灼,注釋信息越靠左管挟,0.5時(shí)在中間
vjust: 調(diào)整上下位置。以紅點(diǎn)為中心弄捕,該數(shù)值越大僻孝,注釋信息越靠下,0.5時(shí)在中間
比如現(xiàn)在守谓,基因名在紅點(diǎn)的右側(cè)偏下穿铆,我想把它調(diào)整到紅點(diǎn)的左上方,那么:
> ggplot(dat,aes(x=geneModuleMembership,y=geneTraitSignificance)) +
geom_point(shape=1,colour="turquoise",size = 2) + xlab("x_axis title")+ ylab("y_axis title")+ ggtitle("picture title") +
theme_classic()+
geom_point(aes(0.9071973,0.89379209),color="red",size=2) +
geom_text(aes(0.9071973,0.89379209,label="Junb"),hjust=1,vjust=-0.5)