Skip to content

GAMs under R?

2 messages · Bill Venables, Brian Ripley

#
At 06:09 AM 5/1/00 +0100, Prof Brian D Ripley wrote:
The restriction "for least-squares fitting" is crucial.  bruto is a form of
additive model fitting, not *generalized* additive model fitting.  If you
have to deal with binomial, poisson, .... models it doesn't help much.
But the effort required should not be under-estimated and the need for a
*good* implementation is pretty strict.  A mediocre or worse implementation
is in many ways worse than useless and even dangerous.  This looks like a
pretty interesting component of an Honours project in Statistics, if some
student is interested in computations and penalized estimation.
.-.-
._._
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-help-request at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
#
On Mon, 1 May 2000, Bill Venables wrote:

            
Whoah!  It helps a great deal.  gam() fits by IWLS just like glm(), so you
can build an outer loop around bruto(), say.  The variable selection
aspects of bruto make that somewhat unstable, but with a good starting
set of weights (e.g. from a large linear fit) it can work well.  (I did
some experiments in ca 1994 when writing my PRNN book.)
(And where do you place the one in S-PLUS?  The main difficulty with
implementing modreg was bugs in the underlying code, bugs which also exist
in the S-PLUS implementation, but Linux compilers/glibc seemed to 
stumble over them where Solaris did not.)

Those who look at the future plans for R (on developer.r-project.org)
will see in my section:

Longer term

...

Something better than gam for additive models. Perhaps along the lines of
BRUTO (in package mda).

[That's not likely to materialize for some months, though.]