找出R中给出3个点的二次方程

Question

R中是否有任何包找到给定3个点的二次方程或者我必须手动完成？

我有一个包含x值的列表和一个包含y值的列表.

x = tail(mylist_x,3)
y = tail(mylist_y,3)

Answer 1

使用线性回归:

x <- 1:3
y <- 0.5 * x^2 + 3 * x - 2

coef(lm(y ~ poly(x, 2, raw = TRUE)))
#(Intercept) poly(x, 2, raw = TRUE)1 poly(x, 2, raw = TRUE)2 
#       -2.0                     3.0                     0.5 

或者如评论中所指出:

solve(cbind(1, x, x^2), y)
#        x      
#-2.0  3.0  0.5 

注意,lm可以处理一些边缘情况(例如,三个非唯一点),而solve不能.