R help - Oct 2005 - GAM and AIC: How can I do??? please

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

   > 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]

======= 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

Please send R and mgcv version numbers, as I can't replicate this
problem. best, Simon

>
>    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
>
>    > 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
>