diff --git a/src/TR.jl b/src/TR.jl
index 33a4db7..2a7d869 100644
--- a/src/TR.jl
+++ b/src/TR.jl
@@ -83,7 +83,7 @@ lambda_min, lambda_min_ind = findmin(rmsecv)
 lambda_min_ind             = lambda_min_ind[1]
 b_lambda_min               = bcoeffs[:,lambda_min_ind]
     
-return b_lambda_min, rmsecv, lambda_min, lambda_min_ind
+return b_lambda_min, rmsecv, lambda_min, lambda_min_ind, bcoeffs
 end
 
 """
@@ -159,7 +159,7 @@ brmsecv = bcoeffs[:, idminrmsecv];
 bgcv    = bcoeffs[:, idmingcv];
 
 
-return brmsecv, bgcv, rmsecv, gcv, idminrmsecv, lambdarmsecv, idmingcv, lambdagcv
+return brmsecv, bgcv, rmsecv, gcv, idminrmsecv, lambdarmsecv, idmingcv, lambdagcv, bcoeffs
 end
 
 
@@ -288,7 +288,7 @@ end
 bpress = bcoeffs[:, idminPRESS]
 bgcv   = bcoeffs[:, idminGCV]
 
-return bpress, bgcv, rmsecv, GCV, idminPRESS, lambdaPRESS, idminGCV, lambdaGCV
+return bpress, bgcv, rmsecv, GCV, idminPRESS, lambdaPRESS, idminGCV, lambdaGCV, bcoeffs
 end
 
 """
@@ -372,7 +372,7 @@ lambda_min, lambda_min_ind = findmin(rmsecv)
 lambda_min_ind             = lambda_min_ind[1]
 b_lambda_min               = bcoeffs[:,lambda_min_ind]
     
-return b_lambda_min, rmsecv, lambda_min, lambda_min_ind
+return b_lambda_min, rmsecv, lambda_min, lambda_min_ind, bcoeffs
 end
 
 """