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.dalga...@biostat.ku.dk) FAX: (+45) 35327907 ______________________________________________ R-devel@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-devel
