mgcv訪問estimate.gam函數

Question

我正在使用R庫mgcv

gam函數返回具有edf1值的對象

> g = gam(y~s(x, k = 10, bs = "cr"), data = data, scale = 1)
> g$edf1
(Intercept)      s(x).1      s(x).2      s(x).3      s(x).4      s(x).5      s(x).6      s(x).7      s(x).8 
  1.0000000   0.2417857   0.4744928   0.4667241   0.4347516   0.4107624   0.4136064   0.4288544   0.5315970 
     s(x).9 
  0.5941861 
> 

我需要了解如何計算此值，因此我正在閱讀gam()函數的代碼。

在gam()內部，結果是由另一個函數獲得的

object <- estimate.gam(G, method, optimizer, control, in.out, 
                         scale, gamma, ...)

問題是我無法訪問estimate.gam()函數

> estimate.gam
Error: object 'estimate.gam' not found
> 

我該如何解決？

Answer 1

您可以下載源代碼並檢查文件。 功能estimate.gam存儲在該文件將R子目錄mgcv.r 。