Added EMSC

src/Ting.jl

@@ -30,6 +30,8 @@ export TRLooCVUpdateNaive
 export TRLooCVUpdateExperimental
 
 # From "variousRegressionFunctions.jl"
+export EMSC
+export EMSCCorrection
 export cVals
 export calculateRMSE
 export oldRegCoeffs

src/variousRegressionFunctions.jl

@@ -1,3 +1,24 @@
+function EMSCCorrection(X, basis)
+    
+    coeffs   = basis \ X';
+    X_Cor   = (X - coeffs[2:end,:]' * basis[:,2:end]') ./ coeffs[1,:];
+    
+    return X_Cor
+end
+
+function EMSC(X, polDeg=2)
+
+n, p    = size(X);
+P       = zeros(p, polDeg+1);
+[P[:,i] = LinRange(-1, 1, p).^(i-1) for i in 1:polDeg+1];
+refSpec = mean(X, dims=1)';
+basis   = [refSpec P];
+
+X_Cor   = EMSCCorrection(X, basis);
+    
+return X_Cor, basis
+end
+
 """
 crit1minInd, crit2minInd, CvalCrit1, CvalCrit2 = function cVals(bcoeffs, rmsecv)
 Crit1 = 'norm(regcoeffs) + cverror'