查找并绘制r中曲线的第二个最陡峭部分

Question

我有一些数据，我已经对使用r进行了7阶多项式回归;该曲线类似于钟形曲线。

我需要找到曲线的两个最陡的部分。到目前为止，我已经设法找到并显示最陡峭的上坡，现在我需要最陡峭的下坡。

这是我到目前为止有：

attach(DATA1) 
names(DATA1) 

"dist" "ohms" 

plot(dist,ohms) 
reg<-lm(ohms~poly(dist,7)) 
summary(reg) 
lines(smooth.spline(dist,predict(reg))) 
xv<-seq(min(dist),max(dist),0.02) 
yv<-predict(reg,list(dist=xv)) 
lines(xv,yv,lwd=2) 
plot(xv,yv,type="l") 
xv[which(abs(diff(yv))==max(abs(diff(yv))))] 
abline(v=xv[which(abs(diff(yv))==max(abs(diff(yv))))]) 

有谁知道如何找坡最陡的下降部分？

感谢

Answer 1

我认为你在这里有一个问题：

abs(diff(yv))==max(abs(diff(yv))) 

在这一行，diff(yv)几乎是yv的衍生物。通过取其绝对值，你忽略了是否增加或减少yv，所以你会得到最陡峭的上升或下降，但你不知道哪一个。

既然你有兴趣在寻找最陡的上升和下降试试这个：

diff(yv) == max(diff(yv)) 
diff(yv) == min(diff(yv)) 

首先应该当YV的变化是最大返回true（最陡的上升），第二个它有最陡的下降。