Per the posting guide linked below:
"If the question relates to a *contributed package* , e.g., one downloaded
from CRAN, try contacting the package maintainer first. You can also use
find("functionname") and packageDescription("packagename")
to find this
information. *Only* send such questions to R-help or R-devel if you get no
reply or need further assistance. This applies to both requests for help
and to bug reports."
So I believe that you should try to contact the maintainer for your
question first. You may get lucky here, of course, but your query seems
rather too technical to expect a useful response on R-Help.
-- Bert
Bert Gunter
"The trouble with having an open mind is that people keep coming along and
sticking things into it."
-- Opus (aka Berkeley Breathed in his "Bloom County" comic strip )
On Wed, Jan 6, 2021 at 12:17 PM varin sacha via R-help <r-help at
r-project.org>
wrote:
> Dear R-Experts,
>
> Here below my reproducible R code.
> How can I get the AIC of my model (robust GAM) ?
>
> Best Regards,
>
>
>
>
y<-c(499,491,500,517,438,495,501,525,516,494,500,453,479,481,505,465,477,520,520,480,477,416,502,503,497,513,492,469,504,482,502,498,463,504,495)
>
>
x<-c(499,496,424,537,480,484,503,575,540,436,486,506,496,481,508,425,501,519,546,507,452,498,471,495,499,522,509,474,502,534,504,466,527,485,525)
> library(robustgam)
> true.family <- poisson()
> fit=robustgam(x,y,sp=0,family=true.family,smooth.basis='ps',K=3)
> AIC(fit)
>
>
>
>
>
> ______________________________________________
> R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide
> http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
>
[[alternative HTML version deleted]]