我已经将4阶多项式曲线拟合到我的数据中,如下所示:
y<-c(-13,16,35,40,28,36,43,33,40,33,22,-5,-27,-31,-29,-25,-26,-31,-26,-24,-25,-29,-23,4)
x<-1:24
#4th order polynomial fit
fit<-lm(y~poly(x,4,raw=TRUE))
plot(x,y,ylim=c(min(y)-10,max(y)+10))
lines(x,predict(fit,data.frame(x=x)),col="red")
abline(h=0,lty=2)
我的最终目标是计算该曲线与零线相交的3个点。
首先,我需要延长曲线拟合的末端,使其第三次超越零线。一旦我完成了这个,我想计算这个等式通过零线的3个点。
答案 0 :(得分:4)
您可以使用predict
功能从拟合模型中获取值。例如
pred <- function(x) predict(fit, newdata=data.frame(x=x))
然后,如果你想要多个根,你可以使用uniroot.all
包中的rootSolve
这样的函数
rootSolve::uniroot.all(pred, c(0,30))
# 1.428859 11.990087 24.420745
将从您的模型中找到0到30之间的根。您也可以多次调用基函数uniroot
。