二次函数的 R 语法

Question

我需要一些帮助来使用mgcv::gam()制定

二次模型

：

到目前为止我的尝试：

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

您认为这是正确的表述吗？我不太确定

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

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