Included constant term in oldregcoeffs

This commit is contained in:
Joakim Skogholt 2024-05-19 12:43:36 +02:00
parent 79657c58e0
commit c92b9c1186

View file

@ -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 regCoeffsInds = findfirst(cverror .< (minimum(cverror) + std(cverror)/sqrt(length(unique(cvfolds)))))[1] # 1 S.E. rule
betas, _ = regFunction(X, y, ncomps) betas, _ = regFunction(X, y, ncomps)
bold = betas[2:end, regCoeffsInds] bold = betas[:, regCoeffsInds]
return bold, regCoeffsInds, cverror return bold, regCoeffsInds, cverror
end end