[R] AIC for robust GAM ?

[Bert Gunter]

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