包装limma是最好的 - 甚至是最好的 - 用于阅读微阵列 - 最近的RNA-seq - 实验结果和其他过程如标准化，拟合线性模型，寻找顶部差异表达基因，制作情节等等。
我在read.ilmn()函数中有一个小问题，用于读取illumina结果：假设有一些注释列 - 即ENTREZ_GENE_ID，REFSEQ_ID等。除了表达式列（即PROBE_ID，SYMBOL，AVG_Signal和Detection.Pval）之外，还需要读取这些列以供以后分析。
我试图将这些列的名称作为annotation的{{1}}和other.columns参数传递，通知此命令也读取这些列，导致{{1}失败}。最后，我必须通过read.ilmn()命令单独阅读这些列，并将其与read.ilmn()结果手动合并。
有没有更好的方法来执行任务？
答案 0 :(得分：0)
我终于通过手动阅读注释列并将其连接到read.ilmn()的列表结果来解决问题。