問：mgcv::gam() 對系數有限制

Question

我正在擬合 logit GAM model 並且需要在平滑條件下對 beta 施加不等式約束。 我試過 mgcv::pcls() 但似乎我無法為不同的樣條系數指定獨立的不等式。 實現這一目標的推薦方法是什么？

Answer 1

自己找到了答案，發布它以防有人會尋找它。 mgcv 中的 Mono.con function 生成約束矩陣 Ain 和 bin。 對於幾個樣條，Ain 矩陣應該只是放在一起成一個塊對角矩陣。