diff --git a/src/variousRegressionFunctions.jl b/src/variousRegressionFunctions.jl
index feb0627..81cc4b2 100644
--- a/src/variousRegressionFunctions.jl
+++ b/src/variousRegressionFunctions.jl
@@ -26,7 +26,7 @@ function oldRegCoeffs(X, y, cvfolds, ncomps, regFunction=bidiag2) # bidiag2 OR P
     regCoeffsInds = findfirst(cverror .< (minimum(cverror) + std(cverror)/sqrt(length(unique(cvfolds)))))[1] # 1 S.E. rule
     
     betas, _ = regFunction(X, y, ncomps)
-    bold = betas[2:end, regCoeffsInds]
+    bold = betas[:, regCoeffsInds]
     
     return bold, regCoeffsInds, cverror
 end