diff --git a/src/TR.jl b/src/TR.jl index b9a924f..e58c21a 100644 --- a/src/TR.jl +++ b/src/TR.jl @@ -19,8 +19,8 @@ for i = 1:n p2 = size(Xdata, 2); for j = 1:length(lambdasu) - betas = [Xs; sqrt(lambdasu[j]) * I(p2)] \ [ys ; sqrt(lambdasu[j]) * bOld]; - rmsecvman[j] += (y[i] - (((X[i,:]' .- mX) * betas)[1] + my))^2; + betas = [Xs; sqrt(lambdasu[j]) * I(p2)] \ [ys ; sqrt(lambdasu[j]) * bOld]; + rmsecv[j] += (y[i] - (((X[i,:]' .- mX) * betas)[1] + my))^2; end end