我画了一个等高线图并添加了一条对角线:
library(MASS)
library(fields)
x <- c(21.06, 28.89, 23.00, 23.61, 23.61, 22.83, 30.44)
y <- c(26.56, 24.00, 13.06, 18.61, 18.61, 14.17, 25.33)
z <- kde2d(x, y, n=32, lims = c(0,32,0,32))
contour(z, col = "red", main = "Density estimation: contour plot",
las=0,
plot.title=
{
title(xlab=expression(alpha),cex.lab=2)
mtext(expression(beta),2,cex=2,line=3,las=1)
}
)
abline(0, 1, col = "red", lwd = 2)
我想我更喜欢填充的轮廓,但现在线条已关闭:
filled.contour(z, plot.title={
title(main = "Density estimation: contour plot")
title(xlab=expression(alpha),cex.lab=2)
mtext(expression(beta),2,cex=2,line=3,las=1)
})
abline(0, 1, col = "red", lwd = 2)
最佳答案
你可以把 abline
扔进去:
library(MASS)
x <- c(21.06, 28.89, 23.00, 23.61, 23.61, 22.83, 30.44)
y <- c(26.56, 24.00, 13.06, 18.61, 18.61, 14.17, 25.33)
z <- kde2d(x, y, n = 32, lims = c(0, 32, 0, 32))
filled.contour(z, plot.title={
title(main = "Density estimation: contour plot")
title(xlab=expression(alpha),cex.lab=2)
mtext(expression(beta),2,cex=2,line=3,las=1)
abline(0, 1, col = "red", lwd = 2)
})
或者(正如评论中所建议的那样),您最好还是执行以下操作:
filled.contour(
z,
plot.title = {
title(main = "Density estimation: contour plot")
title(xlab = expression(alpha), cex.lab = 2)
mtext(expression(beta), 2, cex = 2, line = 3, las = 1)
},
plot.axes = {
abline(0, 1, col = "red", lwd = 2)
}
)
只是因为在 plot.title
中放置一个 abline
感觉有点奇怪。
关于r - 将 abline 添加到填充等高线图,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/53907534/