[R] GAM and AIC: How can I do??? please

[ronggui]

======= 2005-10-24 09:55:32 ÄúÔÚÀ´ÐÅÖÐÐ´µÀ£º=======

>
>   Hello,  I'm a Korean researcher who have been started to learn the "R"
>   package.
>
>   I want to make gam model and AIC value of the model to compare several
>   models.
>
>   I did the GAM model, but there were error for AIC.
>
>   SO, how can I do? pleas help me!!!
>
>
>
>   I did like below;
>
>
>   > a.fit <- gam(pi~ s(t1r), family = gaussian(link="log"))
>   > summary(a.fit)
>
>
>   Family: gaussian
>   Link function: log
>
>   Formula:
>   pi ~ s(t1r)
>
>   Parametric coefficients:
>              Estimate  std. err.    t ratio    Pr(>|t|)
>   constant   0.093105   0.005238      17.77    < 2.22e-16
>
>   Approximate significance of smooth terms:
>                 edf       chi.sq     p-value
>   s(t1r)      1.833       24.153     0.00014213
>
>   R-sq.(adj) =  0.435   Deviance explained = 47.1%
>   GCV score = 0.0010938   Scale est. = 0.00099053  n = 30

are you using the mgcv package?
if you are,just use a.fit$aic to get the aic.

>   > AIC(a.fit)
>   Error in logLik(object) : no applicable method for "logLik"
>
>
>   Eun A Kim, MD, MPH, Ph.D
>   Senior Researcher
>   Occupational safety and Health Research Institute
>   Korea Occupational Safety and Health Agency
>   TEL : +82-32-510-0910, FAX: +82-32-518-0862
>   Address:  34-4 Gusan-dong, Bupyung-Gu, Incheon city, 430-711, Republic
>   of Korea
>   Home Fax +82-(303)3111-0573
>   [receive_check.cgi?sender=euna0 at dreamwiz.com&msgid=%3C20051024015532.0
>   000C00B01D902E9 at pmail1.dreamwiz.com%3E&receiver=r-help at stat.math.ethz.
>   ch&key=bd02e77f5ec97f754394e2adff337f11]
>______________________________________________
>R-help at stat.math.ethz.ch mailing list
>https://stat.ethz.ch/mailman/listinfo/r-help
>PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html

= = = = = = = = = = = = = = = = = = = =
			


 

2005-10-24

------
Deparment of Sociology
Fudan University

My new mail addres is ronggui.huang at gmail.com
Blog:http://sociology.yculblog.com