gam我使用mgcv包中的负二项式系列拟合广义加法模型。我有一个数据框，其中包含我的因变量Y、自变量X、其他自变量Oth和一个因子Fac。我想适合以下模型
Y ~ s(X) + Oth
theta每个因素水平不同。换句话说，我使用
fit <- gam(Y~s(X)+Oth, family=nb())
但这仅给了我theta整个数据集的一个分散参数。相反，我相信各因子的均值是相同的，因此s(X)和只需要一组系数Oth，但方差会随因子而变化，因此我希望theta每个 的水平有一个离散估计Fac。
自然地，每个因子级别拟合一个模型是行不通的，因为我会为每个因子级别的自变量获取一组系数，而不是为整个数据集获取一组系数。