表題通りです。とりあえず、sklearn の load_boston の結果を表示する様にしてあります。
from sklearn.datasets import load_boston
X = load_boston(return_X_y=True)[0]
from sklearn import decomposition
pca=decomposition.PCA()
pca.fit(X)
y = np.add.accumulate(pca.explained_variance_ratio_)
x = np.arange(y.size)+1
plt.plot(x, y, ls='--', marker='o')
plt.ylabel("cumulative contribution ratio")
plt.xlabel("number of principal components")
plt.show()
