biplot中图例不正确,我有两个目标变量,可以在一个包含两个不同图例的PCA中添加它们吗?
创始人
2024-12-19 03:01:13
0

若想在一个PCA biplot中添加两个目标变量并同时显示两个不同的图例,可以采用如下方法:

from sklearn.decomposition import PCA
import matplotlib.pyplot as plt

# 使用sklearn PCA计算主成分
pca = PCA(n_components=2)
pc = pca.fit_transform(X)

# 绘制双变量坐标系 PCA biplot
plt.scatter(pc[:,0], pc[:,1], alpha=0.5)
plt.ylim(-0.3, 0.3)

# 添加第一个目标变量
for i, target in enumerate(targets):
  indicesToKeep = target == y
  plt.scatter(pc[indicesToKeep, 0], pc[indicesToKeep, 1], alpha=0.8, marker='o', c=colors[i], edgecolor='black', label=target)

# 添加第二个目标变量
for i, target2 in enumerate(targets2):
  indicesToKeep2 = target2 == y2
  plt.scatter(pc[indicesToKeep2, 0], pc[indicesToKeep2, 1], alpha=0.8, marker='+', c=colors2[i], edgecolor='black', label=target2)

# 添加图例
plt.legend(loc='best', framealpha=0.5)

# 为每个特征添加箭头
features = ['sepal length (cm)', 'sepal width (cm)', 'petal length (cm)', 'petal width (cm)']
for i, feature in enumerate(features):
    plt.arrow(0, 0, pca.components_[0,i], pca.components_[1,i], head_width=0.05, head_length=0.1, fc='k', ec='k')
    plt.text(pca.components_[0,i]* 1.15, pca.components_[1,i] * 1.15, feature, color='k', fontsize=12)

plt.show()

其中,targets和targets2分别是两个目标变量的名称列表;y和y2是数据集中分类变量的列表(即两个目标变量对应的分类变量);colors和colors2分别是两个目标变量的标记颜色列表;X是数据矩阵,每一行代表一个样本;注意要将目标变量对应的分类变量和标记颜色一一对应,并在散点图中添加相应的图例。

另外,需要注意的是,在绘制双变量坐标系的biplot时,如果两个主成分不是同一单位的,应当对两个主成分做归一化处理,并且设置纵轴的取值范围。

相关内容

热门资讯

安卓系统怎么连不上carlif... 安卓系统无法连接CarLife的原因及解决方法随着智能手机的普及,CarLife这一车载互联功能为驾...
iwatch怎么连接安卓系统,... 你有没有想过,那款时尚又实用的iWatch,竟然只能和iPhone好上好?别急,今天就来给你揭秘,怎...
oppo手机安卓系统换成苹果系... OPPO手机安卓系统换成苹果系统:现实吗?如何操作?随着智能手机市场的不断发展,用户对于手机系统的需...
安卓平板改windows 系统... 你有没有想过,你的安卓平板电脑是不是也能变身成Windows系统的超级英雄呢?想象在同一个设备上,你...
iphone系统与安卓系统更新... 最近是不是你也遇到了这样的烦恼?手机更新系统总是失败,急得你团团转。别急,今天就来给你揭秘为什么iP...
安卓系统上滑按键,便捷生活与高... 你有没有发现,现在手机屏幕越来越大,操作起来却越来越方便了呢?这都得归功于安卓系统上的那些神奇的上滑...
安卓系统连接耳机模式,蓝牙、有... 亲爱的手机控们,你们有没有遇到过这种情况:手机突然变成了“耳机模式”,明明耳机没插,声音却只从耳机孔...
希沃系统怎么装安卓系统,解锁更... 亲爱的读者们,你是否也像我一样,对希沃一体机上的安卓系统充满了好奇呢?想象在教室里,你的希沃一体机不...
安装了Anaconda之后找不... 在安装Anaconda后,如果找不到Jupyter Notebook,可以尝试以下解决方法:检查环境...
安卓换鸿蒙系统会卡吗,体验流畅... 最近手机圈可是热闹非凡呢!不少安卓用户都在议论纷纷,说鸿蒙系统要来啦!那么,安卓手机换上鸿蒙系统后,...