Создайте биплот из первых трех коэффициентов главных компонентов, наблюдений и наблюдаемых переменных для 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