當前位置: 首頁>>代碼示例 >>用法及示例精選 >>正文


R fix.family.link 修改族以用於 GAM 擬合和檢查


R語言 fix.family.link 位於 mgcv 包(package)。

說明

通過 ‘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 估計需要對數飽和似然及其前兩個導數。尺度參數。這些函數將評估這些量的函數添加到一個族中。

如果該係列已經具有函數 dvard2vard3vard2linkd3linkd4link 以及 RE/ML ls ,則這些函數僅返回未修改的係列:這允許非- 使用外部迭代時與 gam 一起使用的標準鏈接(性能迭代與未修改的族一起運行)。請注意,如果您隻需要 Fisher 評分,則 d4linkd3var 可以是虛擬的,因為它們會被忽略。類似的 ls 僅適用於 RE/ML。

dvar 函數是均值向量 mu 的函數,並返回族方差函數的相應一階導數的向量。 d2link 函數也是平均值向量 mu 的函數:它返回鏈接的二階導數向量,在 mu 處求值。高階導數的定義類似。

如果修改您自己的係列,請注意,如果您的係列僅需要標準係列之一中出現的鏈接,則通常可以隻提供 dvard2var 函數。

後兩個函數對於研究殘差的分布很有用,由 qq.gam 使用。如果可能,函數將分位數 (qf) 或隨機偏差 (rd) 生成函數添加到族中。如果一個係列已經具有 qfrd 函數,則保持不變。 qf 函數僅適用於某些係列,而對於準係列,兩種類型的函數均不可用。

具有額外組件函數 dvard2vard2linkd3linkd4linkls 以及可能 qfrd 的係列對象,具體取決於調用的函數。 fix.family.var 還添加了一個設置為負數的變量 scale,以指示該係列具有自由縮放參數。

作者

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

也可以看看

gam.fit3 , qq.gam

相關用法


注:本文由純淨天空篩選整理自R-devel大神的英文原創作品 Modify families for use in GAM fitting and checking。非經特殊聲明,原始代碼版權歸原作者所有,本譯文未經允許或授權,請勿轉載或複製。