diff --git a/src/variousRegressionFunctions.jl b/src/variousRegressionFunctions.jl index bb296f7..b46c1d1 100644 --- a/src/variousRegressionFunctions.jl +++ b/src/variousRegressionFunctions.jl @@ -95,24 +95,16 @@ stdX = std(X, dims=1); mX = mean(X, dims=1); my = mean(y, dims=1); y = y .- my; +X = X .- mX; -if centre - X = X .- mX; -end - -#if standardize -# X = X ./ stdX; -#end U, s, V = svd(X, full=false); q = s[1:kmax].^(-1) .*(U[:,1:kmax]'y); B = cumsum(V[:,1:kmax] .* q', dims=2); -if centre - b0 = my .- mX * B - B = [b0; B]; -end +b0 = my .- mX * B +B = [b0; B]; return B, U, s, V end