代码如下:
library(pacman)
library(SummarizedExperiment)
library(TCGAbiolinks)
library(tidyverse)
library(magrittr)
library(data.table)
library(biomaRt)
projectPath = "D:/Project/sxkong_200415"
dataPath = paste(projectPath, "Data", sep = "/")
if(!dir.exists(dataPath)) dir.create(dataPath)
setwd(dataPath)
TCGAbiolinks:::getProjectSummary("TCGA-LIHC")
LIHC_query <- GDCquery(project = "TCGA-LIHC",
legacy = FALSE,
experimental.strategy = "RNA-Seq",
data.category = "Transcriptome Profiling",
data.type = "Gene Expression Quantification",
workflow.type = "HTSeq - FPKM-UQ")
GDCdownload(LIHC_query,files.per.chunk = 50)
dataAssy = GDCprepare(LIHC_query, summarizedExperiment = F)
R版本4.0.5 下载完数据运行 : dataAssy = GDCprepare(LIHC_query, summarizedExperiment = F)时,
出现以下错误
> dataAssy = GDCprepare(LIHC_query, summarizedExperiment = F,directory = "GDCdata")
| | 0% 错误: 'C:\Users\神的孩~1\AppData\Local\Temp\RtmpMl9OQE\file8a46f293134' does not exist.
Error in if (value == n) { : 参数长度为零
> dataAssy = GDCprepare(LIHC_query, summarizedExperiment = F,directory = "GDCdata")
| | 0% 错误: 'C:\Users\神的孩~1\AppData\Local\Temp\RtmpMl9OQE\file8a45b9d58fd' does not exist.
Error in if (value == n) { : 参数长度为零
> dataAssy = GDCprepare(LIHC_query, summarizedExperiment = F)
| | 0% 错误: 'C:\Users\神的孩~1\AppData\Local\Temp\RtmpMl9OQE\file8a4798e51d6' does not exist.
Error in if (value == n) { : 参数长度为零
求大神指点,小白一个,鼓捣两天了,实在是没辙了T T
可能是下载的数不完整,或者数据存放路径有问题。建议先将dataPath这个部分代码注释掉,这样来写:
BiocManager::install("TCGAbiolinks")
library(TCGAbiolinks)
?GDCquery
x<-getGDCprojects()$project_id
TCGAbiolinks:::getProjectSummary("TCGA-LIHC")
query<-GDCquery(project = "TCGA-LIHC",
legacy=F,
experimental.strategy ="RNA-Seq",
data.category = "Transcriptome Profiling",
data.type="Gene Expression Quantification",
workflow.type="HTSeq - FPKM-UQ")
GDCdownload(query,files.per.chunk = 50)
dataAssy = GDCprepare(query, summarizedExperiment = F)
|=========================================================================|100% Completed after 47 s
别急,我看下这个问题,过会答复你。
您好,我是有问必答小助手,你的问题已经有小伙伴为您解答了问题,您看下是否解决了您的问题,可以追评进行沟通哦~
如果有您比较满意的答案 / 帮您提供解决思路的答案,可以点击【采纳】按钮,给回答的小伙伴一些鼓励哦~~
ps:问答VIP仅需29元,即可享受5次/月 有问必答服务,了解详情>>>https://vip.csdn.net/askvip?utm_source=1146287632