diff --git a/src/TR.jl b/src/TR.jl index 6a321a3..1487c91 100644 --- a/src/TR.jl +++ b/src/TR.jl @@ -133,7 +133,7 @@ function TRPress(TR::TRSVD, y, lambdas) denom = broadcast(.+, broadcast(./, lambdas, TR.s'), TR.s')'; resid = broadcast(.-, y, TR.U * broadcast(./, TR.s .* (TR.U'*y), denom)); H = broadcast(.+, TR.U.^2 * broadcast(./, TR.s, denom), 1/TR.n); - @inline press = TRPress(TR, y, lambdas, denom, H, resid); + @inline press = TRPress(TR, y, lambdas, H, resid); return press end @@ -158,7 +158,7 @@ function TRGCV(TR::TRSVD, y, lambdas) resid = broadcast(.-, y, TR.U * broadcast(./, TR.s .* (TR.U'*y), denom)); H = broadcast(.+, TR.U.^2 * broadcast(./, TR.s, denom), 1/TR.n); - @inline GCV = TRGCV(TR, y, lambdas, denom, H, resid); + @inline GCV = TRGCV(TR, y, lambdas, H, resid); return GCV end