[Rd] bug (PR#13570)

[Peter Dalgaard]

Prof Brian Ripley wrote:
> On Thu, 5 Mar 2009, Benjamin Tyner wrote:
[...]
>>
>> I submitted a bug fix to Eric Grosse, the maintainer of the netlib
>> routines; the fixed lines of fortran are identified in the comments at
>> (just search for my email address):
>>
>> http://www.netlib.org/a/loess
>>
>> These fixes would be relatively simple to incorporate into R's version
>> of loessf.f
> 
> The fixes from dloess even more simply, since R's code is based on
> dloess.  Thank you for the suggestion.
> 
> Given how tricky this is to reproduce, I went back to my example under
> valgrind.  If I use the latest dloess code, it crashes, but by
> selectively importing some of the differences I can get it to work.
> 
> So it looks as if we are on the road to a solution, but something in the
> current version (not necessarily in these changes) is incompatible with
> the current R code and I need to dig further (not for a few days).

What a nice "war story" this is!

Good that it now seems fixable; even though degree=0 is not of much
practical use, it is the sort of thing people like to have available
when explaining how the method works.


-- 
   O__  ---- Peter Dalgaard             Øster Farimagsgade 5, Entr.B
  c/ /'_ --- Dept. of Biostatistics     PO Box 2099, 1014 Cph. K
 (*) \(*) -- University of Copenhagen   Denmark      Ph:  (+45) 35327918
~~~~~~~~~~ - (p.dalgaard at biostat.ku.dk)              FAX: (+45) 35327907