From P.Dalgaard at biostat.ku.dk  Fri Mar  6 15:57:22 2009
From: P.Dalgaard at biostat.ku.dk (Peter Dalgaard)
Date: Fri, 06 Mar 2009 15:57:22 +0100
Subject: [Rd] bug (PR#13570)
In-Reply-To: <alpine.LFD.2.00.0903060718570.24770@auk.stats.ox.ac.uk>
References: <20090304211010.A588D2832180@mail.pubhealth.ku.dk>	<49AF88B7.1000002@biostat.ku.dk>	<20090305174301.20099b11@berwin-nus1>	<49AFB412.2060807@statistik.tu-dortmund.de>	<alpine.LFD.2.00.0903051121300.2646@gannet.stats.ox.ac.uk>	<49AFC711.1030902@biostat.ku.dk>	<alpine.LFD.2.00.0903051237560.3864@gannet.stats.ox.ac.uk>	<E0F590EA-8F52-44A6-AD5D-E050BD6740D7@stat.purdue.edu>	<B37C0A15B8FB3C468B5BC7EBC7DA14CC61CC034A04@LP-EXMBVS10.CO.IHC.COM>	<7A1AA352-D2E8-48A2-BDE3-FF48B71E28D9@stat.purdue.edu>	<49B07B45.3040308@gmail.com>
	<alpine.LFD.2.00.0903060718570.24770@auk.stats.ox.ac.uk>
Message-ID: <49B139D2.5060101@biostat.ku.dk>

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



[Rd] bug (PR#13570)
Peter Dalgaard