[R] Extract components of gam object
    arun 
    smartpink111 at yahoo.com
       
    Sat Sep  7 17:54:54 CEST 2013
    
    
  
Hi,
summary(pres.gam)
Family: gaussian 
Link function: identity 
Formula:
prestige ~ s(income) + s(education)
Parametric coefficients:
            Estimate Std. Error t value Pr(>|t|)    
(Intercept)  47.3276     0.6914   68.45   <2e-16 ***
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Approximate significance of smooth terms:
               edf Ref.df     F  p-value    
s(income)    3.079  3.813 12.18 8.03e-08 ***
s(education) 3.005  3.746 42.68  < 2e-16 ***
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
R-sq.(adj) =   0.84   Deviance explained =   85%
GCV score = 50.493  Scale est. = 46.843    n = 98
Assuming that you want to extract `3.079`.
summary(pres.gam)$edf[1]
#[1] 3.07938
A.K.
I want to extract the components of a gam object. Sample code and data follows: 
library(MASS) 
library(mgcv)#package for GAM 
library(car)#to get the data 
data(Prestige) 
Prestige2<-na.omit(Prestige) 
pres.gam<-gam(prestige~s(income)+s(education),data =Prestige2) 
As an example I wish to extract the edf for income. In addition I
 want to learn how to systematically understand the components of 
structure pres.gam. I have tried str(pres.gam) which produce large 
output which is difficult to understand. Please explain using the edf 
component. 
Thank you
    
    
More information about the R-help
mailing list