我试图计算由数据集中每个基因的表达和甲基化数据之间的第三个变量调节的互信息。
我的数据包含一个包含两个矩阵的列表，对于第一个元素，其内容是在35名患有肺腺癌的患者中测量的10000个基因的RPKM表达值。该列表的第二个元素是由35个患者中相同的10000个基因的β甲基化值组成的另一个矩阵。两个矩阵都有列中的基因和行中的基因。
因此，我想计算两种基质中所有成对基因的甲基化和表达值之间的相互信息。但是，与此同时，这两个变量受甲基化阈值的调节，其值在0和1之间振荡。我已经看到了condinformation（）函数，但我不能将这个函数应用于每个成对甲基化和每个基因的表达信息。
非常感谢你！