Создайте биплот первых трех коэффициентов основных компонентов, наблюдений и наблюдаемых переменных для carsmall набор данных.

Загрузите выборочные данные.

load carsmall

Создайте матрицу, состоящую из переменных Acceleration, Displacement, Horsepower, MPG, и Weight. Удалите строки в матрице, которые имеют отсутствующие значения.

X = [Acceleration Displacement Horsepower MPG Weight];
X = rmmissing(X);

Стандартизируйте X и выполните анализ основного компонента.

Z = zscore(X); % Standardized data
[coefs,score] = pca(Z);

Матрица 5 на 5 coefs содержит коэффициенты основного компонента (по одному столбцу для каждого основного компонента). Матрица score содержит основные счета компонента (наблюдения).

Создайте биплот первых трех коэффициентов главных компонентов. Оси биплота представляют столбцы coefs, и векторы в биплоте представляют строки coefs.

biplot(coefs(:,1:3))

Создайте более подробный биграфик путем маркировки каждой переменной и построения графиков наблюдений в пространстве первых трех основных компонентов.

vbls = {'Accel','Disp','HP','MPG','Wgt'}; % Labels for the variables
biplot(coefs(:,1:3),'Scores',score(:,1:3),'VarLabels',vbls);

Загрузите fisheriris набор данных, стандартизируйте измерения цветов в meas, и выполните анализ основного компонента.

load fisheriris
Z = zscore(meas);
[coefs,scores] = pca(Z);

Создайте рисунок с двумя подграфиками и верните Axes объекты как ax1 и ax2. Создайте биплот в каждом наборе осей путем обращения к соответствующему Axes объект. В верхней подграфике отобразите биплот с помощью первых двух основных компонентов. В нижней подграфике отобразите биплот с использованием третьего и четвертого основных компонентов. Задайте пределы оси X и оси Y путем прохождения соответствующей Axes объекты к xlim и ylim функций. Измените метки оси X и оси Y на нижнем графике путем передачи ax2 на xlabel и ylabel.

figure('Units','normalized','Position',[0.3 0.3 0.3 0.5])
variables = {'SepalLength','SepalWidth','PetalLength','PetalWidth'};
ax1 = subplot(2,1,1); % Top subplot
biplot(ax1,coefs(:,1:2),'Scores',scores(:,1:2),'VarLabels',variables);
xlim(ax1,[-1 1])
ylim(ax1,[-1 1])

ax2 = subplot(2,1,2); % Bottom subplot
biplot(ax2,coefs(:,3:4),'Scores',scores(:,3:4),'VarLabels',variables);
xlim(ax2,[-1 1])
ylim(ax2,[-1 1])
xlabel(ax2,'Component 3')
ylabel(ax2,'Component 4')

Управляйте внешним видом биплота, задавая поддерживаемые имена свойства и значения линий и используя указатели на графические объекты, созданные biplot.

Загрузите выборочные данные.

load carsmall

Создайте матрицу, состоящую из переменных Acceleration, Displacement, и MPG. Удалите строки в матрице, которые имеют отсутствующие значения.

X = [Acceleration Displacement MPG];
X = rmmissing(X);

Стандартизируйте X и выполните анализ основного компонента.

Z = zscore(X); % Standardized data
[coefs,score] = pca(Z);

Матрица 3 на 3 coefs содержит коэффициенты основного компонента (по одному столбцу для каждого основного компонента). Матрица score содержит основные счета компонента (наблюдения).

Создайте биографию наблюдений в пространстве первых двух основных компонентов. Используйте свойства по умолчанию для биплота.

h = biplot(coefs(:,1:2),'Scores',score(:,1:2));

h - вектор указателей для графических объектов. Можно изменять свойства объектов линий, возвращаемых biplot.

Пометьте три переменные для простой идентификации. Задайте круги как символ маркера и синий как цвет линии для всех объектов линии.

vbls = {'Accel','Disp','MPG'}; % Array of variable labels
h1 = biplot(coefs(:,1:2),'Scores',score(:,1:2),...
    'Color','b','Marker','o','VarLabels',vbls);

h1 - вектор указателей для графических объектов. Просмотр первых нескольких элементов h1.

h1(1:10) % First ten object handles

ans = 
  10x1 graphics array:

  Line    (varline)
  Line    (varline)
  Line    (varline)
  Line    (varmarker)
  Line    (varmarker)
  Line    (varmarker)
  Text    (varlabel)
  Text    (varlabel)
  Text    (varlabel)
  Line    (obsmarker)

Указатели на метки переменной (h1(7:9)) являются текстом. Поэтому параметры, заданные для свойств линии, не влияют на эти метки.

Создайте другой биплот наблюдений в пространстве первых двух основных компонентов и пометьте три переменные для легкой идентификации.

h2 = biplot(coefs(:,1:2),'Scores',score(:,1:2),'VarLabels',vbls);

h2 - вектор указателей для графических объектов. Просмотр первых нескольких элементов h2.

h2(1:10) % First ten object handles

ans = 
  10x1 graphics array:

  Line    (varline)
  Line    (varline)
  Line    (varline)
  Line    (varmarker)
  Line    (varmarker)
  Line    (varmarker)
  Text    (varlabel)
  Text    (varlabel)
  Text    (varlabel)
  Line    (obsmarker)

h2 содержит 104 указателя на объекты.

Измените определенные свойства биплота с помощью указателей на графические объекты.

Измените цвет линии переменных (векторов).

for k = 1:3
    h2(k).Color = 'r'; % Specify red as the line color
end

Измените шрифт меток переменных.

for k = 7:9
    h2(k).FontWeight = 'bold';  % Specify bold font
end

Измените цвет маркеров наблюдения.

for k = 10:103
    h2(k).MarkerEdgeColor = 'k';  % Specify black color for the observations
end