二次函数的 R 公式语法 [已关闭]

Question

我有一个二次模型，我正在尝试适应它

mgcv::gam

，我需要一些帮助将其编码为 r

formula

：

到目前为止我的尝试：

model <- mgcv::gam(W ~ 1 + z + I(0.5*z^2))

您认为这是公式的正确语法吗？也就是说，

formula

与我尝试使用的方程匹配吗？我不太确定

0.5

是否放置正确。

Answer 1

公式中的

1

将被解释为要估计的未知截距，而不是常数1。我认为

W ~ -1 + offset(rep(1, length(z))) + z + I(0.5*z^2)

会做你想要的。

-1

指定没有截距，并且

offset(rep(1, length(z)))

充当方程中的 1。这适用于

stats::glm

，因此根据帮助页面，它也应该适用于

mcgv::gam

，但我不熟悉

mcgv

，因此可能存在我不知道的怪癖。