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

Матричный coefs 5 на 5 содержит коэффициенты основного компонента (один столбец для каждого основного компонента). Матричный 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