Re: [R] "gam()" in "gam" package

From: Kunio takezawa <r.otasuke_at_gmail.com>
Date: Wed, 19 Dec 2007 10:10:44 +0900

R-users
E-mail: r-help_at_r-project.org

> This iteration seems to be for "iteratively reweighted least squares" not
>for backfitting. And lm.wfit may solve multiple linear equation using
>QR decomposition; but I am not sure.

   Let me tell you something about my guess above. The iteration below is from 1 to maxit.

  for (iter in 1:maxit) {


        fit <- eval(bf.call)
        ------

  }

The help of "gam.control()" says:

   maxit: maximum number of local scoring iterations    bf.maxit: maximum number of backfitting iterations

Therefore the iteration above is local scoring iteration not backfitting iteration.

   The first line of "lm.wfit()" is:
function (x, y, w, offset = NULL, method = "qr", tol = 1e-07,  singular.ok = TRUE, ...)
There is no "bf.maxit" or its equivalent in these arguments. And it contains 'method = "qr"'; it may mean that this routine solves a multiple linear equation using QR decomposition.

-- 
*****    r.otasuke_at_gmail.com    *****
http://cse.naro.affrc.go.jp/takezawa/intro.html

	[[alternative HTML version deleted]]

______________________________________________
R-help_at_r-project.org mailing list
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.
Received on Wed 19 Dec 2007 - 01:14:10 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Thu 20 Dec 2007 - 01:30:20 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.