From 82740963ec33303f14bea07e2df5a85ce5b09a55 Mon Sep 17 00:00:00 2001 From: Joakim Skogholt Date: Sun, 26 May 2024 11:05:16 +0200 Subject: [PATCH] Rettet feil etter forenkling av PCR --- src/variousRegressionFunctions.jl | 14 +++----------- 1 file changed, 3 insertions(+), 11 deletions(-) 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