查找並繪製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（最陡的上升），第二個它有最陡的下降。