神技能|一个基因的突变，会不会影响某两个基因表达量的相关性？（代码+短视频）

 今天是生信星球陪你的第632天

   大神一句话，菜鸟跑半年。我不是大神，但我可以缩短你走弯路的半年~

   就像歌儿唱的那样，如果你不知道该往哪儿走，就留在这学点生信好不好~

   这里有豆豆和花花的学习历程，从新手到进阶，生信路上有你有我！

这是一个突变数据和表达量数据关联起来的例子，数据来自TCGA的KIRC。

按照A基因是否发生了突变，将病人分为两组（T，F），将病人表达矩阵的所有基因（或miRNA）两两组合，计算其相关性，如果突变组（T）和未突变组（F）之间的某两个基因（或miRNA）相关性不同，是否意味着A基因的突变与这两个基因的相关关系有故事可讲？

两两组合计算相关性是不是听起来很难。我已经把它写成函数了，可以直接使用：批量计算表达矩阵任意两个基因的相关性。100个基因会约有5000种组合，500个基因约有15万种组合，约运行了30s，只需要等一等即可~基因数量更多时，组合数量是指数式增长的，相应的计算时间也会边长~

代码已打包好，公众号后台回复“mut632”即可获得，如果对你有用，转发、在看，支持一下可好~