NMR Batch Correction Example Script - Revision history

Mjeppesen at 05:31, 20 January 2022

2022-01-20T05:31:21Z

← Older revision		Revision as of 05:31, 20 January 2022
Line 1:		Line 1:
	== PLS Fit Model ==		== PLS Fit Model ==
	X = ~~whole dataset~~ (first & second);clsD = ~~class info~~ (~~dummy matrix~~).		first_value = your_choice1 %%You choose here
			second_value = your_choice2 %%You choose here too
			X = whole_dataset ([first_value:second_value],:); %%This selects a subset of data from your first choice to your second choice
			clsD = class_info (dummy_matrix);
	result = pls(X, clsD, @suv, [7,1], 1);		result = pls(X, clsD, @suv, [7,1], 1);
	mn = result.mean.X;sd = result.scale.X;		mn = result.mean.X;
			sd = result.scale.X;
	Xscaled = (X-(repmat(mn,[size(X,1),1])))./(repmat(sd,[size(X,1),1]));%size(Xscaled)		Xscaled = (X-(repmat(mn,[size(X,1),1])))./(repmat(sd,[size(X,1),1]));%size(Xscaled)
	XPLSfit = result.T*result.P';		XPLSfit = result.T*result.P';
Line 21:		Line 25:
	Note: mean can be replaced with median in a median of median model.		Note: mean can be replaced with median in a median of median model.

	~~[[category:Protocols]]~~
	[[category:Data_Processing_and_Analysis]]		[[category:Data_Processing_and_Analysis]]

Mjeppesen: LearningSorting

2021-09-11T04:52:18Z

LearningSorting

← Older revision		Revision as of 04:52, 11 September 2021
Line 21:		Line 21:
	Note: mean can be replaced with median in a median of median model.		Note: mean can be replaced with median in a median of median model.

	[[category:~~Protocols\|~~Protocols]]		[[category:Protocols]]
	[[category:Data_Processing_and_Analysis]]		[[category:Data_Processing_and_Analysis]]

Mjeppesen: Added to Category

2021-09-11T03:26:14Z

Added to Category

← Older revision		Revision as of 03:26, 11 September 2021
Line 20:		Line 20:

	Note: mean can be replaced with median in a median of median model.		Note: mean can be replaced with median in a median of median model.

			[[category:Protocols\|Protocols]]
			[[category:Data_Processing_and_Analysis]]

Acrook: Created page with "== PLS Fit Model == X = whole dataset (first & second);clsD = class info (dummy matrix). result = pls(X, clsD, @suv, [7,1], 1); mn = result.mean.X;sd = result.scale...."

2021-01-07T17:09:22Z

Created page with "== PLS Fit Model == X = whole dataset (first & second);clsD = class info (dummy matrix). result = pls(X, clsD, @suv, [7,1], 1); mn = result.mean.X;sd = result.scale...."

New page

== PLS Fit Model ==
X = whole dataset (first & second);clsD = class info (dummy matrix).
result = pls(X, clsD, @suv, [7,1], 1);
mn = result.mean.X;sd = result.scale.X;
Xscaled = (X-(repmat(mn,[size(X,1),1])))./(repmat(sd,[size(X,1),1]));%size(Xscaled)
XPLSfit = result.T*result.P';
XscaledNew = Xscaled - XPLSfit;
XNew = (XscaledNew.*(repmat(sd,[size(XscaledNew,1),1]))) + (repmat(mn,[size(XscaledNew,1),1]));

== Mean of Mean Model ==

Xall = [Xtrain;Xtest];
refX = median(median(Xall));
mnXtrain = mean(mean(Xtrain));
newXtrain = Xtrain/repmat(mnXtrain, [size(Xtrain, 1), size(Xtrain, 2)]);
newXtrain = newXtrain*repmat(refX, [size(Xtrain, 1), size(Xtrain, 2)]);
mnXtest = mean(mean(Xtest));
newXtest = Xtest/repmat(mnXtest, [size(Xtest, 1), size(Xtest, 2)]);
newXtest = newXtest*repmat(refX, [size(Xtest, 1), size(Xtest, 2)]);

Note: mean can be replaced with median in a median of median model.