R：具有多个负二项式thetas的GAM

Question

我使用gam包中的mgcv在负二项式族中拟合广义加法模型。我有一个数据框，其中包含我的因变量Y，一个独立变量X，其他独立变量Oth和一个因子Fac。我想适合以下模型

Y ~ s(X) + Oth

每个因素级别不同theta。换句话说，我使用

fit <- gam(Y~s(X)+Oth, family=nb())

但这只为我提供了整个数据集的一个分散参数theta。相反，我认为平均值在各个因素之间是相同的，因此s(X)和Oth只需要一组系数，但方差会因子而变化，所以我想要一个色散估计值{{ 1}}每级theta。

当然，每个因子水平拟合一个模型是行不通的，因为我会得到每个因子水平的自变量的一组系数，而不是整个数据集的一个系数。

Answer 1

解决问题的最佳方法是使用gamlss包。您可以按mean变量对X进行建模，variance和X变量对Fac进行建模，这样您就可以获得所有levels变量的参数你的因素。