如何根据每列的频率将np.array(多个基因型)中的字符串替换为浮点数?

我有一个np.array矩阵(1826*5000),其中行是我的样本,列是特征。这意味着每行有一个基因型,个体核苷酸以字符串的形式表示,如下所示:

[['G' 'G' 'G' ... 'T' 'T' 'A'] ['G' 'G' 'G' ... 'A' 'T' 'A'] ['A' 'G' 'A' ... 'A' 'T' 'A'] ... ['G' 'A' 'G' ... 'T' 'T' 'A'] ['G' 'G' 'A' ... 'A' 'T' 'A'] ['G' 'G' 'G' ... 'A' 'T' 'C']]

每列中只出现两种不同的核苷酸。

现在我想用0和2替换这些个体字符串,使得每列中出现频率较高的核苷酸被替换为0,出现频率较低的核苷酸被替换为2。

这意味着在第一列中,“G”应该被替换为0,“A”应该被替换为2,因为“G”出现的频率更高。

最终结果应该如下所示:

[['0' '0' '0' ... '2' '0' '0'] ['0' '0' '0' ... '0' '0' '0'] ['2' '0' '2' ... '0' '0' '0'] ... ['0' '2' '0' ... '2' '0' '0'] ['0' '0' '2' ... '0' '0' '0'] ['0' '0' '0' ... '0' '0' '2']]

有人能告诉我如何做到这一点(使用Sklearn和Numpy函数)吗?


回答:

给定一个数组arr,最简单的解决方法是:

解释:首先,你将数组转换为Pandas数据框。然后,对于每一列,你用“2”替换众数,用“0”替换其他值。最后,你将数据框转换回我们命名为arr1的数组。

Related Posts

L1-L2正则化的不同系数

我想对网络的权重同时应用L1和L2正则化。然而,我找不…

使用scikit-learn的无监督方法将列表分类成不同组别,有没有办法?

我有一系列实例,每个实例都有一份列表,代表它所遵循的不…

f1_score metric in lightgbm

我想使用自定义指标f1_score来训练一个lgb模型…

通过相关系数矩阵进行特征选择

我在测试不同的算法时,如逻辑回归、高斯朴素贝叶斯、随机…

可以将机器学习库用于流式输入和输出吗?

已关闭。此问题需要更加聚焦。目前不接受回答。 想要改进…

在TensorFlow中,queue.dequeue_up_to()方法的用途是什么?

我对这个方法感到非常困惑,特别是当我发现这个令人费解的…

发表回复

您的邮箱地址不会被公开。 必填项已用 * 标注