[R] help interpreting output?

Michael F. Palopoli mpalopol at bowdoin.edu
Tue Jan 7 02:52:03 CET 2003


Dear R experts,

I'm hoping someone can help me to interpret the results of building 
gam's with mgcv in R.

Below are summaries of two gam's based on the same dataset.  The first 
gam (named "gam.mod") has six predictor variables.  The second gam 
(named "gam.mod2") is exactly the same except it is missing one of the 
predictor variables.  What is confusing me is the estimated defrees of 
freedom for each of the splines in the second model....

________________

 > summary.gam(mod.gam)

Family: gaussian
Link function: identity

Formula:
INT ~ s(IGS) + s(L2E) + s(TED) + s(PSD) + s(OPD) + s(GED)

Parametric coefficients:
           Estimate  std. err.    t ratio    Pr(>|t|)
constant     302.32      5.192      58.23    < 2.22e-16

Approximate significance of smooth terms:
              edf       chi.sq     p-value
s(IGS)      4.254       58.308     9.5524e-12
s(L2E)          1       8.7673     0.0030668
s(TED)          1       8.3915     0.0037697
s(PSD)          1       6.0234     0.014118
s(OPD)      2.289       12.745     0.0024349
s(GED)      3.791       152.68     < 2.22e-16

R-sq.(adj) = 0.885   Deviance explained = 91.1%
GCV score = 2124.9   Scale est. = 1617.3    n = 60

________________

 >summary.gam(mod.gam2)

Family: gaussian
Link function: identity

Formula:
INT ~ s(IGS) + s(L2E) + s(TED) + s(PSD) + s(OPD)

Parametric coefficients:
           Estimate  std. err.    t ratio    Pr(>|t|)
constant     302.32  4.736e-14  6.384e+15    < 2.22e-16

Approximate significance of smooth terms:
              edf       chi.sq     p-value
s(IGS)  1.757e-05   1.3524e+09     < 2.22e-16
s(L2E)   0.009991      0.21394     0.6437
s(TED)  2.945e-05   1.4913e+07     < 2.22e-16
s(PSD)  2.566e-05   6.5495e+06     < 2.22e-16
s(OPD)  5.023e-05   3.2332e+07     < 2.22e-16

R-sq.(adj) = 0.645   Deviance explained = 64.5%
GCV score = 7489.7   Scale est. = 6069.7    n = 60


________________


Any suggestions about either (1) what went wrong with the second model? 
 or (2) how the heck do I interpet these results?

Thanks,

Mike.




More information about the R-help mailing list