Re: [Rd] LOESS (PR#9064)

From: <>
Date: Fri 07 Jul 2006 - 09:39:34 GMT

On Fri, 7 Jul 2006, wrote:

> Hello,
> I found a little BUG in loess <stats>. It does not receive the iterations
> parameter.

Please provide an example demonstrating this 'BUG', as we do ask you to do, along with basic information like your name, the version of R, the platform ....

> It can be debugged in the following way:

(Not in any sense of `debugged' that my dictionary mentions.)

> ....
> fit <- simpleLoess(y, x, w, span, degree, parametric, drop.square,
> normalize, control$statistics, control$surface, control$cell,
> iterations, control$trace.hat)
> Replace argument iterations with control$iterations.

Why? The code contains a few lines earlier

     iterations <- if(family=="gaussian") 1 else control$iterations

which is the intention. ?loess.control does say

iterations: the number of iterations used in robust fitting.

> I found also remarkable differences when comparing loess with PROC LOESS
> (in SAS) using the same parameters and variables. It really confuses me.
> Where can I report this?

R-bugs is not a place to ask questions, but bugs in SAS can be reported to SAS inc.

Brian D. Ripley,        
Professor of Applied Statistics,
University of Oxford,             Tel:  +44 1865 272861 (self)
1 South Parks Road,                     +44 1865 272866 (PA)
Oxford OX1 3TG, UK                Fax:  +44 1865 272595

______________________________________________ mailing list
Received on Fri Jul 07 19:43:01 2006

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.1.8, at Fri 07 Jul 2006 - 10:25:42 GMT.

Mailing list information is available at Please read the posting guide before posting to the list.