R fix.family.link 修改族以用于 GAM 拟合和检查

说明

通过 ‘outer’ 迭代进行广义加法模型拟合，需要将方差和链接函数的额外导数添加到族对象中。前 3 个函数添加了所需的内容。可以通过向族中添加分位数和随机偏差生成函数来辅助模型检查。最后两个函数执行此操作。

用法

fix.family.link(fam)
fix.family.var(fam)
fix.family.ls(fam)
fix.family.qf(fam)
fix.family.rd(fam)

参数

fam

一个family。

细节

首先考虑前 3 个函数。

外迭代 GAM 估计需要 GCV、UBRE/gAIC、GACV、REML 或 ML 分数的导数，这些导数是通过查找模型系数 w.r.t. 的导数获得的。使用隐函数定理对数平滑参数。导数的表达式需要链接 w.r.t 的二阶和三阶导数。平均值(如果不使用 Fisher 评分，则为四阶导数)。还需要方差函数 w.r.t 的一阶和二阶导数。平均值(如果不使用 Fisher 评分，则加上三阶导数)。最后，平滑参数的 REML 或 ML 估计需要对数饱和似然及其前两个导数。尺度参数。这些函数将评估这些量的函数添加到一个族中。

如果该系列已经具有函数 dvar 、 d2var 、 d3var 、 d2link 、 d3link 、 d4link 以及 RE/ML ls ，则这些函数仅返回未修改的系列：这允许非- 使用外部迭代时与 gam 一起使用的标准链接(性能迭代与未修改的族一起运行)。请注意，如果您只需要 Fisher 评分，则 d4link 和 d3var 可以是虚拟的，因为它们会被忽略。类似的 ls 仅适用于 RE/ML。

dvar 函数是均值向量 mu 的函数，并返回族方差函数的相应一阶导数的向量。 d2link 函数也是平均值向量 mu 的函数：它返回链接的二阶导数向量，在 mu 处求值。高阶导数的定义类似。

如果修改您自己的系列，请注意，如果您的系列仅需要标准系列之一中出现的链接，则通常可以只提供 dvar 和 d2var 函数。

后两个函数对于研究残差的分布很有用，由 qq.gam 使用。如果可能，函数将分位数 (qf) 或随机偏差 (rd) 生成函数添加到族中。如果一个系列已经具有 qf 或 rd 函数，则保持不变。 qf 函数仅适用于某些系列，而对于准系列，两种类型的函数均不可用。

值

具有额外组件函数 dvar 、 d2var 、 d2link 、 d3link 、 d4link 、 ls 以及可能 qf 和 rd 的系列对象，具体取决于调用的函数。 fix.family.var 还添加了一个设置为负数的变量 scale，以指示该系列具有自由缩放参数。

作者

Simon N. Wood simon.wood@r-project.org

也可以看看

gam.fit3 , qq.gam