R 提取GAM模型对象

假设我通过以下操作创建GAM模型：

a <- runif(10)
b <- runif(10)
gm <- gam(a ~ ns(b, df=2))
plot(gm, all.terms=T, shade=T)

及

检查（而不是使用

摘要

）-它提供对象的结构

我认为

gm$model

正是您想要的

gm$model
            a ns(b, df = 2).1 ns(b, df = 2).2
1  0.69342149      0.07841860     -0.05184526
2  0.23538533      0.52006793      0.20238728
3  0.47125666      0.24808303     -0.15840080
4  0.04405890      0.00000000      0.00000000
5  0.54696387      0.34211652      0.77302788

---

方向很好。但我仍然无法重现以红色突出显示的功能

> names(gm)
 [1] "coefficients"      "residuals"         "fitted.values"     "family"            "linear.predictors"
 [6] "deviance"          "null.deviance"     "iter"              "weights"           "prior.weights"    
[11] "df.null"           "y"                 "converged"         "sig2"              "edf"              
[16] "edf1"              "hat"               "R"                 "boundary"          "sp"               
[21] "nsdf"              "Ve"                "Vp"                "rV"                "mgcv.conv"        
[26] "gcv.ubre"          "aic"               "rank"              "gcv.ubre.dev"      "scale.estimated"  
[31] "method"            "smooth"            "formula"           "var.summary"       "cmX"              
[36] "model"             "control"           "terms"             "pred.formula"      "pterms"           
[41] "assign"            "xlevels"           "offset"            "df.residual"       "min.edf"          
[46] "optimizer"         "call" 

gm$model
            a ns(b, df = 2).1 ns(b, df = 2).2
1  0.69342149      0.07841860     -0.05184526
2  0.23538533      0.52006793      0.20238728
3  0.47125666      0.24808303     -0.15840080
4  0.04405890      0.00000000      0.00000000
5  0.54696387      0.34211652      0.77302788