1

我正在使用 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
> 

我该如何解决这个问题?

4

1 回答 1

1

您可以下载源代码并检查文件。该函数estimate.gam存储在文件的 R 子目录中mgcv.r

于 2015-07-03T11:14:14.340 回答