我有两个光栅图像的散点图。我需要的是在散点图周围创建一个三角形。到目前为止我拥有它
minveg <- quantile(veg, 0.01, na.rm=TRUE)
mintemp <- quantile(temps, 0.01, na.rm=TRUE)
plot(veg, temps, xlab='veg', ylab='temp')
abline(h=mintemp, v=minveg)
这给我以下图片
我需要的是沿着散点图的第三条线,并用另外两条线完成三角形。我确定我可以添加一个方程式,只是看着它并符合我的最佳想法,但我必须为数百张图像做这个,并且正在寻找一种更自动化的方法。任何建议都有帮助。谢谢!
答案 0 :(得分:1)
请在提出R问题时提供可重现的示例!我们没有'veg'和'temps'。
我认为您可以使用分位数回归。
library(quantreg)
# example data
set.seed(0)
x = rnorm(1000)
y = rnorm(1000) - x
minx <- quantile(x, 0.01, na.rm=TRUE)
miny <- quantile(y, 0.01, na.rm=TRUE)
plot(x,y, xlim=c(-6,6), ylim=c(-6,6))
abline(h=minx, v=miny, col='red')
abline(rq(y ~ x, tau=.99), col='blue')