re-instate lag fdHess sse

Author: rsbivand

inst/tinytest/test_ml_sse_hess.R

@@ -5,11 +5,11 @@ boston_nb <- spdep::poly2nb(boston.tr)
 lw <- spdep::nb2listw(boston_nb)
 e <- eigenw(lw)
 f <- log(MEDV) ~ CRIM + ZN + INDUS + CHAS + I(NOX^2) + I(RM^2) + AGE + log(DIS) + log(RAD) + TAX + PTRATIO + B + log(LSTAT)
-error.eig <- errorsarlm(f, boston.tr, lw, control=list(pre_eig=e, fdHess=TRUE))
+error.eig <- errorsarlm(f, boston.tr, lw, control=list(pre_eig=e, compiled_sse=FALSE, fdHess=TRUE))
 error.eig.alt <- errorsarlm(f, boston.tr, lw, control=list(pre_eig=e, compiled_sse=TRUE, fdHess=TRUE))
 expect_true(isTRUE(all.equal(error.eig$lambda, error.eig.alt$lambda, scale=1, tolerance=3e-8)))
 expect_true(isTRUE(all.equal(error.eig$fdHess, error.eig.alt$fdHess, scale=1, tolerance=1e-5)))
-# lag.eig <- lagsarlm(f, boston.tr, lw, control=list(pre_eig=e, fdHess=TRUE))
-# lag.eig.alt <- lagsarlm(f, boston.tr, lw, control=list(pre_eig=e, compiled_sse=TRUE, fdHess=TRUE))
-# expect_true(isTRUE(all.equal(lag.eig$fdHess, lag.eig.alt$fdHess, scale=1, tolerance=1e-6)))
+lag.eig <- lagsarlm(f, boston.tr, lw, control=list(pre_eig=e, compiled_sse=FALSE, fdHess=TRUE))
+lag.eig.alt <- lagsarlm(f, boston.tr, lw, control=list(pre_eig=e, compiled_sse=TRUE, fdHess=TRUE))
+expect_true(isTRUE(all.equal(lag.eig$fdHess, lag.eig.alt$fdHess, scale=1, tolerance=1e-6)))