matlab - 線形判別分析で得られた投影からデータを再構築する方法

Question

LDA で得られた投影からデータを再構築しようとしています。アイデアは、縮小された LDA 係数のセットから得られた再構成エラーを評価することです。次の matlab コードでは、射影されたデータ p と固有ベクトル LTrans を使用して再構成を取得する方法が問題になります。

L = fitcdiscr(data,class);
[LTrans,Lambda] = eig(L.BetweenSigma,L.Sigma,'chol');
[Lambda,sorted] = sort(diag(Lambda),'descend'); % sort by eigenvalues
LTrans = LTrans(:,sorted);
xc = L.X;
mu = mean(xc);
Xm = bsxfun(@minus, xc, mu);
for i_fact=1:size(L.Sigma,2)
  z = Xm*LTrans(:,1:i_fact);
  p = z*LTrans(:,1:i_fact)';
  p = bsxfun(@plus, p, mu);
end