diff --git a/src/TR.jl b/src/TR.jl
index 02d1cdd..951fa2d 100644
--- a/src/TR.jl
+++ b/src/TR.jl
@@ -131,7 +131,7 @@ end
 # Calculating rmsecv and regression coefficients
 press   = sum(rescv.^2, dims=1)';
 rmsecv  = sqrt.(1/n .* press);
-bcoeffs      = V * broadcast(./, (U' * ys), denom) .+ bold .- V * broadcast(./, V' * bold, denom2);
+bcoeffs      = V * broadcast(./, (U' * y), denom) .+ bold .- V * broadcast(./, V' * bold, denom2);
 bcoeffs = regMat \ bcoeffs;
 
 # Creating regression coefficients for uncentred data
@@ -209,7 +209,7 @@ lambdarmsecv = lambdas[idminrmsecv];
 lambdagcv    = lambdas[idmingcv];
     
 # Calculating regression coefficients
-bcoeffs      = V * broadcast(./, (U' * ys), denom) .+ bold .- V * broadcast(./, V' * bold, denom2);   
+bcoeffs      = V * broadcast(./, (U' * y), denom) .+ bold .- V * broadcast(./, V' * bold, denom2);   
 bcoeffs      = regMat \ bcoeffs;
     
 if my != 0