mfig('Bi-Plot: Standardized Data');
hbi = biplot(coeff(:,1:3), ...
    'Scores', score(:,1:3), ...
    'VarLabels', Labels, ...
    'ObsLabels', num2str(Statuses), ...
    'markersize', 15 ...
);
title('Bi-Plot: Standardized Data');
xlabel('Bi-Plot: Standardized Data');
xlabel('PCA1');
ylabel('PCA2');

% Manipulate plot colors
for ii = 1:length(hbi)-length(Data)
    set(hbi(ii), 'Color', [0.2 0.2 0.2]);
end
for ii = length(hbi)-length(Data):length(hbi)
    userdata = get(hbi(ii), 'UserData');
    if ~isempty(userdata)
        if Statuses(userdata) == 1
            set(hbi(ii), 'Color', [0 0.75 0]);
        elseif Statuses(userdata) == 0
            set(hbi(ii), 'Color', 'b');
        end
    else
        set(hbi(ii), 'Color', 'k');
    end
end