> I think I've found a bug in mfittrend in the version of kappa
> (1.9-11) which was part of the lehuakona starlink release. It seems
> that straight line fitting can give very wrong fits for certain map
> regions in this release but works in the previous humu and, interestingly
> the latest nanahope release.
> This bug could potentially have serious repercussions for any NGLS
> data reduced in the 9 month time period where lehuakona was the main
> release.
Whatever it was has been fixed. Are you asking that the offending
bug be relocated and a patched lehuakona be made? Off the top of my
head I don't know what bug was introduced and subsequently fixed.
I'd seen something like this with higher-order polynomials behaving
badly at the spectrum ends, but not as extreme, when trying to automate
the masking and fitting for some of Remo's galaxies. Thus I recommended
visual checking and why there was a warning message. You could also
test for anomalies in the residuals.
> Details for reproduction are below.
I actually found it easier to output to B_smooth_bsl and B_smooth_msk
for the baseline and mask NDFs, and use DATACUBE/trendview, for example.
% trendview -i B_smooth"(-5:3,7:12,)"
to check that the fits are suitable. The fit lines should approximately
cover the original spectra apart from noise and features. If you can
see the yellow spectrum and separate green fit curve, something has gone
wrong.
> I do not know if there is any platform dependence - the test machine
> for the buggy lehuakona results was an Ubuntu Linux Hardy Heron machine.
lehuakona reproduced the problem on 64-bit OpenSuSE 11.1. It appears
mostly frequently in the largest and smallest but it's not a rule. I
can see well-fit low- and high-signal spectra.
Malcolm
|