"M. Edward (Ed) Borasky" <[EMAIL PROTECTED]> writes: [snip] > >Hmm, could you replace the a1 == a2 with all.equal(a1, a2) instead? > >(inside reg-tests-1.R of course) > > > >Asking for identity up to machine precision does look a bit optimistic... > > > > > That worked ... it got through reg-tests-1.R fine. However, it failed > a little further down in the NA handling tests: > > running tests of NA handling functions > make[3]: Entering directory `/home/znmeb/R-beta/tests' > running code in 'nafns.R' ...make[3]: *** [nafns.Rout] Error 1 > make[3]: Leaving directory `/home/znmeb/R-beta/tests' > make[2]: *** [test-Nafns] Error 2 > make[2]: Leaving directory `/home/znmeb/R-beta/tests' > make[1]: *** [test-all-devel] Error 1 > make[1]: Leaving directory `/home/znmeb/R-beta/tests' > make: *** [check-all] Error 2 > > > sm(fitted(gfit2)) > length 153 with 42 NAs > > sm(resid(gfit2)) > length 153 with 42 NAs > > sm(predict(gfit2)) > length 153 with 42 NAs > > (pp2 <- predict(gfit2, nd)) > 6 25 26 27 > NA -16.177404 1.688479 NA > > stopifnot(all.equal(pp, pp2)) > > > > ## more precise tests. > > f1 <- fitted(gfit) > > f2 <- fitted(gfit2) > > common <- match(names(f1), names(f2)) > > stopifnot(max(abs(f1 - f2[common])) < 100*.Machine$double.eps) > Error: max(abs(f1 - f2[common])) < 100 * .Machine$double.eps is not TRUE > Execution halted > "tests/nafns.Rout.fail"
This looks more serious. 100 times machine precision is quite a large margin in these matters. Could you perhaps stick in a printout of the two terms and their difference? I have an ATLAS build on AMD64 and it passes all the checks, but it is using ATLAS 3.7.8, so you might want to try an upgrade. -- O__ ---- Peter Dalgaard Blegdamsvej 3 c/ /'_ --- Dept. of Biostatistics 2200 Cph. N (*) \(*) -- University of Copenhagen Denmark Ph: (+45) 35327918 ~~~~~~~~~~ - ([EMAIL PROTECTED]) FAX: (+45) 35327907 ______________________________________________ R-devel@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-devel
