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

Question

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

Y ~ s(X) + Oth

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

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

但这仅给我一个分散参数theta 用于整个数据集。相反，我相信各因子的平均值是相同的，因此s(X) 和 Oth 只需要一组系数，但方差会随着因子的变化而变化，所以我希望每个水平的离散度估计 theta Fac.

当然，每个因子级别拟合一个模型是行不通的，因为我会为每个因子级别的自变量获取一组系数，而不是为整个数据集获取一组系数。

Answer 1

解决问题的最佳方法是使用gamlss 包。您将能够通过X 变量和variance 通过X 和Fac 变量对mean 建模，因此您将获得所有levels 因子的参数。