具有零膨脹數據的 GAMM

是否可以為 R 中的零膨脹數據擬合 GAMM（廣義加法混合模型）？

如果不是，是否可以為 R 中具有負二項式或準泊松分佈的零膨脹數據擬合 GAM（廣義加法模型）？（我發現了用於泊松分佈的COZIGAM::zigam和mgcv:ziP函數）

除了mgcv及其零膨脹泊松族 (ziP()和ziplss())，您還可以查看Paul-Christian Bürkner的 brms包。它可以擬合分佈模型（您不僅可以建模平均值，在您的情況下，可以將模型的零通脹分量建模為協變量的函數，就像計數函數一樣）。

您可以分別通過簡單的 1-d 或各向同性 2-d 樣條或各向異性張量積樣條的術語s()在任何線性預測變量（平均值/計數、零膨脹部分等）中包含平滑。t2()它支持零膨脹二項式、泊松、負二項式和 beta 分佈，以及零一膨脹 beta 分佈。它還具有泊松和負二項式響應的障礙模型（其中模型的計數部分是截斷分佈，以免產生更多的零計數）。

brms使用****STAN擬合這些模型，因此它們是完全貝葉斯的，但這需要您學習一組新的接口來提取相關信息。也就是說，有幾個包只為這個任務提供支持功能，並且brms編寫了利用這些輔助包的輔助函數。您需要安裝 STAN，並且需要 C++ 編譯器，因為brms將使用 R 定義的模型編譯為 STAN 代碼以進行評估。

引用自：https://stats.stackexchange.com/questions/295998