问题:矩阵A存储10个基因和20个样本的基因表达值,矩阵B存储着20个样本的分类情况(包括三类),B矩阵第一列为样本名,第二列为样本类型。想根据样本的类别,计算10个基因的表达平均值。
思路:
将A矩阵转置后,确保A的行名(样本名)和B矩阵的行名一致,然后将B矩阵的Type添加到A矩阵中,然后用aggregate函数:
gene_means <- aggregate(. ~ Type, data = combined_data, mean)
# 打印结果
print(gene_means)