雖說簡(jiǎn)單，但對(duì)R語(yǔ)言不熟悉的人來說還是需要指點(diǎn)一下，才會(huì)用。首先是GPL下載的文件，以GPL96為例。
提取最感興趣的3列。
好的，接下來需要得到每個(gè)probe 對(duì)應(yīng)單個(gè)symbol name的數(shù)據(jù)框。比如，10073_s_at ——DDR1；10073_s_at ——MIR4640；1053_at——RFC2……
上代碼：
symbola=data.frame(probe=selectAnno$probe,symbol=selectAnno$symbol,
stringsAsFactors = FALSE)#22283
symbola=symbola[which(symbola$symbol!= ""),] #21225 #刪去symbol為空的
symbolb=strsplit(symbola$symbol,split = " /// ") # 分割字符串
names(symbolb) <- symbola$probe #這一步很關(guān)鍵
symbolc <- ldply(symbolb,data.frame) #23749 #list轉(zhuǎn)換為data.frame
colnames(symbolc) <- c("probe","symbol")
save(symbolc,file = "GPL96_probe_symbol.Rdata")
symbola
symbolb
symbolc
大功告成！
