plotnirdesign.m

PCA and PLS aims:to get some insight into the bilinear factor models Principal Component Analysis

function plotNIRDesign(T);
% 020228 FvdB

pure = [1 18 22 1];
purein = [2 4 13 21 19 9 2];
pureinin = [6 7 12 16 15 10 6];
mix = [1 3 12 19 17 8 1];
mixin = [5 6 11 15 14 9 5];

if (size(T,2) ~= 2) & (size(T,2) ~= 3)
    error('ERROR: incorrect number of columns in "T"');
end
if size(T,2) == 2
    if size(T,1) == 22
        plot(T(:,1),T(:,2),'o',T(pure,1),T(pure,2),'g',T(purein,1),T(purein,2),'g',T(pureinin,1),T(pureinin,2),'g');
        text(T([1 18 22],1),T([1 18 22],2),['E' 'W' 'I']');
    elseif size(T,1) == 19
        plot(T(:,1),T(:,2),'o',T(mix,1),T(mix,2),'g',T(mixin,1),T(mixin,2),'g');
    else
        error('ERROR: Incorect matrix "T"')
    end
else
    if size(T,1) == 22
        plot3(T(:,1),T(:,2),T(:,3),'o',T(pure,1),T(pure,2),T(pure,3),'g',T(purein,1),T(purein,2),T(purein,3),'g',T(pureinin,1),T(pureinin,2),T(pureinin,3),'g');
        text(T([1 18 22],1),T([1 18 22],2),T([1 18 22],3),['E' 'W' 'I']');
    elseif size(T,1) == 19
        plot3(T(:,1),T(:,2),T(:,3),'o',T(mix,1),T(mix,2),T(mix,3),'g',T(mixin,1),T(mixin,2),T(mixin,3),'g');
    else
        error('ERROR: Incorect matrix "T"')
    end
end
grid