+ geom_segment(data = coor, aes(x = x, y = y, xend = xend, yend =yend),+ lty =1, colour ="red")
平常我们利用GEODatasets的表达量数据进行基因表达分析的时候,常常是下载表达矩阵,利用limma包进行分析,如果能找到注释包,那就把探针注释成Entrez_id或者是Symbol_id,如果找不到,那就下载对应GPL提供的注释文件,但是很可能还是不好注释,比如今天拿来做例子的GSE15222对应的GPL2700: 我们看到这个GPL只提供了GB_ACC也就是...