我有两个人的生物学数据,我用R作为散点图使用ggplot绘制图形,如下所示:
p1<-ggplot(data, aes(meth_matrix$sample1, meth_matrix$sample3)) +
geom_point() +
theme_minimal()
哪作得很完美,但我想为它添加一些行:将散点图分成两半的abline:
p1 + geom_abline(color="blue")
我的问题是:如何绘制两条与该对角线平行的红线(y截距为0.2,斜率与蓝线相同)??
另外:如何使用ggplot在类似的散点图(它看起来像水平散点图)中绘制两个样本的差异?现在我只能用以下情节来做:
dif_samples<-meth_matrix$sample1- meth_matrix$sample3
plot(dif_samples, main="difference",
xlab="CpGs ", ylab="Methylation ", pch=19)
(我也想把水平蓝线和红线paralllel加到蓝线上)
请帮助!!!
非常感谢。
答案 0 :(得分:2)
您可以在geom_abline()
函数中指定斜率和截距。我将使用ggplot2
附带的iris数据集来说明:
# I'll use the iris dataset. I normalise by dividing the variables by their max so that
# a line through the origin will be visible
library(ggplot2)
p1 <- ggplot(iris, aes(Sepal.Length/max(Sepal.Length), Sepal.Width/max(Sepal.Width))) +
geom_point() + theme_minimal()
# Draw lines by specifying their slopes and intercepts. since all lines
# share a slope I just give one insted of a vector of slopes
p1 + geom_abline(intercept = c(0, .2, -.2), slope = 1,
color = c("blue", "red", "red"))
我不清楚你想要第二个情节的确切内容,但你可以直接在ggplot()
的调用中绘制差异,你可以添加geom_hline()
的水平线:< / p>
# Now lets plot the difference between sepal length and width
# for each observation
p2 <- ggplot(iris, aes(x = 1:nrow(iris),
y = (Sepal.Length - Sepal.Width) )) +
geom_point() + theme_minimal()
# we'll add horizontal lines -- you can pick values that make sense for your problem
p2 + geom_hline(yintercept = c(3, 3.2, 2.8),
color = c("blue", "red", "red"))
由reprex package(v0.2.0)创建于2018-03-21。