From ripley at stats.ox.ac.uk  Fri Jul  7 13:08:09 2006
From: ripley at stats.ox.ac.uk (ripley at stats.ox.ac.uk)
Date: Fri,  7 Jul 2006 13:08:09 +0200 (CEST)
Subject: [Rd] User Error (was LOESS (PR#9064))
Message-ID: <20060707110809.4B0AA549E9@slim.kubism.ku.dk>

Please do as we ask (repeatedly) and study the help page before posting.

'family' is a separate argument, not part of loess.control, as the help 
page correctly documents.  If you use

cars.lo2 <- loess(dist ~ speed, cars, family = "symmetric",
  control = loess.control(surface = "direct",  iterations = 20))
cars.lo2$pars$iterations

it prints *20*, as it is documented to do.


On Fri, 7 Jul 2006, lamp at few.eur.nl wrote:

> Dear Mr Ripley,

> Thank you for your help.
>
> Version of R: 2.3.1
> Platform: Windows XP
>
> Example:
>
> cars.lo2 <- loess(dist ~ speed, cars,
>  control = loess.control(surface = "direct" , family = 'symmetric' ,
> iterations = 20))
>  cars.lo2$pars$iterations
>
> Returns 1, where it should return 10.

Why 10? It should and does return 1, but you seem to think you asked for 
20!

> If one makes the alterations this will not happen.

It would then be incorrect on the example on the help page.

>
> Best regards,
>
> Felix Lamp
>
>
>
>
>> On Fri, 7 Jul 2006, lamp at few.eur.nl 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.)
>>
>>> THIS IS AN EXCERPT FROM THE CODE:
>>>
>>> ....
>>>
>>> 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,                  ripley at stats.ox.ac.uk
>> Professor of Applied Statistics,  http://www.stats.ox.ac.uk/~ripley/
>> University of Oxford,             Tel:  +44 1865 272861 (self)
>> 1 South Parks Road,                     +44 1865 272866 (PA)
>> Oxford OX1 3TG, UK                Fax:  +44 1865 272595
>>
>
>
>

-- 
Brian D. Ripley,                  ripley at stats.ox.ac.uk
Professor of Applied Statistics,  http://www.stats.ox.ac.uk/~ripley/
University of Oxford,             Tel:  +44 1865 272861 (self)
1 South Parks Road,                     +44 1865 272866 (PA)
Oxford OX1 3TG, UK                Fax:  +44 1865 272595



[Rd] User Error (was LOESS (PR#9064))
ripley at stats.ox.ac.uk