我正在使用 R 中的 mgcv::gam 模型来生成预测，其中时间 (year) 和结果变量 (out) 之间的关系会发生变化。例如，在一种情况下，我想强制时间以线性方式影响结果变量，在另一种情况下以略微减少的方式，在另一种情况下，我想指定时间-结果交互的特定斜率。我不确定如何强制预测以特定方式处理时间和结果变量之间的交互：
res <- gam(out ~ s(time) + s(GEOID, bs='re'), data = df, method = "REML")
pred <- predict(gam, newdata = ndf, type="response", se=T)
【问题讨论】：
标签： r prediction gam
time 和out 之间没有交互；这里time 对out 有潜在的非线性影响。
我们是在谈论尝试为time 的功能强制使用某些形状吗？如果是这样，您将需要估计不同的模型；如果您想要线性效果，请使用time：
res_lin <- gam(out ~ time + s(GEOID, bs='re'), data = df, method = "REML")
并查看形状约束的 p 样条线以强制执行单调性或凹/凸关系。
scam 包具有这些类型的约束，并使用带有 GCV 平滑度选择的 mgcv 来拟合形状约束模型。
至于为时间的线性效应指定特定斜率，我认为您需要在模型中包含time 作为偏移量。所以说你想要的斜率是 0.5 我认为你需要做 + offset(I(0.5*time)) 因为偏移量的系数定义为 1。我会仔细检查一下，因为我可能在这里搞砸了我的想法。
【讨论】：
?smooth.construct.bs.smooth.spec 的帮助，因为您需要使用 b 样条基进行估计，并适当设置惩罚范围以包括预测期。