朋友问到以下问题:
假设有100个人,分别检查了10万个miRNA的表达情况和1000个基因的突变情况,然后想基于这100个人的上述检测结果,分别分析每个miRNA与每个突变基因之间两两的相关性,并将所有的结果输出到一个新的数据库,应该如何用R来写程序?
看上去是一件非常复杂的事情,自己并非生物信息学专业的,所以不清楚Bioconductor等软件是否有更成熟快速的方法来完成这个程序。
自己简单写了一个示例程序,并验证结果没有问题,但是如果有学生物信息学的老师熟悉采用更加快捷方便的程序来完成该过程,非常欢迎指点。