F改成T试试:
if(F){
library(GEOquery)
eSet <- getGEO('GSE112093', destdir=".",
AnnotGPL = F,
getGPL = F)
save(eSet,file='GSE112093_eSet.Rdata')
}
load('GSE112093_eSet.Rdata')
b = eSet[[1]]
raw_exprSet=exprs(b)
phe=pData(b)
library(stringr)
group_list= str_split(as.character(phe$title),' ',simplify = T)[,1]
save(raw_exprSet,group_list,
file='GSE112093_raw_exprSet.Rdata')
这是生信大佬的代码 我只改了GSE数据集,但是就提示下载不了,但是原代码上的数据集是可以下载的。错误提示如下:Error in readChar(con, 5L, useBytes = TRUE) : 无法打开链结 此外: Warning message: In readChar(con, 5L, useBytes = TRUE) : 无法打开压缩文件'GSE112093_eSet.Rdata',可能是因为'No such file or directory'
恳请大佬解答!