3、 比较1和2中的基因组位置关系，并找出来未被基因组注释的EST。
这时候应该还剩下上千条，其中绝大部分都是蛋白编码基因的反义RNA，当然一小部分是基因间区的非编码RNA。
4.、将affimetrix的probe sequence用blast比对到这些非编码RNA中。建议先找找关心的事情是否有tiling-array data。
另外要注意非编码RNA和它可能有关的基因一起分析，比如附近基因，或者互补的基因。
做芯片一般不去找表达或非表达，而往往是在找差异表达。
比如某非编码RNA和它附近的基因共同在肿瘤细胞中特异表达。
seq函数是产生序列用的他的用法是seq(from,to,by)或者是seq（下界，by=，length=） 下面是用r运行的结果 seq(2,6,2) [1] 2 4 6 seq(10,by=2,length=5) [1] 10 12 14 16 18看注释前的字母。R语言初学指南可在脚本中加入注释。在脚本中，任何以“#”（sharp/numbersymbol）开头的命令行都会被R忽略。
同样，若“#”出现在某行的中间，则该行中“#”后面的语句都会被忽略。可利用这一特性对脚本添加注释，以便用户或他人日后查阅。
例如，作者每次查看前一天编写的脚本时，都要重新梳理并回忆每条脚本语句的作用。
