【問(wèn)題】

library(GEOquery)

Gset = getGEO('GSE2034',destdir = ".", GSEMatrix =TRUE)#文件夾中有會(huì)直接讀入

a=Gset[[1]]

dat=exprs(a)

pd=pData(a)

跑出來(lái)↓

Found 1 file(s)

GSE2034_series_matrix.txt.gz

Using locally cached version: ./GSE2034_series_matrix.txt.gz

─ Column specification ─────────────────────────────────

cols(

? .default = col_double(),

? ID_REF = col_character()

)

? Use `spec()` for the full column specifications.

↓?

問(wèn)題

可以看到行名probe 沒(méi)有出來(lái)。。。

【解決方法】加上getGPL = FALSE

Gset = getGEO('GSE2034',destdir = ".", GSEMatrix =TRUE, getGPL = FALSE)

Found 1 file(s)

GSE2034_series_matrix.txt.gz

Using locally cached version: ./GSE2034_series_matrix.txt.gz

─ Column specification ─────────────────────────────────

cols(

? .default = col_double(),

? ID_REF = col_character()

)

? Use `spec()` for the full column specifications.

|=============================================================================| 100% 35 MB

還是會(huì)有提醒，但是probe name已經(jīng)出來(lái)了。至于GPL可以另外下載哦。

問(wèn)題解決

【參考資料】https://www.biostars.org/p/396021/ （Kevin老師yyds