r - 拟合具有已知截距的多项式

Question

我正在使用 lm(y~poly(x,2)) 来将二阶多项式拟合到我的数据中。但我只是找不到一种方法来指定已知的截距值。如何使用 lm 拟合具有已知截距值(例如“k”)的多项式模型？

Answer 1

 lm(y~-1+x+I(x^2)+offset(k))

应该这样做。

• -1 抑制自动添加的拦截项
• x 添加线性项
• I(x^2) 添加一个二次项；需要 I() ，以便 R 将 ^2 解释为平方，而不是在 x 与其自身之间进行交互(根据公式规则相当于单独的 x)
• offset(k) 添加已知常数截距

我不知道poly(x,2)-1是否可以消除截距；你可以尝试一下看看。从数据中减去偏移量应该可以正常工作，但 offset(k) 可能会稍微更明确一些。您可能必须将 k 设为向量(即将其复制到数据集的长度，或者更好地将其作为数据集中的一列并使用 data=.. .