我正在拟合 logit GAM 模型,需要对平滑项中的 beta 施加不等式约束。我已经尝试过 mgcv::pcls()
,但似乎我无法指定不同样条系数的独立不等式。实现这一目标的推荐方法是什么?
最佳答案
我自己找到了答案,发布它以防有人会寻找它。 mgcv
中的 mono.con
函数生成约束矩阵 Ain 和 bin。对于多个样条,Ain 矩阵应该放在一起形成一个 block 对角矩阵。
关于r - 问:mgcv::gam() 对系数有限制,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/58186927/