最初の 3 つの主成分の係数、観測値、およびcarsmallデータセットの観測された変数からバイプロットを作成します。

標本データを読み込みます。

loadcarsmall

変数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には、主成分の係数が格納されます (各主成分に対して 1 列)。行列scoreには、主成分のスコア (観測値) が格納されます。

最初の 3 つの主成分の係数からバイプロットを作成します。バイプロットの座標軸はcoefsの列を、バイプロットのベクトルはcoefsの行を表します。

biplot(coefs(:,1:3))

各変数にラベルを付け、最初の 3 つの主成分からなる空間に観測値をプロットすることにより、より詳細なバイプロットを作成します。

vbls = {'Accel','Disp','HP',“英里”,'Wgt'};% Labels for the variablesbiplot(coefs(:,1:3),'Scores',score(:,1:3),'VarLabels',vbls);

fisheririsデータセットを読み込み、花の測定値をmeasで標準化し、主成分分析を実行します。

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

2 つのサブプロットを使用して Figure を作成し、ax1およびax2としてAxesオブジェクトを返します。対応するAxesオブジェクトを参照して、各座標軸のセットでバイプロットを作成します。上のサブプロットに、最初の 2 つの主成分を使用してバイプロットを表示します。下のサブプロットに、3 番目と 4 番目の主成分を使用してバイプロットを表示します。対応するAxesオブジェクトを関数xlimとylimに渡して、X 軸と Y 軸の範囲を指定します。ax2をxlabelとylabelに渡して、下のプロットの X 軸と Y 軸のラベルを変更します。

figure('Units','normalized','Position',[0.3 0.3 0.3 0.5]) variables = {'SepalLength','SepalWidth','PetalLength','PetalWidth'}; ax1 = subplot(2,1,1);% Top subplotbiplot(ax1,coefs(:,1:2),'Scores',scores(:,1:2),'VarLabels',variables); xlim(ax1,[-1 1]) ylim(ax1,[-1 1]) ax2 = subplot(2,1,2);% Bottom subplotbiplot(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')

サポートされている line プロパティの名前と値を指定し、biplotによって作成されたグラフィックス オブジェクトのハンドルを使用することにより、バイプロットの外観を制御します。

標本データを読み込みます。

loadcarsmall

変数Acceleration、DisplacementおよびMPGから構成される行列を作成します。この行列から、欠損値が含まれている行を削除します。

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

Xを標準化し、主成分分析を実行します。

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

3 行 3 列の行列coefsには、主成分の係数が格納されます (各主成分に対して 1 列)。行列scoreには、主成分のスコア (観測値) が格納されます。

最初の 2 つの主成分からなる平面に、観測値のバイプロットを作成します。バイプロットには、既定のプロパティを使用します。

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

hは、グラフィックス オブジェクトのハンドルのベクトルです。biplotによって返された line オブジェクトのプロパティを変更できます。

識別を容易にするため、3 つの変数にラベルを付けます。すべての line オブジェクトに対して、マーカー記号として円を、ラインの色として青を指定します。

vbls = {'Accel','Disp',“英里”};% Array of variable labelsh1 = 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)) はテキストです。したがって、line プロパティについて指定した設定はこれらのラベルに影響を与えません。

最初の 2 つの主成分からなる平面に別の観測値のバイプロットを作成し、識別を容易にするため 3 つの変数にラベルを付けます。

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 個のオブジェクト ハンドルが含まれています。

グラフィックスオブジェクトのハンドルを使用して,バイプロットの特定のプロパティを変更します。

変数 (ベクトル) のラインの色を変更します。

fork = 1:3 h2(k).Color ='r';% Specify red as the line colorend

変数ラベルのフォントを変更します。

fork = 7:9 h2(k).FontWeight ='bold';% Specify bold fontend

観測値マーカーの色を変更します。

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